! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! Sparse Stoichiometric Data Structures File ! ! Generated by KPP-2.2 symbolic chemistry Kinetics PreProcessor ! (http://www.cs.vt.edu/~asandu/Software/KPP) ! KPP is distributed under GPL, the general public licence ! (http://www.gnu.org/copyleft/gpl.html) ! (C) 1995-1997, V. Damian & A. Sandu, CGRER, Univ. Iowa ! (C) 1997-2005, A. Sandu, Michigan Tech, Virginia Tech ! With important contributions from: ! M. Damian, Villanova University, USA ! R. Sander, Max-Planck Institute for Chemistry, Mainz, Germany ! ! File : gckpp_StoichiomSP.f90 ! Time : Tue Mar 16 10:29:04 2010 ! Working directory : /mnt/lstr04/srv/home/c/ccarouge/KPP/geoschem_kppfiles/v8-03-01/Paulot_scheme ! Equation file : gckpp.kpp ! Output root filename : gckpp ! ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODULE gckpp_StoichiomSP USE gckpp_Precision PUBLIC SAVE ! Row-compressed sparse data for the Jacobian of reaction products JVRP INTEGER, PARAMETER, DIMENSION(360) :: CROW_JVRP_0 = (/ & 1, 3, 5, 7, 9, 11, 12, 13, 15, 17, 19, 20, & 21, 23, 24, 26, 28, 29, 30, 32, 34, 36, 38, 40, & 42, 44, 46, 48, 50, 51, 53, 55, 57, 59, 61, 62, & 64, 65, 67, 69, 71, 73, 75, 76, 78, 80, 82, 84, & 86, 88, 90, 92, 94, 96, 98,100,102,104,106,108, & 110,112,114,116,117,119,121,123,125,127,129,130, & 132,133,135,137,139,141,143,145,147,149,151,153, & 155,157,159,161,163,165,167,169,171,173,175,177, & 179,181,183,185,186,188,190,192,194,196,198,200, & 201,203,204,205,207,209,211,213,215,217,219,221, & 223,225,227,229,231,232,234,236,238,240,242,244, & 246,248,250,252,254,256,258,260,262,264,266,268, & 270,272,274,276,278,280,282,284,286,288,290,292, & 294,296,298,300,302,304,306,308,310,312,314,316, & 318,320,322,324,326,328,330,332,334,335,337,339, & 341,343,345,347,348,350,352,354,356,358,360,362, & 364,366,368,370,372,374,376,378,380,382,384,386, & 388,390,392,394,396,398,400,402,404,406,408,409, & 409,409,409,409,409,409,409,409,409,409,409,409, & 409,409,409,409,409,409,410,411,412,413,414,415, & 416,417,418,419,420,421,422,423,424,425,427,429, & 430,431,432,433,434,435,436,437,439,441,443,445, & 447,449,451,453,455,457,459,461,463,465,467,469, & 471,473,475,477,479,481,482,484,486,488,490,492, & 494,496,498,500,502,504,506,508,510,512,514,516, & 518,520,522,524,526,528,529,530,531,532,533,534, & 535,536,537,538,539,540,541,542,543,544,545,546, & 547,548,549,550,551,552,553,554,555,556,557,558, & 559,560,561,562,563,564,565,566,567,568,569,570, & 571,572,573,574,575,576,577,578,579,580,581,582 /) INTEGER, PARAMETER, DIMENSION(2) :: CROW_JVRP_1 = (/ & 583,584 /) INTEGER, PARAMETER, DIMENSION(362) :: CROW_JVRP = (/& CROW_JVRP_0, CROW_JVRP_1 /) INTEGER, PARAMETER, DIMENSION(360) :: ICOL_JVRP_0 = (/ & 134,135,132,135,135,136,128,135,133,135,132,132, & 132,136, 46,132,134,136,136,132, 91,132,132,133, & 134,133,136,133,133, 62,132, 62,132, 63,132, 63, & 132,114,132,128,132, 99,132,132,134, 53,132,128, & 136, 59, 59,132,129,136,129,134,129,132,128,129, & 55, 72,132,132,128,129,114,129,115,132,115,129, & 128,130, 49,130,134, 52,132,112,134, 44,132, 44, & 132, 94,134,109,134, 51,132,126,134,126,134,116, & 134,108,134,106,134,123,134,101,134,105,134, 97, & 134, 98,134,117,134,117,129,117,110,134, 51,129, & 132,137, 89,132,127,132,113,128, 45,122,128, 79, & 113,134,122,134,127,129, 85,132, 85,132, 94,133, & 109,133,126,136,116,136,108,136,106,136,123,136, & 101,136,105,136, 97,136, 98,136,117,136,110,136, & 131,132,112,133,129,131,126,133,116,133,108,133, & 106,133,123,133,123,111,133,101,133,105,133, 97, & 133, 98,133,117,133,110,133,132, 88,132,112,112, & 112,136, 94,136,109,136,130,136,113,136,122,136, & 86,132, 86,135, 79,132, 95,134, 95,136, 79,135, & 75,128, 65, 75,134, 75,136,100,132,100,132, 86, & 129, 90,132,118,132, 90,129,118,129, 81,132,120, & 132,121,132,119,132,119,132, 94,130,109,130, 94, & 130,109,130, 81,135,120,135,121,135, 81,129,121, & 129,121,129,113,133,122,133,113,133,122,133, 71, & 132, 71,132, 67,132, 61,132, 56,132, 57,132, 66, & 132, 58,132, 69,132, 74,132, 74,132, 47,132, 84, & 136, 84,134, 76,132, 73,132, 70,132, 70,132, 64, & 132, 64,132,102,130,102,130,102,133,102,102,136, & 102,134,102,134, 54,132, 52,129,124,135,130,130, & 133,130,133,126,130,108,130,106,130,123,130,111 /) INTEGER, PARAMETER, DIMENSION(223) :: ICOL_JVRP_1 = (/ & 130,101,130,105,130, 98,130,116,130, 97,130,117, & 130,110,130,126,130,108,130,106,130,123,130,111, & 130,101,130,105,130,116,130, 97,130,117,130,110, & 130, 98,130,112,130,112,130,113,130,122,130,129, & 128, 50,135, 49, 65, 99,114, 55, 46, 79, 45,137, & 47, 74,125, 82, 82,135,125,135, 96,104, 93, 72, & 89,128,129, 55, 60,132, 60,132, 60,129, 48,132, & 124,132,111,134,111,134,111,136,125,132, 82,132, & 103,134,103,134,103,136,107,134,107,134,107,136, & 80,132, 93,132,104,132, 92,134, 92,136, 92,128, & 50, 77,132, 83,134, 83,134, 83,136, 78,132, 87, & 134, 87,134, 87,136, 78,135, 68,132, 96,132, 18, & 132, 22,132, 26,132, 17,136, 17,134, 21,136, 21, & 134, 25,136, 25,134, 30,132, 29,136, 29,134,135, & 128, 46, 62,114,114, 99, 53, 59,129,129, 55, 55, & 59,115,115, 49,127, 85, 85,131,100, 90, 90,118, & 118,120,120,120,121,121,119, 71, 67, 61, 56, 57, & 66, 69, 58, 74, 76, 80, 73, 70, 64,137, 54,104, & 93, 82,125, 96, 63, 65, 64 /) INTEGER, PARAMETER, DIMENSION(583) :: ICOL_JVRP = (/& ICOL_JVRP_0, ICOL_JVRP_1 /) INTEGER, PARAMETER, DIMENSION(360) :: IROW_JVRP_0 = (/ & 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 7, & 8, 8, 9, 9, 10, 10, 11, 12, 13, 13, 14, 15, & 15, 16, 16, 17, 18, 19, 19, 20, 20, 21, 21, 22, & 22, 23, 23, 24, 24, 25, 25, 26, 26, 27, 27, 28, & 28, 29, 30, 30, 31, 31, 32, 32, 33, 33, 34, 34, & 35, 36, 36, 37, 38, 38, 39, 39, 40, 40, 41, 41, & 42, 42, 43, 44, 44, 45, 45, 46, 46, 47, 47, 48, & 48, 49, 49, 50, 50, 51, 51, 52, 52, 53, 53, 54, & 54, 55, 55, 56, 56, 57, 57, 58, 58, 59, 59, 60, & 60, 61, 61, 62, 62, 63, 63, 64, 65, 65, 66, 66, & 67, 67, 68, 68, 69, 69, 70, 70, 71, 72, 72, 73, & 74, 74, 75, 75, 76, 76, 77, 77, 78, 78, 79, 79, & 80, 80, 81, 81, 82, 82, 83, 83, 84, 84, 85, 85, & 86, 86, 87, 87, 88, 88, 89, 89, 90, 90, 91, 91, & 92, 92, 93, 93, 94, 94, 95, 95, 96, 96, 97, 97, & 98, 98, 99, 99,100,101,101,102,102,103,103,104, & 104,105,105,106,106,107,107,108,109,109,110,111, & 112,112,113,113,114,114,115,115,116,116,117,117, & 118,118,119,119,120,120,121,121,122,122,123,123, & 124,124,125,126,126,127,127,128,128,129,129,130, & 130,131,131,132,132,133,133,134,134,135,135,136, & 136,137,137,138,138,139,139,140,140,141,141,142, & 142,143,143,144,144,145,145,146,146,147,147,148, & 148,149,149,150,150,151,151,152,152,153,153,154, & 154,155,155,156,156,157,157,158,158,159,159,160, & 160,161,161,162,162,163,163,164,164,165,165,166, & 166,167,167,168,168,169,169,170,170,171,171,172, & 172,173,173,174,174,175,175,176,176,177,178,178, & 179,179,180,180,181,181,182,182,183,183,184,185, & 185,186,186,187,187,188,188,189,189,190,190,191 /) INTEGER, PARAMETER, DIMENSION(223) :: IROW_JVRP_1 = (/ & 191,192,192,193,193,194,194,195,195,196,196,197, & 197,198,198,199,199,200,200,201,201,202,202,203, & 203,204,204,205,205,206,206,207,207,208,208,209, & 209,210,210,211,211,212,212,213,213,214,214,215, & 234,235,236,237,238,239,240,241,242,243,244,245, & 246,247,248,249,250,250,251,251,252,253,254,255, & 256,257,258,259,260,260,261,261,262,262,263,263, & 264,264,265,265,266,266,267,267,268,268,269,269, & 270,270,271,271,272,272,273,273,274,274,275,275, & 276,276,277,277,278,278,279,279,280,280,281,281, & 282,283,283,284,284,285,285,286,286,287,287,288, & 288,289,289,290,290,291,291,292,292,293,293,294, & 294,295,295,296,296,297,297,298,298,299,299,300, & 300,301,301,302,302,303,303,304,304,305,305,306, & 307,308,309,310,311,312,313,314,315,316,317,318, & 319,320,321,322,323,324,325,326,327,328,329,330, & 331,332,333,334,335,336,337,338,339,340,341,342, & 343,344,345,346,347,348,349,350,351,352,353,354, & 355,356,357,358,359,360,360 /) INTEGER, PARAMETER, DIMENSION(583) :: IROW_JVRP = (/& IROW_JVRP_0, IROW_JVRP_1 /) ! Stoichiometric Matrix in Compressed Column Sparse Format INTEGER, PARAMETER, DIMENSION(360) :: CCOL_STOICM_0 = (/ & 1, 4, 7, 10, 13, 17, 19, 21, 23, 26, 30, 32, & 34, 38, 40, 45, 48, 50, 53, 56, 58, 61, 63, 67, & 70, 73, 76, 79, 82, 85, 88, 92, 95, 99,102,105, & 109,112,114,119,125,129,132,135,140,143,148,151, & 154,159,165,168,180,183,190,196,202,213,222,231, & 239,244,250,256,261,266,270,273,277,280,283,286, & 289,292,296,304,308,311,314,320,328,331,334,340, & 345,352,362,370,378,381,386,389,392,397,401,413, & 422,429,437,445,452,462,471,479,488,494,504,512, & 515,519,522,524,527,530,533,540,546,556,559,567, & 570,580,585,590,593,596,601,606,613,619,622,626, & 630,635,640,644,647,651,655,661,666,672,676,681, & 691,699,707,710,713,717,722,727,730,733,736,740, & 744,748,752,756,760,763,767,771,773,776,787,800, & 804,808,811,814,817,820,826,830,838,842,847,853, & 856,859,863,872,874,877,881,892,895,898,908,920, & 927,935,940,948,958,967,973,977,981,985,989,993, & 997,1001,1006,1011,1016,1021,1025,1030,1034,1038,1041,1043, & 1044,1045,1046,1047,1048,1049,1050,1051,1052,1053,1054,1055, & 1056,1057,1058,1059,1060,1061,1064,1067,1070,1073,1076,1079, & 1082,1085,1088,1091,1094,1097,1100,1103,1106,1109,1115,1122, & 1125,1128,1131,1134,1137,1140,1142,1144,1149,1154,1160,1164, & 1167,1179,1182,1188,1191,1194,1207,1210,1224,1233,1236,1246, & 1251,1258,1261,1271,1282,1285,1288,1294,1302,1305,1311,1314, & 1320,1323,1329,1335,1340,1344,1346,1348,1350,1352,1354,1356, & 1358,1360,1362,1364,1366,1368,1370,1373,1375,1379,1382,1384, & 1387,1390,1393,1396,1398,1401,1405,1408,1412,1414,1419,1423, & 1426,1429,1434,1438,1441,1444,1448,1451,1454,1459,1462,1465, & 1470,1474,1479,1484,1488,1492,1496,1500,1505,1509,1516,1523, & 1528,1535,1541,1545,1556,1559,1566,1570,1574,1578,1583,1587 /) INTEGER, PARAMETER, DIMENSION(2) :: CCOL_STOICM_1 = (/ & 1592,1596 /) INTEGER, PARAMETER, DIMENSION(362) :: CCOL_STOICM = (/& CCOL_STOICM_0, CCOL_STOICM_1 /) INTEGER, PARAMETER, DIMENSION(360) :: IROW_STOICM_0 = (/ & 128,134,135,132,135,136,132,135,136,128,129,135, & 114,133,135,136,132,135, 46,132,132,136, 46,132, & 136,128,132,134,136, 46,136,132,136, 42, 91,132, & 136,132,133,114,128,133,134,136, 62,133,136,114, & 133,114,133,136, 62,132,133, 62,114, 63,108,132, & 63,118, 91,114,132,136, 99,128,132, 99,129,132, & 53,132,134, 53,128,132, 59,128,136, 59,128,136, & 59,128,132,128,129,132,136,128,129,134,128,129, & 132,136, 55,128,129, 55,128,129, 42, 72,132,136, & 114,132,136,129,134, 91, 99,114,129,136, 91,114, & 115,130,132,136, 99,115,129,130, 49,128,130, 49, & 128,130, 42,128,130,133,134, 52,112,132,112,115, & 128,134,136, 44, 98,132, 44, 94,132, 94,127,128, & 134,136,109,114,115,128,134,136, 51,126,132, 85, & 94, 98,112,115,126,127,128,131,133,134,136,126, & 134,137,114,115,116,126,127,128,134,108,114,128, & 130,134,137,106,115,128,130,134,137, 82, 83,114, & 120,121,123,124,125,128,134,136, 93,100,101,114, & 118,128,130,134,136, 91,104,105,114,118,119,128, & 134,136, 91, 96, 97,114,118,128,132,134, 85, 98, & 128,134,136,117,124,125,128,134,136,117,124,125, & 128,129,136,117,124,125,128,136,110,114,115,128, & 134, 51, 99,126,129,116,132,137, 42, 89,132,133, & 113,127,132, 45,113,128, 45,113,128, 79,122,128, & 79,122,128,112,113,128,134, 42, 91,114,122,128, & 130,133,134, 99,113,127,129, 85,108,132, 85,108, & 132, 88, 94,114,127,133,136, 88,109,114,115,119, & 127,133,136, 66,126,136,116,136,137, 63,108,114, & 130,132,136,106,115,130,132,136, 74,114,120,121, & 123,132,136, 73,100,101,114,118,127,130,131,132 /) INTEGER, PARAMETER, DIMENSION(360) :: IROW_STOICM_1 = (/ & 136, 70, 91,105,114,118,119,132,136, 91, 96, 97, & 114,118,128,132,136, 57, 98,136, 71,117,125,132, & 136, 67,110,136,106,131,132,112,114,115,133,136, & 99,106,129,131, 85, 88, 94, 98,112,114,115,126, & 127,131,133,136, 88,114,115,116,126,127,128,133, & 136,108,114,118,119,130,133,136, 88,106,114,115, & 130,131,133,136, 88,114,120,121,123,124,133,136, & 88,114,120,121,123,124,136, 88, 91,100,111,114, & 118,119,131,133,136, 88,100,101,114,118,130,131, & 133,136, 88, 91,105,114,118,119,133,136, 91, 96, & 97,114,118,128,133,136,137, 85, 88, 98,114,133, & 136, 88, 99,114,117,120,121,127,128,133,136, 88, & 110,114,115,127,128,133,136,115,132,136, 88,127, & 132,136,112,115,136,112,115, 61,112,136, 56, 94, & 136, 69,109,136, 54, 89,130,132,133,135,136, 58, & 112,113,132,135,136, 42, 64,108,114,119,122,132, & 133,135,136, 86,109,132, 86, 91,114,115,132,133, & 135,136, 79, 95,132, 42, 50, 65, 95,114,119,128, & 129,134,136, 42, 95,119,128,136, 79,114,128,135, & 136, 65, 75,128, 65, 75,128, 42, 75,128,130,134, & 42, 75,130,132,136, 42, 90, 91,100,114,132,136, & 42, 72, 91,100,132,136, 86,110,129, 90, 91,132, & 136, 91,118,130,132, 90, 91, 99,129,136, 91, 99, & 118,129,130, 14, 81,123,132,101,120,132,105,121, & 122,132,118,119,132,136, 42, 72, 89, 91,119,133, & 94,127,130,133,136,109,114,115,130,133,136, 89, & 94,127,130, 89,109,119,127,130, 42, 81, 86, 91, & 114,120,121,132,135,136, 91,114,115,118,120,132, & 135,136, 42, 91,114,118,121,132,135,136, 81,117, & 129, 97,121,129, 99,121,122,129,112,113,114,133 /) INTEGER, PARAMETER, DIMENSION(360) :: IROW_STOICM_2 = (/ & 136,114,122,130,133,136,113,114,133,114,122,133, & 71,128,131, 71,117,132,137, 67, 96,110,132, 61, & 112,115,132, 56, 94,127,132, 57, 85, 98,132, 66, & 126,127,132, 58,113,132, 69,109,119,132, 74,123, & 124,132, 47, 74, 47, 84,132, 42, 72, 84, 90, 91, & 100,114,118,119,132,136, 42, 72, 84, 90, 91,100, & 114,118,119,128,132,134,136, 76, 77,111,132, 73, & 101,131,132, 70,105,132, 42, 70,119, 64,122,132, & 64,102,132, 42,102,119,130,132,133, 89,102,130, & 131, 42, 88,102,114,119,132,133,136, 42,102,119, & 132, 42,102,119,132,136, 42,102,119,128,132,134, & 99,102,134, 54,130,132, 52, 99,112,129, 72, 91, & 100,114,118,119,124,132,135,130,133,114,130,136, & 89,114,130,133, 85, 94, 98,112,115,126,127,130, & 131,133,136,108,114,133,106,115,133, 42, 83,114, & 120,121,123,124,130,133,136, 77, 78, 90, 91,100, & 111,118,119,127,130,133,136,100,101,114,118,130, & 133,136, 91,105,114,118,119,130,133,136, 85, 98, & 130,133,136,114,115,116,126,127,128,130,133, 42, & 91, 96, 97,114,118,128,130,133,136, 99,114,117, & 120,121,128,130,133,136,110,114,115,128,130,133, & 89,126,130,131, 89,108,118,130, 89,106,130,131, & 89,123,130,131, 89,111,130,131, 89,101,130,131, & 89,105,130,131, 89,116,127,128,130, 89, 97,127, & 128,130, 89,117,127,128,130, 89,110,127,128,130, & 85, 89, 98,130,112,115,130,133,136, 89,112,115, & 130,112,113,130,133,114,122,133,128,129,134,128, & 91, 51, 81, 85, 86, 44, 52,131,115,114,135, 99, & 18, 22, 26, 30, 5, 43,128, 9, 43, 50, 6, 43, & 135, 7, 43, 49, 41, 43, 65, 3, 43, 99, 1, 43 /) INTEGER, PARAMETER, DIMENSION(360) :: IROW_STOICM_3 = (/ & 114, 4, 43, 55, 2, 43, 46, 8, 43, 79, 10, 43, & 45, 11, 43,137, 35, 43, 47, 34, 43, 74, 32, 43, & 125, 33, 43, 82, 82, 93,104,114,132,135, 68, 96, & 100,119,125,132,135, 38, 43, 96, 36, 43,104, 37, & 43, 93, 39, 43, 72, 40, 43, 89, 53, 99,128, 99, & 129, 55, 99, 48, 60,114,132,133, 13, 48, 60,132, & 133, 48, 60, 99,114,129,133, 12, 48,132,136,111, & 124,132, 77, 78, 90, 91,100,111,118,119,127,128, & 134,136, 99,111,134, 76,100,111,118,132,136,103, & 125,132, 82,107,132, 68, 72, 96,100,103,104,114, & 119,128,129,131,134,136, 99,103,134, 68, 72, 80, & 96,100,103,104,114,119,128,129,131,132,136, 93, & 100,104,107,114,119,128,134,136, 99,107,134, 80, & 93,100,104,107,114,119,128,132,136, 80,103,107, & 116,132, 31, 72, 93,114,118,129,132, 92,104,132, & 42, 72, 89, 92,114,118,119,128,129,134, 42, 72, & 89, 92,114,118,119,128,129,132,136, 50, 92,128, & 50, 92,128, 77, 91,119,131,132,136, 83, 90,100, & 118,119,128,134,136, 83, 99,134, 83, 90,100,118, & 119,132, 78, 87,132, 42, 87,127,128,134,136, 87, & 99,134, 42, 66, 87,127,132,136, 42, 78,131,132, & 135,136, 42, 68,114,128,132, 96,118,128,132, 17, & 18, 21, 22, 25, 26, 15, 17, 16, 17, 19, 21, 20, & 21, 23, 25, 24, 25, 29, 30, 27, 29, 28, 29,132, & 135,128,134,135, 46,132, 62,114,132,136, 91,114, & 136, 91,114, 99,128,132, 53,132,134, 59,129,132, & 128,129,135,129,134, 55,128,129, 55,129,134,135, & 59,128,136, 91,115,133,136, 91,115, 49,128,129, & 130,133, 91,112,127,136, 85,130,133, 85, 91,133, & 112,113,130,131,133, 91,100,114,136, 90, 91,114 /) INTEGER, PARAMETER, DIMENSION(155) :: IROW_STOICM_4 = (/ & 90, 91,136, 91,118,130,136, 91,115,118, 86, 91, & 120, 91,114,120,130,136,120,122,133,121,122,136, & 91,114,121,130,136,114,119,130,136, 71,127,128, & 132,136, 67,127,128,132,136, 61,115,132,136, 56, & 127,132,136, 57, 85,132,136, 66,127,132,136, 69, & 114,115,132,136, 58,115,132,136, 74,114,120,121, & 124,132,136, 76, 91,100,118,119,132,136, 80,127, & 128,132,136, 73,100,114,118,130,132,136, 70, 91, & 114,119,132,136, 64,114,130,132, 85, 94, 98,112, & 115,127,128,131,133,136,137, 54,132,133, 91,104, & 114,118,119,128,136, 93,100,128,130, 82,124,128, & 136,124,125,128,136, 91, 96,114,128,133, 63,114, & 130,132, 42, 65,128,129,130, 42, 65,129,130 /) INTEGER, PARAMETER, DIMENSION(1595) :: IROW_STOICM = (/& IROW_STOICM_0, IROW_STOICM_1, IROW_STOICM_2, IROW_STOICM_3, IROW_STOICM_4 /) INTEGER, PARAMETER, DIMENSION(360) :: ICOL_STOICM_0 = (/ & 1, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, & 5, 5, 5, 5, 6, 6, 7, 7, 8, 8, 9, 9, & 9, 10, 10, 10, 10, 11, 11, 12, 12, 13, 13, 13, & 13, 14, 14, 15, 15, 15, 15, 15, 16, 16, 16, 17, & 17, 18, 18, 18, 19, 19, 19, 20, 20, 21, 21, 21, & 22, 22, 23, 23, 23, 23, 24, 24, 24, 25, 25, 25, & 26, 26, 26, 27, 27, 27, 28, 28, 28, 29, 29, 29, & 30, 30, 30, 31, 31, 31, 31, 32, 32, 32, 33, 33, & 33, 33, 34, 34, 34, 35, 35, 35, 36, 36, 36, 36, & 37, 37, 37, 38, 38, 39, 39, 39, 39, 39, 40, 40, & 40, 40, 40, 40, 41, 41, 41, 41, 42, 42, 42, 43, & 43, 43, 44, 44, 44, 44, 44, 45, 45, 45, 46, 46, & 46, 46, 46, 47, 47, 47, 48, 48, 48, 49, 49, 49, & 49, 49, 50, 50, 50, 50, 50, 50, 51, 51, 51, 52, & 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 52, 53, & 53, 53, 54, 54, 54, 54, 54, 54, 54, 55, 55, 55, & 55, 55, 55, 56, 56, 56, 56, 56, 56, 57, 57, 57, & 57, 57, 57, 57, 57, 57, 57, 57, 58, 58, 58, 58, & 58, 58, 58, 58, 58, 59, 59, 59, 59, 59, 59, 59, & 59, 59, 60, 60, 60, 60, 60, 60, 60, 60, 61, 61, & 61, 61, 61, 62, 62, 62, 62, 62, 62, 63, 63, 63, & 63, 63, 63, 64, 64, 64, 64, 64, 65, 65, 65, 65, & 65, 66, 66, 66, 66, 67, 67, 67, 68, 68, 68, 68, & 69, 69, 69, 70, 70, 70, 71, 71, 71, 72, 72, 72, & 73, 73, 73, 74, 74, 74, 74, 75, 75, 75, 75, 75, & 75, 75, 75, 76, 76, 76, 76, 77, 77, 77, 78, 78, & 78, 79, 79, 79, 79, 79, 79, 80, 80, 80, 80, 80, & 80, 80, 80, 81, 81, 81, 82, 82, 82, 83, 83, 83, & 83, 83, 83, 84, 84, 84, 84, 84, 85, 85, 85, 85, & 85, 85, 85, 86, 86, 86, 86, 86, 86, 86, 86, 86 /) INTEGER, PARAMETER, DIMENSION(360) :: ICOL_STOICM_1 = (/ & 86, 87, 87, 87, 87, 87, 87, 87, 87, 88, 88, 88, & 88, 88, 88, 88, 88, 89, 89, 89, 90, 90, 90, 90, & 90, 91, 91, 91, 92, 92, 92, 93, 93, 93, 93, 93, & 94, 94, 94, 94, 95, 95, 95, 95, 95, 95, 95, 95, & 95, 95, 95, 95, 96, 96, 96, 96, 96, 96, 96, 96, & 96, 97, 97, 97, 97, 97, 97, 97, 98, 98, 98, 98, & 98, 98, 98, 98, 99, 99, 99, 99, 99, 99, 99, 99, & 100,100,100,100,100,100,100,101,101,101,101,101, & 101,101,101,101,101,102,102,102,102,102,102,102, & 102,102,103,103,103,103,103,103,103,103,104,104, & 104,104,104,104,104,104,104,105,105,105,105,105, & 105,106,106,106,106,106,106,106,106,106,106,107, & 107,107,107,107,107,107,107,108,108,108,109,109, & 109,109,110,110,110,111,111,112,112,112,113,113, & 113,114,114,114,115,115,115,115,115,115,115,116, & 116,116,116,116,116,117,117,117,117,117,117,117, & 117,117,117,118,118,118,119,119,119,119,119,119, & 119,119,120,120,120,121,121,121,121,121,121,121, & 121,121,121,122,122,122,122,122,123,123,123,123, & 123,124,124,124,125,125,125,126,126,126,126,126, & 127,127,127,127,127,128,128,128,128,128,128,128, & 129,129,129,129,129,129,130,130,130,131,131,131, & 131,132,132,132,132,133,133,133,133,133,134,134, & 134,134,134,135,135,135,135,136,136,136,137,137, & 137,137,138,138,138,138,139,139,139,139,139,139, & 140,140,140,140,140,141,141,141,141,141,141,142, & 142,142,142,143,143,143,143,143,144,144,144,144, & 144,144,144,144,144,144,145,145,145,145,145,145, & 145,145,146,146,146,146,146,146,146,146,147,147, & 147,148,148,148,149,149,149,149,150,150,150,150 /) INTEGER, PARAMETER, DIMENSION(360) :: ICOL_STOICM_2 = (/ & 150,151,151,151,151,151,152,152,152,153,153,153, & 154,154,154,155,155,155,155,156,156,156,156,157, & 157,157,157,158,158,158,158,159,159,159,159,160, & 160,160,160,161,161,161,162,162,162,162,163,163, & 163,163,164,164,165,165,165,166,166,166,166,166, & 166,166,166,166,166,166,167,167,167,167,167,167, & 167,167,167,167,167,167,167,168,168,168,168,169, & 169,169,169,170,170,170,171,171,171,172,172,172, & 173,173,173,174,174,174,174,174,174,175,175,175, & 175,176,176,176,176,176,176,176,176,177,177,177, & 177,178,178,178,178,178,179,179,179,179,179,179, & 180,180,180,181,181,181,182,182,182,182,183,183, & 183,183,183,183,183,183,183,184,184,185,185,185, & 186,186,186,186,187,187,187,187,187,187,187,187, & 187,187,187,188,188,188,189,189,189,190,190,190, & 190,190,190,190,190,190,190,191,191,191,191,191, & 191,191,191,191,191,191,191,192,192,192,192,192, & 192,192,193,193,193,193,193,193,193,193,194,194, & 194,194,194,195,195,195,195,195,195,195,195,196, & 196,196,196,196,196,196,196,196,196,197,197,197, & 197,197,197,197,197,197,198,198,198,198,198,198, & 199,199,199,199,200,200,200,200,201,201,201,201, & 202,202,202,202,203,203,203,203,204,204,204,204, & 205,205,205,205,206,206,206,206,206,207,207,207, & 207,207,208,208,208,208,208,209,209,209,209,209, & 210,210,210,210,211,211,211,211,211,212,212,212, & 212,213,213,213,213,214,214,214,215,215,216,217, & 218,219,220,221,222,223,224,225,226,227,228,229, & 230,231,232,233,234,234,234,235,235,235,236,236, & 236,237,237,237,238,238,238,239,239,239,240,240 /) INTEGER, PARAMETER, DIMENSION(360) :: ICOL_STOICM_3 = (/ & 240,241,241,241,242,242,242,243,243,243,244,244, & 244,245,245,245,246,246,246,247,247,247,248,248, & 248,249,249,249,250,250,250,250,250,250,251,251, & 251,251,251,251,251,252,252,252,253,253,253,254, & 254,254,255,255,255,256,256,256,257,257,257,258, & 258,259,259,260,260,260,260,260,261,261,261,261, & 261,262,262,262,262,262,262,263,263,263,263,264, & 264,264,265,265,265,265,265,265,265,265,265,265, & 265,265,266,266,266,267,267,267,267,267,267,268, & 268,268,269,269,269,270,270,270,270,270,270,270, & 270,270,270,270,270,270,271,271,271,272,272,272, & 272,272,272,272,272,272,272,272,272,272,272,273, & 273,273,273,273,273,273,273,273,274,274,274,275, & 275,275,275,275,275,275,275,275,275,276,276,276, & 276,276,277,277,277,277,277,277,277,278,278,278, & 279,279,279,279,279,279,279,279,279,279,280,280, & 280,280,280,280,280,280,280,280,280,281,281,281, & 282,282,282,283,283,283,283,283,283,284,284,284, & 284,284,284,284,284,285,285,285,286,286,286,286, & 286,286,287,287,287,288,288,288,288,288,288,289, & 289,289,290,290,290,290,290,290,291,291,291,291, & 291,291,292,292,292,292,292,293,293,293,293,294, & 294,295,295,296,296,297,297,298,298,299,299,300, & 300,301,301,302,302,303,303,304,304,305,305,306, & 306,307,307,307,308,308,309,309,309,309,310,310, & 310,311,311,312,312,312,313,313,313,314,314,314, & 315,315,315,316,316,317,317,317,318,318,318,318, & 319,319,319,320,320,320,320,321,321,322,322,322, & 322,322,323,323,323,323,324,324,324,325,325,325, & 326,326,326,326,326,327,327,327,327,328,328,328 /) INTEGER, PARAMETER, DIMENSION(155) :: ICOL_STOICM_4 = (/ & 329,329,329,330,330,330,330,331,331,331,332,332, & 332,333,333,333,333,333,334,334,334,335,335,335, & 336,336,336,336,336,337,337,337,337,338,338,338, & 338,338,339,339,339,339,339,340,340,340,340,341, & 341,341,341,342,342,342,342,343,343,343,343,344, & 344,344,344,344,345,345,345,345,346,346,346,346, & 346,346,346,347,347,347,347,347,347,347,348,348, & 348,348,348,349,349,349,349,349,349,349,350,350, & 350,350,350,350,351,351,351,351,352,352,352,352, & 352,352,352,352,352,352,352,353,353,353,354,354, & 354,354,354,354,354,355,355,355,355,356,356,356, & 356,357,357,357,357,358,358,358,358,358,359,359, & 359,359,360,360,360,360,360,361,361,361,361 /) INTEGER, PARAMETER, DIMENSION(1595) :: ICOL_STOICM = (/& ICOL_STOICM_0, ICOL_STOICM_1, ICOL_STOICM_2, ICOL_STOICM_3, ICOL_STOICM_4 /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_0 = (/ & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -2.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -2.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -2.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -2.000000e+00_dp, 2.000000e+00_dp, & -2.000000e+00_dp, 2.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 2.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 5.000000e-02_dp, 5.000000e-02_dp, & -1.000000e+00_dp, 9.500000e-01_dp, -1.000000e+00_dp, 5.000000e-02_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_1 = (/ & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 3.200000e-01_dp, 5.000000e-02_dp, 1.800000e-01_dp, & 3.200000e-01_dp, 3.200000e-01_dp, -1.000000e+00_dp, 1.300000e-01_dp, 1.000000e+00_dp, & 1.900000e-01_dp, 1.800000e-01_dp, -1.000000e+00_dp, 2.700000e-01_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 3.900000e-01_dp, 7.500000e-01_dp, -1.000000e+00_dp, & 3.000000e-01_dp, 5.700000e-01_dp, 2.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 9.600000e-01_dp, 9.600000e-01_dp, 9.600000e-01_dp, -1.000000e+00_dp, 4.000000e-02_dp, & -1.000000e+00_dp, 9.300000e-01_dp, 9.300000e-01_dp, 9.300000e-01_dp, -1.000000e+00_dp, & 7.000000e-02_dp, 4.700000e-02_dp, 1.000000e-01_dp, 6.600000e-01_dp, 4.000000e-01_dp, & 2.600000e-01_dp, -1.000000e+00_dp, 1.230000e-01_dp, 7.000000e-02_dp, 8.830000e-01_dp, & -1.000000e+00_dp, 7.830000e-01_dp, 1.200000e-01_dp, 5.300000e-01_dp, -1.000000e+00_dp, & 3.500000e-01_dp, 3.500000e-01_dp, 8.800000e-01_dp, 5.300000e-01_dp, -1.000000e+00_dp, & 3.500000e-01_dp, 8.570000e-01_dp, 1.500000e-01_dp, -1.000000e+00_dp, 1.430000e-01_dp, & 1.430000e-01_dp, 8.570000e-01_dp, 8.500000e-01_dp, -1.000000e+00_dp, 8.500000e-01_dp, & 5.000000e-01_dp, 5.000000e-01_dp, -1.000000e+00_dp, 5.000000e-01_dp, 5.000000e-01_dp, & 1.500000e+00_dp, 5.000000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 2.000000e-01_dp, & 8.000000e-01_dp, 1.200000e+00_dp, -1.000000e+00_dp, 8.000000e-01_dp, -1.000000e+00_dp, & 2.000000e-01_dp, 8.000000e-01_dp, 1.200000e+00_dp, -1.000000e+00_dp, 8.000000e-01_dp, & -2.000000e+00_dp, 2.800000e-01_dp, 1.720000e+00_dp, 2.800000e-01_dp, 1.120000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 2.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 5.000000e-01_dp, 5.000000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_2 = (/ & 5.000000e-01_dp, 5.000000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 2.500000e-01_dp, -1.000000e+00_dp, & 7.500000e-01_dp, 7.500000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, 2.500000e-01_dp, & -1.000000e+00_dp, 1.250000e+00_dp, 5.000000e-01_dp, 1.600000e-01_dp, 9.000000e-02_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 8.500000e-01_dp, -1.000000e+00_dp, & 1.500000e-01_dp, 1.500000e-01_dp, 1.500000e-01_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 8.800000e-01_dp, & 1.200000e-01_dp, 7.300000e-02_dp, 4.700000e-02_dp, -1.000000e+00_dp, 1.200000e-01_dp, & -8.800000e-01_dp, 1.000000e-01_dp, 5.780000e-01_dp, -1.000000e+00_dp, 1.020000e-01_dp, & 1.020000e-01_dp, 3.300000e-02_dp, 5.780000e-01_dp, 1.870000e-01_dp, 6.800000e-01_dp, & -8.980000e-01_dp, 2.000000e-02_dp, 6.860000e-01_dp, -1.000000e+00_dp, 2.940000e-01_dp, & 2.940000e-01_dp, 6.860000e-01_dp, 9.800000e-01_dp, -2.000000e-02_dp, 5.000000e-01_dp, & 5.000000e-01_dp, -1.000000e+00_dp, 5.000000e-01_dp, 5.000000e-01_dp, 5.000000e-01_dp, & 1.000000e+00_dp, -5.000000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 5.000000e-01_dp, -1.000000e+00_dp, 5.000000e-01_dp, 5.000000e-01_dp, -5.000000e-01_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 7.500000e-01_dp, 7.500000e-01_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.600000e-01_dp, 2.500000e-01_dp, 3.000000e-02_dp, 9.000000e-02_dp, 1.600000e-01_dp, & 7.500000e-01_dp, 1.600000e-01_dp, -1.000000e+00_dp, 7.000000e-02_dp, 3.500000e-01_dp, & -9.100000e-01_dp, 6.400000e-01_dp, 2.500000e-01_dp, 9.500000e-01_dp, 3.800000e-01_dp, & -1.000000e+00_dp, 1.500000e-01_dp, 5.400000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 5.000000e-01_dp, -1.000000e+00_dp, 5.000000e-01_dp, 5.000000e-01_dp, 2.000000e-01_dp, & 3.000000e-01_dp, -1.000000e+00_dp, 3.000000e-01_dp, 2.500000e-01_dp, -1.000000e+00_dp, & 7.500000e-01_dp, 5.000000e-01_dp, 5.000000e-01_dp, 2.500000e-01_dp, -1.000000e+00_dp, & 5.000000e-01_dp, 2.400000e-01_dp, 1.220000e+00_dp, 2.800000e-01_dp, 1.800000e-01_dp, & -1.000000e+00_dp, 3.000000e-01_dp, -1.000000e+00_dp, 1.100000e+00_dp, 4.800000e-01_dp, & 9.200000e-01_dp, 5.600000e-01_dp, 3.600000e-01_dp, -2.000000e+00_dp, 5.000000e-01_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_3 = (/ & 1.280000e+00_dp, 2.500000e-01_dp, 3.300000e-01_dp, 1.300000e-01_dp, -1.000000e+00_dp, & 9.500000e-01_dp, 2.900000e-01_dp, 1.800000e-01_dp, 2.500000e-01_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 2.500000e-01_dp, 3.600000e-01_dp, -1.000000e+00_dp, 8.900000e-01_dp, & 1.400000e-01_dp, 3.600000e-01_dp, 2.500000e-01_dp, -1.000000e+00_dp, 6.400000e-01_dp, & 1.500000e-01_dp, 5.950000e-01_dp, -1.000000e+00_dp, 1.255000e+00_dp, 2.550000e-01_dp, & 5.950000e-01_dp, -1.000000e+00_dp, 1.700000e+00_dp, 3.750000e-01_dp, 3.750000e-01_dp, & -1.000000e+00_dp, 6.250000e-01_dp, 3.750000e-01_dp, 3.750000e-01_dp, -1.000000e+00_dp, & 3.750000e-01_dp, 2.500000e-01_dp, 7.500000e-01_dp, 2.500000e-01_dp, -1.000000e+00_dp, & 7.500000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, 2.500000e-01_dp, 4.250000e-01_dp, & 8.300000e-01_dp, -1.000000e+00_dp, 3.000000e-02_dp, 5.000000e-02_dp, 2.500000e-01_dp, & 5.750000e-01_dp, -1.000000e+00_dp, 9.000000e-01_dp, 2.500000e-01_dp, -1.000000e+00_dp, & 1.250000e+00_dp, 5.000000e-01_dp, 2.500000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 5.000000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -2.000000e+00_dp, 2.000000e+00_dp, & 2.000000e+00_dp, -2.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 4.100000e-01_dp, 1.500000e-01_dp, -1.000000e+00_dp, & 4.400000e-01_dp, 4.400000e-01_dp, 1.500000e-01_dp, -1.000000e+00_dp, 4.100000e-01_dp, & 4.400000e-01_dp, -1.000000e+00_dp, 4.400000e-01_dp, 1.500000e-01_dp, -1.000000e+00_dp, & 8.700000e-01_dp, 1.000000e-01_dp, 1.520000e-01_dp, 4.350000e-01_dp, 2.830000e-01_dp, & -1.000000e+00_dp, 8.700000e-01_dp, 4.350000e-01_dp, 3.000000e-02_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 4.200000e-01_dp, & 5.350000e-01_dp, 5.000000e-01_dp, 1.350000e-01_dp, 3.050000e-01_dp, -1.000000e+00_dp, & 3.000000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 5.400000e-01_dp, & 1.000000e-01_dp, 3.600000e-01_dp, -1.000000e+00_dp, 5.400000e-01_dp, 3.600000e-01_dp, & 9.000000e-01_dp, 5.400000e-01_dp, -1.000000e+00_dp, 5.400000e-01_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 6.000000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_4 = (/ & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 2.600000e-01_dp, 1.550000e-01_dp, 5.800000e-01_dp, -1.000000e+00_dp, 8.500000e-01_dp, & -7.400000e-01_dp, 7.400000e-01_dp, 8.100000e-01_dp, 1.000000e+00_dp, 1.900000e-01_dp, & -1.000000e+00_dp, -8.100000e-01_dp, 8.100000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 2.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 2.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 4.700000e-01_dp, -1.000000e+00_dp, 5.300000e-01_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 5.000000e-01_dp, & 5.000000e-01_dp, 5.000000e-01_dp, 5.000000e-01_dp, -1.000000e+00_dp, 5.000000e-01_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 6.500000e-01_dp, 3.500000e-01_dp, -1.000000e+00_dp, & 1.500000e-01_dp, -1.000000e+00_dp, 7.000000e-02_dp, 5.000000e-02_dp, 9.000000e-01_dp, & 1.590000e-01_dp, 3.870000e-01_dp, 2.700000e-01_dp, -9.000000e-01_dp, 6.000000e-02_dp, & 5.000000e-02_dp, 8.000000e-01_dp, 4.000000e-02_dp, 8.200000e-01_dp, -1.000000e+00_dp, & 8.000000e-02_dp, -8.000000e-01_dp, 6.000000e-02_dp, 1.600000e-01_dp, 2.000000e-01_dp, & 7.000000e-01_dp, 8.000000e-01_dp, -1.000000e+00_dp, 2.150000e-01_dp, -8.000000e-01_dp, & 2.750000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 2.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 3.600000e-01_dp, -3.600000e-01_dp, 6.400000e-01_dp, -1.000000e+00_dp, & 7.910000e-01_dp, 2.090000e-01_dp, -2.090000e-01_dp, -1.000000e+00_dp, 3.600000e-01_dp, & 6.400000e-01_dp, -3.600000e-01_dp, -1.000000e+00_dp, 3.600000e-01_dp, 6.400000e-01_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_5 = (/ & -3.600000e-01_dp, -1.000000e+00_dp, 7.910000e-01_dp, 2.090000e-01_dp, -2.090000e-01_dp, & -1.000000e+00_dp, 2.090000e-01_dp, 7.910000e-01_dp, -2.090000e-01_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 2.090000e-01_dp, 7.910000e-01_dp, & -2.090000e-01_dp, -1.000000e+00_dp, 3.870000e-01_dp, 6.130000e-01_dp, -3.870000e-01_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 2.000000e-01_dp, 7.400000e-02_dp, -1.000000e+00_dp, 2.750000e-01_dp, 2.510000e-01_dp, & 2.750000e-01_dp, 3.750000e-01_dp, 2.750000e-01_dp, 7.250000e-01_dp, 1.125000e+00_dp, & -1.750000e-01_dp, 2.000000e-01_dp, 7.400000e-02_dp, -1.000000e+00_dp, 2.750000e-01_dp, & 2.510000e-01_dp, 2.750000e-01_dp, 3.750000e-01_dp, 2.750000e-01_dp, 7.250000e-01_dp, & 1.000000e+00_dp, 1.250000e-01_dp, -1.000000e+00_dp, 8.250000e-01_dp, -1.000000e+00_dp, & 6.540000e-01_dp, 3.460000e-01_dp, -3.460000e-01_dp, -1.000000e+00_dp, 2.090000e-01_dp, & 7.910000e-01_dp, -2.090000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 2.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 7.000000e-01_dp, & 3.000000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, 7.000000e-01_dp, 7.000000e-01_dp, & -1.000000e+00_dp, 7.000000e-01_dp, 2.000000e+00_dp, -2.000000e+00_dp, 2.000000e+00_dp, & 2.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 2.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 2.000000e-01_dp, 4.000000e-01_dp, 2.800000e-01_dp, & 1.200000e-01_dp, 6.000000e-01_dp, 2.000000e-01_dp, -1.000000e+00_dp, 1.000000e-01_dp, & -7.000000e-01_dp, -2.000000e+00_dp, 2.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 3.200000e-01_dp, 5.000000e-02_dp, 1.800000e-01_dp, 3.200000e-01_dp, 3.200000e-01_dp, & -1.000000e+00_dp, 1.300000e-01_dp, -1.000000e+00_dp, 1.900000e-01_dp, 1.000000e+00_dp, & 2.700000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.130000e-01_dp, 7.470000e-01_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_6 = (/ & 4.530000e-01_dp, 2.940000e-01_dp, -1.000000e+00_dp, 1.400000e-01_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 8.870000e-01_dp, 2.900000e-01_dp, 1.700000e-01_dp, 2.160000e-01_dp, & 9.000000e-02_dp, 2.340000e-01_dp, -1.000000e+00_dp, 2.340000e-01_dp, 2.160000e-01_dp, & 9.000000e-02_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 6.000000e-01_dp, & -1.000000e+00_dp, 4.000000e-01_dp, 4.000000e-01_dp, -4.000000e-01_dp, 1.000000e+00_dp, & 4.000000e-01_dp, 8.570000e-01_dp, -1.000000e+00_dp, 1.430000e-01_dp, 1.430000e-01_dp, & 8.570000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, 8.500000e-01_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 3.900000e-01_dp, & 7.500000e-01_dp, -1.000000e+00_dp, 3.000000e-01_dp, 5.700000e-01_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 5.000000e-01_dp, 5.000000e-01_dp, & -1.000000e+00_dp, 5.000000e-01_dp, 5.000000e-01_dp, 5.000000e-01_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 5.000000e-01_dp, 8.500000e-01_dp, 1.500000e-01_dp, -1.000000e+00_dp, & 5.000000e-02_dp, 1.000000e-01_dp, 1.500000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 8.000000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 2.000000e+00_dp, -2.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_7 = (/ & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 6.100000e-01_dp, & 3.900000e-01_dp, 1.000000e+00_dp, 2.700000e-01_dp, -1.000000e+00_dp, 5.000000e-01_dp, & 5.000000e-01_dp, 5.000000e-01_dp, 5.000000e-01_dp, -1.000000e+00_dp, 2.700000e-01_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 5.000000e-01_dp, 5.000000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 2.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 2.500000e-01_dp, 7.500000e-01_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 2.900000e-01_dp, 1.700000e-01_dp, 2.160000e-01_dp, 9.000000e-02_dp, & 2.340000e-01_dp, -1.000000e+00_dp, 2.340000e-01_dp, 2.160000e-01_dp, 9.000000e-02_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e-01_dp, 9.000000e-01_dp, -1.000000e+00_dp, 9.000000e-01_dp, & 9.000000e-01_dp, -1.000000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.300000e-01_dp, 3.100000e-01_dp, & 1.500000e-01_dp, 1.500000e-01_dp, -1.000000e+00_dp, 7.000000e-02_dp, 6.500000e-01_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_8 = (/ & 4.400000e-01_dp, 1.340000e+00_dp, 3.100000e-01_dp, 3.400000e-01_dp, -1.000000e+00_dp, & 3.500000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 6.500000e-02_dp, & 1.550000e-01_dp, 5.000000e-01_dp, 7.500000e-02_dp, 7.500000e-02_dp, -1.000000e+00_dp, & 3.500000e-02_dp, 3.250000e-01_dp, 2.200000e-01_dp, 1.700000e-01_dp, 1.550000e-01_dp, & 1.700000e-01_dp, 5.000000e-01_dp, -8.250000e-01_dp, 1.400000e-01_dp, 6.000000e-01_dp, & 2.600000e-01_dp, -1.000000e+00_dp, 4.000000e-01_dp, 6.000000e-01_dp, 1.600000e+00_dp, & -1.000000e+00_dp, 4.000000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 5.000000e-01_dp, 7.000000e-02_dp, 3.000000e-01_dp, 1.300000e-01_dp, -1.000000e+00_dp, & 2.000000e-01_dp, 3.000000e-01_dp, 3.000000e-01_dp, 5.000000e-01_dp, -8.000000e-01_dp, & -1.000000e+00_dp, 1.930000e-01_dp, 1.930000e-01_dp, 6.120000e-01_dp, -3.880000e-01_dp, & 3.500000e-01_dp, 6.500000e-01_dp, -1.000000e+00_dp, 3.500000e-01_dp, 6.500000e-01_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 9.300000e-01_dp, 7.000000e-02_dp, 8.000000e-02_dp, -1.000000e+00_dp, 8.000000e-02_dp, & 7.000000e-02_dp, 8.500000e-01_dp, 1.850000e+00_dp, 1.500000e-01_dp, -1.000000e+00_dp, & 9.300000e-01_dp, 7.000000e-02_dp, 8.000000e-02_dp, -1.000000e+00_dp, 8.000000e-02_dp, & 7.000000e-02_dp, 8.500000e-01_dp, 8.500000e-01_dp, 1.500000e-01_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.500000e+00_dp, 5.000000e-01_dp, & 5.000000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 4.800000e-01_dp, & 5.200000e-01_dp, 5.200000e-01_dp, 4.800000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 4.800000e-01_dp, 5.200000e-01_dp, 5.200000e-01_dp, 4.800000e-01_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 5.000000e-01_dp, 5.000000e-01_dp, -1.000000e+00_dp, & 5.000000e-01_dp, 5.000000e-01_dp, -5.000000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_9 = (/ & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 2.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 2.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 2.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 6.000000e-01_dp, & 4.000000e-01_dp, 6.000000e-01_dp, 4.000000e-01_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 2.000000e+00_dp, 8.500000e-01_dp, 1.500000e-01_dp, & 8.500000e-01_dp, -1.000000e+00_dp, 1.500000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 2.000000e+00_dp, -1.000000e+00_dp, 1.500000e+00_dp, 5.000000e-01_dp, & -1.000000e+00_dp, 2.000000e+00_dp, 2.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp /) REAL(kind=dp), PARAMETER, DIMENSION(95) :: STOICM_10 = (/ & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 7.100000e-01_dp, & 4.250000e-01_dp, 2.850000e-01_dp, 2.900000e-01_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 6.700000e-01_dp, 2.600000e-01_dp, 5.800000e-01_dp, 3.600000e-01_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 7.000000e-01_dp, 3.000000e-01_dp, & 3.000000e-01_dp, 7.000000e-01_dp, 1.000000e+00_dp, 3.000000e-01_dp, -1.000000e+00_dp, & 5.000000e-01_dp, 5.000000e-01_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 3.200000e-01_dp, & 5.000000e-02_dp, 1.800000e-01_dp, 3.200000e-01_dp, 3.200000e-01_dp, 1.300000e-01_dp, & 1.000000e+00_dp, 1.900000e-01_dp, 1.800000e-01_dp, 2.700000e-01_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 5.000000e-01_dp, -1.000000e+00_dp, & 5.000000e-01_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 7.000000e-01_dp, 3.000000e-01_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp /) REAL(kind=dp), PARAMETER, DIMENSION(1595) :: STOICM = (/& STOICM_0, STOICM_1, STOICM_2, STOICM_3, STOICM_4, & STOICM_5, STOICM_6, STOICM_7, STOICM_8, STOICM_9, & STOICM_10 /) END MODULE gckpp_StoichiomSP