! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! 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 : Wed Sep 15 16:00:44 2010 ! Working directory : /mnt/lstr04/srv/home/c/ccarouge/KPP/geoschem_kppfiles/v8-03-02/SOA ! 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(325) :: CROW_JVRP = (/ & 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, 47, 49, 51, 53, 55, 57, 58, 59, 60, & 62, 64, 66, 68, 70, 71, 73, 75, 77, 79, 81, 83, & 85, 87, 89, 91, 93, 95, 97, 99,101,103,105,107, & 109,111,113,115,117,119,121,123,125,127,129,130, & 132,134,135,137,138,140,141,142,142,144,146,148, & 149,151,153,155,157,159,161,163,165,167,169,171, & 173,175,177,179,181,183,185,187,189,191,193,195, & 197,199,201,203,205,207,209,211,213,215,217,219, & 221,223,225,226,227,228,229,231,233,235,237,239, & 241,243,244,246,248,250,252,254,256,257,259,260, & 262,264,266,268,270,272,274,276,278,280,282,284, & 286,288,290,292,294,296,297,299,301,303,304,306, & 308,310,312,314,316,318,320,322,323,325,327,329, & 331,333,335,337,339,340,342,344,345,347,349,351, & 353,355,357,359,361,363,365,367,369,371,373,375, & 377,379,381,383,385,387,389,391,393,395,397,399, & 401,403,405,407,409,411,413,415,417,419,420,421, & 421,421,421,421,421,421,421,421,421,421,421,421, & 421,421,421,421,421,421,422,423,424,425,426,427, & 428,429,430,431,432,433,434,435,437,439,441,443, & 445,447,449,450,452,454,456,458,460,462,463,464, & 465,466,467,468,469,470,471,472,473,474,475,476, & 477,478,479,480,481,482,483,483,484,484,484,485, & 486,487,488,489,490,491,492,493,494,495,496,497, & 498,499,500,501,501,502,503,504,505,506,507,508, & 509 /) INTEGER, PARAMETER, DIMENSION(360) :: ICOL_JVRP_0 = (/ & 98,103, 98,102, 96, 98, 98, 99, 95, 98,102,102, & 96,102, 29,102, 96,103, 96,102, 60,102,102, 95, & 103, 95, 96, 95, 95, 39,102, 39,102, 82,102, 99, & 102, 67,102,102,103, 37,102, 96, 99, 49, 49,102, & 96,100,100,103,100,102, 99,100, 38,102,102, 99, & 100, 82,100, 84,102, 84,100, 99,101, 34,101,103, & 36,102, 87,103, 26,102, 26,102, 65,103, 83,103, & 35,102, 85,103, 85,103, 68,103, 74,103, 73,103, & 92,103, 70,103, 70,103, 79,103, 80,103, 75,103, & 75,103, 72,103, 72,103, 64,103, 69,103, 66,103, & 78,103, 76,103, 35,100, 86,102,102, 81,102, 90, & 99, 27, 94, 99, 28, 97, 99, 58, 99, 90,103, 94, & 103, 97,103,103, 81,100, 62,102, 62,102, 65, 95, & 83, 95, 85, 96, 68, 96, 74, 96, 73, 96, 92, 96, & 70, 96, 79, 96, 80, 96, 75, 96, 72, 96, 64, 96, & 69, 96, 66, 96, 78, 96, 76, 96, 89,102, 87, 95, & 89,100, 85, 95, 68, 95, 74, 95, 73, 95, 92, 95, & 70, 95, 79, 95, 80, 95, 75, 95, 72, 95, 64, 95, & 69, 95, 66, 95, 78, 95, 76, 95,102,102, 87, 87, & 87, 96, 65, 96, 83, 96, 96,101, 90, 96, 94, 96, & 96, 97, 96, 61,102, 61, 98, 58,102, 58, 98, 63, & 102, 61,100,102, 88,102,100, 88,100, 59,102, 93, & 102, 91,102, 77,102, 65,101, 83,101, 65,101, 83, & 101, 59, 98, 93, 98, 91, 98, 59,100, 91,100, 91, & 100, 90, 95, 94, 95, 95, 97, 95, 90, 95, 94, 95, & 95, 97, 95, 52,102, 51,102, 46,102, 41,102, 42, & 102, 40,102, 43,102, 50,102, 47,102,102, 57,102, & 54,102, 56,102, 55,102, 53,102, 45,102, 48,102, & 36,100,102, 71,102, 71, 98,101, 95,101, 95,101, & 85,101, 74,101, 73,101, 92,101, 70,101, 79,101 /) INTEGER, PARAMETER, DIMENSION(148) :: ICOL_JVRP_1 = (/ & 80,101, 75,101, 72,101, 66,101, 68,101, 64,101, & 69,101, 78,101, 76,101, 85,101, 74,101, 73,101, & 92,101, 70,101, 79,101, 75,101, 72,101, 68,101, & 80,101, 64,101, 69,101, 78,101, 76,101, 66,101, & 87,101, 87,101, 90,101, 94,101, 97,101,101,100, & 99, 98, 34, 67, 82, 38, 29, 58, 27, 86, 96, 99, & 100, 38, 44,102, 44,102, 44,100, 30,102, 25,102, & 23,102, 24,102, 1, 31, 96, 31,103, 32, 96, 32, & 103, 33, 96, 33,103, 98, 99, 29, 39, 82, 82, 67, & 37, 49,100,100, 38, 38, 49, 84, 84, 34, 81, 62, & 62, 89, 63, 88, 88, 93, 93, 93, 91, 91, 77, 52, & 51, 46, 41, 42, 40, 50, 43, 47, 57, 54, 56, 55, & 53, 45, 86, 47 /) INTEGER, PARAMETER, DIMENSION(508) :: 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, 28, 28, & 29, 29, 30, 30, 31, 31, 32, 32, 33, 34, 35, 36, & 36, 37, 37, 38, 38, 39, 39, 40, 40, 41, 42, 42, & 43, 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, 64, 65, 65, 66, 66, & 67, 67, 68, 68, 69, 69, 70, 70, 71, 72, 72, 73, & 73, 74, 75, 75, 76, 77, 77, 78, 79, 81, 81, 82, & 82, 83, 83, 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,100, & 101,101,102,102,103,103,104,104,105,105,106,106, & 107,107,108,108,109,109,110,110,111,111,112,112, & 113,113,114,114,115,115,116,116,117,117,118,118, & 119,119,120,120,121,121,122,122,123,124,125,126, & 127,127,128,128,129,129,130,130,131,131,132,132, & 133,133,134,135,135,136,136,137,137,138,138,139, & 139,140,140,141,142,142,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,163,163,164,164, & 165,165,166,167,167,168,168,169,169,170,170,171, & 171,172,172,173,173,174,174,175,175,176,177,177, & 178,178,179,179,180,180,181,181,182,182,183,183, & 184,184,185,186,186,187,187,188,189,189,190,190, & 191,191,192,192,193,193,194,194,195,195,196,196 /) INTEGER, PARAMETER, DIMENSION(148) :: IROW_JVRP_1 = (/ & 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,215,216,216,217,217,218,218,219,219,220,220, & 221,221,222,222,223,223,224,224,225,225,226,227, & 246,247,248,249,250,251,252,253,254,255,256,257, & 258,259,260,260,261,261,262,262,263,263,264,264, & 265,265,266,266,267,268,268,269,269,270,270,271, & 271,272,272,273,273,274,275,276,277,278,279,280, & 281,282,283,284,285,286,287,288,289,290,291,292, & 293,294,296,299,300,301,302,303,304,305,306,307, & 308,309,310,311,312,313,314,315,317,318,319,320, & 321,322,323,323 /) INTEGER, PARAMETER, DIMENSION(508) :: IROW_JVRP = (/& IROW_JVRP_0, IROW_JVRP_1 /) ! Stoichiometric Matrix in Compressed Column Sparse Format INTEGER, PARAMETER, DIMENSION(325) :: CCOL_STOICM = (/ & 1, 4, 7, 10, 13, 17, 19, 21, 23, 26, 30, 32, & 34, 38, 40, 45, 48, 50, 53, 56, 58, 62, 65, 68, & 71, 74, 77, 80, 83, 87, 90, 94, 97,100,103,106, & 108,113,119,123,126,129,134,137,142,145,148,153, & 159,162,174,177,184,190,196,204,210,213,223,230, & 238,241,247,250,259,264,269,277,282,286,289,292, & 295,298,301,304,307,310,313,314,315,319,324,329, & 333,337,340,343,348,355,358,361,365,369,372,375, & 378,381,384,387,390,393,396,399,402,405,410,414, & 425,433,440,447,456,462,471,479,487,493,501,508, & 513,522,529,532,535,538,540,543,546,549,555,559, & 564,568,570,573,581,587,592,597,600,603,607,611, & 616,620,623,627,631,636,642,645,649,659,667,675, & 678,681,685,690,694,699,703,706,709,712,714,717, & 720,724,728,732,736,740,743,746,748,752,755,760, & 764,767,770,774,778,781,786,794,796,799,802,813, & 819,822,831,837,846,852,859,867,872,880,887,893, & 902,908,911,914,917,920,923,926,929,932,936,940, & 944,948,952,956,959,964,967,971,976,979,983,985, & 986,987,988,989,990,991,992,993,994,995,996,997, & 998,999,1000,1001,1002,1003,1006,1009,1012,1015,1018,1021, & 1024,1027,1030,1033,1035,1038,1040,1042,1047,1052,1058,1062, & 1065,1068,1071,1073,1076,1079,1082,1085,1088,1091,1093,1096, & 1098,1102,1105,1107,1110,1113,1116,1119,1121,1124,1128,1131, & 1135,1137,1142,1146,1149,1152,1157,1159,1163,1165,1167,1171, & 1174,1177,1182,1185,1188,1193,1197,1202,1207,1211,1215,1219, & 1223,1228,1232,1236,1239,1246,1253,1258,1265,1271,1275,1286, & 1289 /) INTEGER, PARAMETER, DIMENSION(360) :: IROW_STOICM_0 = (/ & 98, 99,103, 96, 98,102, 96, 98,102, 98, 99,100, & 82, 95, 96, 98, 98,102, 29,102, 96,102, 29, 96, & 102, 96, 99,102,103, 29, 96, 96,102, 21, 60, 96, & 102, 95,102, 82, 95, 96, 99,103, 39, 95, 96, 82, & 95, 82, 95, 96, 39, 95,102, 39, 82, 60, 82, 96, & 102, 67, 99,102, 67,100,102, 37,102,103, 37, 99, & 102, 49, 96, 99, 49, 96, 99, 49, 99,102, 96, 99, & 100,102, 99,100,103, 96, 99,100,102, 38, 99,100, & 38, 99,100, 21, 96,102, 82, 96,102,100,103, 60, & 67, 82, 96,100, 60, 82, 84, 96,101,102, 67, 84, & 100,101, 34, 99,101, 34, 99,101, 21, 95, 99,101, & 103, 36, 87,102, 84, 87, 96, 99,103, 26, 66,102, & 26, 65,102, 65, 81, 96, 99,103, 82, 83, 84, 96, & 99,103, 35, 85,102, 62, 65, 66, 81, 84, 85, 87, & 89, 95, 96, 99,103, 85, 86,103, 68, 81, 82, 84, & 85, 99,103, 74, 82, 86, 99,101,103, 73, 84, 86, & 99,101,103, 71, 82, 91, 92, 93, 96, 99,103, 70, & 71, 82, 96, 99,103, 67, 70,103, 60, 63, 67, 77, & 79, 82, 88, 96, 99,103, 63, 67, 77, 80, 96, 99, & 103, 63, 75, 82, 88, 96, 99,101,103, 67, 75,103, & 72, 77, 82, 96, 99,103, 67, 72,103, 63, 64, 67, & 82, 88, 96, 99,101,103, 69, 82, 88, 99,103, 62, & 66, 96, 99,103, 67, 78, 82, 91, 93, 96, 99,103, & 76, 82, 84, 99,103, 35, 67, 85,100, 68, 86,102, & 21, 95,102, 81, 90,102, 27, 90, 99, 27, 90, 99, & 28, 94, 99, 28, 94, 99, 58, 97, 99, 58, 97, 99, & 99, 99, 87, 90, 99,103, 82, 94, 96, 99,103, 82, & 97, 99,101,103, 60, 96, 99,103, 67, 81, 90,100, & 62, 74,102, 62, 74,102, 65, 81, 82, 95, 96, 77, & 81, 82, 83, 84, 95, 96, 40, 85, 96, 68, 86, 96 /) INTEGER, PARAMETER, DIMENSION(360) :: IROW_STOICM_1 = (/ & 74, 95, 96,101, 73, 88, 95, 96, 57, 92, 96, 57, & 70, 96, 54, 79, 96, 56, 80, 96, 55, 75, 96, 53, & 72, 96, 56, 64, 96, 56, 69, 96, 42, 66, 96, 52, & 78, 96, 51, 76, 96, 73, 89,102, 82, 84, 87, 95, & 96, 67, 73, 89,100, 62, 65, 66, 81, 82, 84, 85, & 87, 89, 95, 96, 68, 81, 82, 84, 85, 95, 96, 99, & 74, 77, 82, 88, 95, 96,101, 73, 82, 84, 89, 95, & 96,101, 70, 71, 82, 89, 91, 92, 93, 95, 96, 70, & 71, 82, 89, 95, 96, 60, 63, 77, 79, 82, 88, 89, & 95, 96, 63, 77, 80, 81, 82, 95, 96, 99, 63, 75, & 82, 88, 89, 95, 96,101, 60, 72, 77, 82, 95, 96, & 64, 81, 82, 88, 95, 96, 99,101, 69, 81, 82, 88, & 95, 96, 99, 62, 66, 82, 95, 96, 67, 78, 81, 82, & 91, 93, 95, 96, 99, 76, 81, 82, 84, 95, 96, 99, & 84, 96,102, 81, 96,102, 84, 87, 96, 84, 87, 46, & 87, 96, 41, 65, 96, 50, 83, 96, 48, 95, 96, 98, & 101,102, 43, 90, 96, 98, 47, 82, 94, 96, 98, 45, & 96, 97, 98, 96, 98, 61, 83,102, 60, 61, 82, 84, & 95, 96, 98,102, 58, 77, 82, 96, 99,102, 58, 82, & 96, 98, 99, 60, 63, 94, 96,102, 61, 76,100, 60, & 96,102, 60, 88,101,102, 60, 67, 96,100, 60, 67, & 88,100,101, 14, 59, 92,102, 75, 93,102, 72, 91, & 97,102, 77, 88, 96,102, 65, 81, 95, 96,101, 82, & 83, 84, 95, 96,101, 65, 81,101, 77, 81, 83,101, & 21, 59, 60, 61, 82, 91, 93, 96, 98,102, 60, 82, & 84, 88, 93, 96, 98,102, 21, 60, 82, 88, 91, 96, & 98,102, 59, 78,100, 69, 91,100, 67, 91, 97,100, & 82, 87, 90, 95, 96, 82, 94, 95, 96, 82, 95, 96, & 97,101, 60, 82, 95, 96, 82, 90, 95, 82, 94, 95, & 82, 95, 97, 82, 95, 52, 78,102, 51, 76,102, 46 /) INTEGER, PARAMETER, DIMENSION(360) :: IROW_STOICM_2 = (/ & 84, 87,102, 41, 65, 81,102, 42, 62, 66,102, 40, & 81, 85,102, 43, 84, 90,102, 50, 83,102, 47, 94, & 102, 60,102, 57, 71, 92,102, 54, 79,102, 56, 80, & 81, 99,102, 55, 75, 81,102, 53, 72,102, 45, 97, & 102, 48, 82,101,102, 36, 67, 87,100, 82, 99,102, & 71, 79, 96, 97,102, 60, 63, 71, 77, 82, 88, 98, & 102, 95,101, 82, 96,101, 82, 95,101, 62, 65, 66, & 81, 84, 85, 87, 89, 95, 96,101, 74, 82, 88, 95, & 96,101, 73, 84, 95, 70, 71, 82, 91, 92, 93, 95, & 96,101, 70, 71, 82, 95, 96,101, 60, 63, 77, 79, & 82, 88, 95, 96,101, 63, 77, 80, 95, 99,101, 63, & 75, 82, 88, 95, 96,101, 60, 72, 77, 82, 88, 95, & 96,101, 62, 66, 95, 96,101, 68, 81, 82, 84, 85, & 95, 99,101, 64, 82, 88, 95, 96, 99,101, 69, 82, & 88, 95, 99,101, 67, 78, 82, 91, 93, 95, 96, 99, & 101, 76, 82, 84, 95, 99,101, 85, 89,101, 74, 89, & 101, 73, 89,101, 89, 92,101, 70, 89,101, 79, 89, & 101, 75, 89,101, 72, 89,101, 68, 81, 99,101, 80, & 81, 99,101, 64, 81, 99,101, 69, 81, 99,101, 78, & 81, 99,101, 76, 81, 99,101, 62, 66,101, 84, 87, & 95, 96,101, 84, 87,101, 87, 90, 95,101, 82, 94, & 95, 96,101, 82, 95, 97, 60, 95, 96,101, 99,100, & 103, 99, 60, 35, 59, 62, 61, 26, 36, 89, 84, 82, & 98, 67, 25, 23, 24, 1, 6, 22, 99, 7, 22, 98, & 8, 22, 34, 4, 22, 67, 2, 22, 82, 5, 22, 38, & 3, 22, 29, 9, 22, 58, 10, 22, 27, 11, 22, 86, & 29, 96, 37, 67, 99, 67,100, 38, 67, 30, 44, 82, & 95,102, 13, 30, 44, 95,102, 30, 44, 67, 82, 95, & 100, 12, 30, 96,102, 25, 31,102, 23, 32,102, 24, & 33,102, 1, 82, 15, 31, 96, 16, 31,103, 17, 32 /) INTEGER, PARAMETER, DIMENSION(208) :: IROW_STOICM_3 = (/ & 96, 18, 32,103, 19, 33, 96, 20, 33,103, 98,102, & 98, 99,103, 29,102, 39, 82, 96,102, 60, 82, 96, & 60, 82, 67, 99,102, 37,102,103, 49,100,102, 98, & 99,100,100,103, 38, 99,100, 38, 98,100,103, 49, & 96, 99, 60, 84, 95, 96, 60, 84, 34, 95, 99,100, & 101, 60, 81, 87, 96, 62, 95,101, 60, 62, 95, 87, & 89, 90, 95,101, 82, 99, 60, 63, 82, 96, 60, 82, & 60, 96, 60, 88, 96,101, 60, 84, 88, 60, 61, 93, & 60, 82, 93, 96,101, 93, 95, 97, 91, 96, 97, 60, & 82, 91, 96,101, 77, 82, 96,101, 52, 81, 96, 99, & 102, 51, 81, 96, 99,102, 46, 84, 96,102, 41, 81, & 96,102, 42, 62, 96,102, 40, 81, 96,102, 50, 82, & 84, 96,102, 43, 84, 96,102, 47, 82, 96,102, 60, & 96,102, 57, 71, 82, 91, 93, 96,102, 54, 60, 63, & 77, 88, 96,102, 56, 81, 96, 99,102, 55, 63, 82, & 88, 96,101,102, 53, 60, 77, 82, 96,102, 45, 82, & 101,102, 62, 65, 66, 81, 84, 86, 87, 89, 95, 96, & 99, 48, 95,102 /) INTEGER, PARAMETER, DIMENSION(1288) :: IROW_STOICM = (/& IROW_STOICM_0, IROW_STOICM_1, IROW_STOICM_2, IROW_STOICM_3 /) 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, & 21, 22, 22, 22, 23, 23, 23, 24, 24, 24, 25, 25, & 25, 26, 26, 26, 27, 27, 27, 28, 28, 28, 29, 29, & 29, 29, 30, 30, 30, 31, 31, 31, 31, 32, 32, 32, & 33, 33, 33, 34, 34, 34, 35, 35, 35, 36, 36, 37, & 37, 37, 37, 37, 38, 38, 38, 38, 38, 38, 39, 39, & 39, 39, 40, 40, 40, 41, 41, 41, 42, 42, 42, 42, & 42, 43, 43, 43, 44, 44, 44, 44, 44, 45, 45, 45, & 46, 46, 46, 47, 47, 47, 47, 47, 48, 48, 48, 48, & 48, 48, 49, 49, 49, 50, 50, 50, 50, 50, 50, 50, & 50, 50, 50, 50, 50, 51, 51, 51, 52, 52, 52, 52, & 52, 52, 52, 53, 53, 53, 53, 53, 53, 54, 54, 54, & 54, 54, 54, 55, 55, 55, 55, 55, 55, 55, 55, 56, & 56, 56, 56, 56, 56, 57, 57, 57, 58, 58, 58, 58, & 58, 58, 58, 58, 58, 58, 59, 59, 59, 59, 59, 59, & 59, 60, 60, 60, 60, 60, 60, 60, 60, 61, 61, 61, & 62, 62, 62, 62, 62, 62, 63, 63, 63, 64, 64, 64, & 64, 64, 64, 64, 64, 64, 65, 65, 65, 65, 65, 66, & 66, 66, 66, 66, 67, 67, 67, 67, 67, 67, 67, 67, & 68, 68, 68, 68, 68, 69, 69, 69, 69, 70, 70, 70, & 71, 71, 71, 72, 72, 72, 73, 73, 73, 74, 74, 74, & 75, 75, 75, 76, 76, 76, 77, 77, 77, 78, 78, 78, & 79, 80, 81, 81, 81, 81, 82, 82, 82, 82, 82, 83, & 83, 83, 83, 83, 84, 84, 84, 84, 85, 85, 85, 85, & 86, 86, 86, 87, 87, 87, 88, 88, 88, 88, 88, 89, & 89, 89, 89, 89, 89, 89, 90, 90, 90, 91, 91, 91 /) INTEGER, PARAMETER, DIMENSION(360) :: ICOL_STOICM_1 = (/ & 92, 92, 92, 92, 93, 93, 93, 93, 94, 94, 94, 95, & 95, 95, 96, 96, 96, 97, 97, 97, 98, 98, 98, 99, & 99, 99,100,100,100,101,101,101,102,102,102,103, & 103,103,104,104,104,105,105,105,106,106,106,106, & 106,107,107,107,107,108,108,108,108,108,108,108, & 108,108,108,108,109,109,109,109,109,109,109,109, & 110,110,110,110,110,110,110,111,111,111,111,111, & 111,111,112,112,112,112,112,112,112,112,112,113, & 113,113,113,113,113,114,114,114,114,114,114,114, & 114,114,115,115,115,115,115,115,115,115,116,116, & 116,116,116,116,116,116,117,117,117,117,117,117, & 118,118,118,118,118,118,118,118,119,119,119,119, & 119,119,119,120,120,120,120,120,121,121,121,121, & 121,121,121,121,121,122,122,122,122,122,122,122, & 123,123,123,124,124,124,125,125,125,126,126,127, & 127,127,128,128,128,129,129,129,130,130,130,130, & 130,130,131,131,131,131,132,132,132,132,132,133, & 133,133,133,134,134,135,135,135,136,136,136,136, & 136,136,136,136,137,137,137,137,137,137,138,138, & 138,138,138,139,139,139,139,139,140,140,140,141, & 141,141,142,142,142,142,143,143,143,143,144,144, & 144,144,144,145,145,145,145,146,146,146,147,147, & 147,147,148,148,148,148,149,149,149,149,149,150, & 150,150,150,150,150,151,151,151,152,152,152,152, & 153,153,153,153,153,153,153,153,153,153,154,154, & 154,154,154,154,154,154,155,155,155,155,155,155, & 155,155,156,156,156,157,157,157,158,158,158,158, & 159,159,159,159,159,160,160,160,160,161,161,161, & 161,161,162,162,162,162,163,163,163,164,164,164, & 165,165,165,166,166,167,167,167,168,168,168,169 /) INTEGER, PARAMETER, DIMENSION(360) :: ICOL_STOICM_2 = (/ & 169,169,169,170,170,170,170,171,171,171,171,172, & 172,172,172,173,173,173,173,174,174,174,175,175, & 175,176,176,177,177,177,177,178,178,178,179,179, & 179,179,179,180,180,180,180,181,181,181,182,182, & 182,183,183,183,183,184,184,184,184,185,185,185, & 186,186,186,186,186,187,187,187,187,187,187,187, & 187,188,188,189,189,189,190,190,190,191,191,191, & 191,191,191,191,191,191,191,191,192,192,192,192, & 192,192,193,193,193,194,194,194,194,194,194,194, & 194,194,195,195,195,195,195,195,196,196,196,196, & 196,196,196,196,196,197,197,197,197,197,197,198, & 198,198,198,198,198,198,199,199,199,199,199,199, & 199,199,200,200,200,200,200,201,201,201,201,201, & 201,201,201,202,202,202,202,202,202,202,203,203, & 203,203,203,203,204,204,204,204,204,204,204,204, & 204,205,205,205,205,205,205,206,206,206,207,207, & 207,208,208,208,209,209,209,210,210,210,211,211, & 211,212,212,212,213,213,213,214,214,214,214,215, & 215,215,215,216,216,216,216,217,217,217,217,218, & 218,218,218,219,219,219,219,220,220,220,221,221, & 221,221,221,222,222,222,223,223,223,223,224,224, & 224,224,224,225,225,225,226,226,226,226,227,227, & 228,229,230,231,232,233,234,235,236,237,238,239, & 240,241,242,243,244,245,246,246,246,247,247,247, & 248,248,248,249,249,249,250,250,250,251,251,251, & 252,252,252,253,253,253,254,254,254,255,255,255, & 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,266, & 266,266,267,267,268,268,268,269,269,269,270,270 /) INTEGER, PARAMETER, DIMENSION(208) :: ICOL_STOICM_3 = (/ & 270,271,271,271,272,272,272,273,273,273,274,274, & 275,275,275,276,276,277,277,277,277,278,278,278, & 279,279,280,280,280,281,281,281,282,282,282,283, & 283,283,284,284,285,285,285,286,286,286,286,287, & 287,287,288,288,288,288,289,289,290,290,290,290, & 290,291,291,291,291,292,292,292,293,293,293,294, & 294,294,294,294,295,295,296,296,296,296,297,297, & 298,298,299,299,299,299,300,300,300,301,301,301, & 302,302,302,302,302,303,303,303,304,304,304,305, & 305,305,305,305,306,306,306,306,307,307,307,307, & 307,308,308,308,308,308,309,309,309,309,310,310, & 310,310,311,311,311,311,312,312,312,312,313,313, & 313,313,313,314,314,314,314,315,315,315,315,316, & 316,316,317,317,317,317,317,317,317,318,318,318, & 318,318,318,318,319,319,319,319,319,320,320,320, & 320,320,320,320,321,321,321,321,321,321,322,322, & 322,322,323,323,323,323,323,323,323,323,323,323, & 323,324,324,324 /) INTEGER, PARAMETER, DIMENSION(1288) :: ICOL_STOICM = (/& ICOL_STOICM_0, ICOL_STOICM_1, ICOL_STOICM_2, ICOL_STOICM_3 /) 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, 1.000000e+00_dp, -2.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, 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, 5.000000e-02_dp, 5.000000e-02_dp, -1.000000e+00_dp, & 5.000000e-02_dp, 9.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, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -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, 3.200000e-01_dp, 5.000000e-02_dp, 1.800000e-01_dp, 1.300000e-01_dp, & 3.200000e-01_dp, -1.000000e+00_dp, 3.200000e-01_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, & -1.000000e+00_dp, -1.000000e+00_dp, 5.700000e-01_dp, 3.900000e-01_dp, 7.500000e-01_dp, & 3.000000e-01_dp, 2.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 9.600000e-01_dp, & 4.000000e-02_dp, 9.600000e-01_dp, 9.600000e-01_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 9.300000e-01_dp, 7.000000e-02_dp, 9.300000e-01_dp, 9.300000e-01_dp, -1.000000e+00_dp, & 3.400000e-01_dp, 5.600000e-01_dp, 2.200000e-01_dp, -1.000000e+00_dp, 3.400000e-01_dp, & 9.000000e-01_dp, 9.000000e-01_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_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, 6.100000e-01_dp, 2.400000e-01_dp, 8.000000e-02_dp, & 3.300000e-01_dp, -1.000000e+00_dp, 3.500000e-01_dp, 5.300000e-01_dp, 9.200000e-01_dp, & 9.200000e-01_dp, -1.000000e+00_dp, 9.500000e-01_dp, 5.000000e-02_dp, 9.500000e-01_dp, & -1.000000e+00_dp, 5.000000e-02_dp, 1.950000e+00_dp, -1.000000e+00_dp, 7.200000e-01_dp, & -1.000000e+00_dp, 2.800000e-01_dp, 2.800000e-01_dp, 2.800000e-01_dp, 1.000000e+00_dp, & 7.200000e-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, 6.000000e-01_dp, & -1.000000e+00_dp, 1.000000e-01_dp, 3.000000e-01_dp, 3.000000e-01_dp, 3.000000e-01_dp, & 1.900000e+00_dp, 6.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, -1.000000e+00_dp, 8.500000e-01_dp, -1.000000e+00_dp, & 1.500000e-01_dp, 1.000000e-01_dp, 5.000000e-02_dp, 8.000000e-01_dp, 1.150000e+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 /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_2 = (/ & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 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.600000e-01_dp, 9.000000e-02_dp, 1.250000e+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, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 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.000000e+00_dp, 1.600000e-01_dp, 3.000000e-02_dp, & 9.000000e-02_dp, 7.000000e-02_dp, 7.500000e-01_dp, 1.600000e-01_dp, -1.000000e+00_dp, & 1.600000e-01_dp, 3.500000e-01_dp, -9.100000e-01_dp, 6.400000e-01_dp, -1.000000e+00_dp, & 5.400000e-01_dp, 9.500000e-01_dp, 3.800000e-01_dp, 1.500000e-01_dp, -1.000000e+00_dp, & 5.000000e-01_dp, 1.000000e+00_dp, -1.000000e+00_dp, 2.000000e-01_dp, 5.000000e-01_dp, & 5.000000e-01_dp, -1.000000e+00_dp, 3.000000e-01_dp, 3.000000e-01_dp, -1.000000e+00_dp, & 7.500000e-01_dp, 5.000000e-01_dp, 2.500000e-01_dp, -1.000000e+00_dp, 5.000000e-01_dp, & 5.000000e-01_dp, 7.000000e-02_dp, 6.000000e-02_dp, 1.100000e+00_dp, 2.500000e-01_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_3 = (/ & 1.400000e-01_dp, -1.000000e+00_dp, 2.000000e-01_dp, -1.000000e+00_dp, 9.200000e-01_dp, & -1.000000e+00_dp, 5.000000e-01_dp, 1.130000e+00_dp, 2.500000e-01_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 3.300000e-01_dp, 1.300000e-01_dp, 1.800000e-01_dp, -1.000000e+00_dp, & 9.500000e-01_dp, 2.900000e-01_dp, 2.500000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 5.000000e-01_dp, 5.000000e-01_dp, -1.000000e+00_dp, 2.500000e-01_dp, 7.500000e-01_dp, & -1.000000e+00_dp, 5.000000e-01_dp, 1.000000e+00_dp, 3.600000e-01_dp, -1.000000e+00_dp, & 8.900000e-01_dp, 1.400000e-01_dp, 2.500000e-01_dp, -1.000000e+00_dp, 6.400000e-01_dp, & 3.600000e-01_dp, 1.500000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, 8.500000e-01_dp, & -1.000000e+00_dp, 1.150000e+00_dp, -1.000000e+00_dp, 2.500000e-01_dp, 1.250000e+00_dp, & 2.500000e-01_dp, -1.000000e+00_dp, 7.500000e-01_dp, 1.000000e+00_dp, 2.500000e-01_dp, & -1.000000e+00_dp, 2.500000e-01_dp, 1.250000e+00_dp, 5.000000e-01_dp, -1.000000e+00_dp, & 5.000000e-01_dp, 1.000000e+00_dp, 7.500000e-01_dp, -1.000000e+00_dp, 7.500000e-01_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 4.250000e-01_dp, -1.000000e+00_dp, 2.500000e-01_dp, & 8.300000e-01_dp, 5.000000e-02_dp, 3.000000e-02_dp, -1.000000e+00_dp, 4.500000e-01_dp, & 5.750000e-01_dp, -1.000000e+00_dp, 2.500000e-01_dp, 1.250000e+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, 2.000000e+00_dp, & -2.000000e+00_dp, 2.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, 4.100000e-01_dp, 4.400000e-01_dp, & -1.000000e+00_dp, 1.500000e-01_dp, -1.000000e+00_dp, 4.400000e-01_dp, 7.000000e-01_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 3.000000e-01_dp, 7.100000e-01_dp, 2.900000e-01_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 2.900000e-01_dp, 7.000000e-01_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 3.000000e-01_dp, -1.000000e+00_dp, 3.000000e-01_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 4.200000e-01_dp, -1.000000e+00_dp, 5.350000e-01_dp, & 5.000000e-01_dp, 3.050000e-01_dp, 3.000000e-01_dp, -1.000000e+00_dp, 1.350000e-01_dp, & -1.000000e+00_dp, 5.900000e-01_dp, 2.230000e+00_dp, 2.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, 4.000000e-01_dp, -1.000000e+00_dp, 8.000000e-01_dp, 2.000000e-01_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 2.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, & -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.300000e-01_dp, -1.000000e+00_dp, 5.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, -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.000000e+00_dp, 1.500000e-01_dp, -1.000000e+00_dp, & 5.000000e-02_dp, 7.000000e-02_dp, 9.000000e-01_dp, 3.870000e-01_dp, 1.590000e-01_dp, & 6.000000e-02_dp, -9.000000e-01_dp, 2.700000e-01_dp, 5.000000e-02_dp, 8.000000e-01_dp, & 4.000000e-02_dp, 8.200000e-01_dp, -1.000000e+00_dp, 6.000000e-02_dp, -8.000000e-01_dp, & 8.000000e-02_dp, 1.600000e-01_dp, 2.000000e-01_dp, 7.000000e-01_dp, 8.000000e-01_dp, & -1.000000e+00_dp, 2.750000e-01_dp, -8.000000e-01_dp, 2.150000e-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, 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, 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, & 5.000000e-01_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, 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, 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, & 5.000000e-01_dp, -5.000000e-01_dp, -1.000000e+00_dp, 5.090000e-01_dp, 4.910000e-01_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_5 = (/ & -4.910000e-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, 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, & 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, 4.400000e-01_dp, 1.500000e-01_dp, 4.100000e-01_dp, -1.000000e+00_dp, & 4.000000e-01_dp, 2.800000e-01_dp, -1.000000e+00_dp, 2.000000e-01_dp, 1.200000e-01_dp, & 6.000000e-01_dp, -7.000000e-01_dp, 1.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, 3.200000e-01_dp, 5.000000e-02_dp, 1.800000e-01_dp, 1.300000e-01_dp, & 3.200000e-01_dp, -1.000000e+00_dp, 3.200000e-01_dp, 1.900000e-01_dp, 1.180000e+00_dp, & 2.700000e-01_dp, -1.000000e+00_dp, -1.000000e+00_dp, 2.000000e-01_dp, 8.000000e-01_dp, & 1.000000e+00_dp, 8.000000e-01_dp, -8.000000e-01_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.360000e-01_dp, 1.270000e-01_dp, 6.900000e-01_dp, 2.880000e-01_dp, & -1.000000e+00_dp, 4.020000e-01_dp, 1.000000e+00_dp, 8.640000e-01_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 7.500000e-01_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 6.500000e-01_dp, 2.600000e-01_dp, 3.600000e-01_dp, -1.000000e+00_dp, & 4.000000e-01_dp, 5.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, 7.200000e-01_dp, -1.000000e+00_dp, 2.800000e-01_dp, 2.800000e-01_dp, & 1.000000e+00_dp, 2.800000e-01_dp, -2.800000e-01_dp, 8.300000e-01_dp, -1.000000e+00_dp, & 8.300000e-01_dp, 1.700000e-01_dp, 1.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, 5.700000e-01_dp, 3.900000e-01_dp, 7.500000e-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, 5.000000e-01_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, 8.500000e-01_dp, -1.000000e+00_dp, 1.500000e-01_dp, & 1.000000e-01_dp, 5.000000e-02_dp, 1.000000e+00_dp, 8.000000e-01_dp, 1.500000e-01_dp /) REAL(kind=dp), PARAMETER, DIMENSION(150) :: STOICM_6 = (/ & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 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, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 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, -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 /) 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, 3.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & -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, 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, 4.000000e-01_dp, 6.000000e-01_dp, 4.000000e-01_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, & 2.000000e+00_dp, 8.500000e-01_dp, -1.000000e+00_dp, 1.500000e-01_dp, 1.500000e-01_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, 2.000000e+00_dp, 1.500000e+00_dp, 5.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, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 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(88) :: STOICM_8 = (/ & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 1.000000e+00_dp, & 1.000000e+00_dp, -1.000000e+00_dp, 1.000000e+00_dp, 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.730000e-01_dp, & 6.270000e-01_dp, 2.590000e-01_dp, 3.680000e-01_dp, 1.000000e+00_dp, 1.000000e+00_dp, & -1.000000e+00_dp, 6.700000e-01_dp, 2.600000e-01_dp, 3.600000e-01_dp, 5.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, 7.000000e-01_dp, 3.000000e-01_dp, & 3.000000e-01_dp, 3.000000e-01_dp, 7.000000e-01_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, 3.200000e-01_dp, & 5.000000e-02_dp, 1.800000e-01_dp, 1.300000e-01_dp, 3.200000e-01_dp, -1.000000e+00_dp, & 3.200000e-01_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 /) REAL(kind=dp), PARAMETER, DIMENSION(1288) :: STOICM = (/& STOICM_0, STOICM_1, STOICM_2, STOICM_3, STOICM_4, & STOICM_5, STOICM_6, STOICM_7, STOICM_8 /) END MODULE gckpp_StoichiomSP