! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! Sparse Jacobian 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_JacobianSP.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_JacobianSP PUBLIC SAVE ! Sparse Jacobian Data INTEGER, PARAMETER, DIMENSION(360) :: LU_IROW_0 = (/ & 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 6, 6, & 7, 7, 8, 8, 9, 9, 10, 10, 11, 11, 12, 12, & 12, 13, 13, 13, 14, 14, 14, 15, 15, 15, 16, 16, & 16, 17, 17, 17, 17, 17, 18, 18, 19, 19, 19, 20, & 20, 20, 21, 21, 21, 21, 21, 22, 22, 23, 23, 23, & 24, 24, 24, 25, 25, 25, 25, 25, 26, 26, 27, 27, & 27, 28, 28, 28, 29, 29, 29, 29, 29, 30, 30, 31, & 31, 31, 32, 32, 33, 33, 34, 34, 35, 35, 36, 36, & 37, 37, 38, 38, 39, 39, 40, 40, 41, 41, 42, 42, & 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, & 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, & 42, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, & 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 44, & 44, 45, 45, 45, 46, 46, 46, 47, 47, 47, 48, 48, & 48, 48, 49, 49, 49, 50, 50, 50, 50, 50, 51, 51, & 51, 52, 52, 52, 53, 53, 53, 53, 54, 54, 54, 54, & 55, 55, 55, 56, 56, 56, 56, 57, 57, 57, 57, 58, & 58, 58, 58, 59, 59, 59, 59, 60, 60, 60, 61, 61, & 61, 61, 62, 62, 62, 62, 63, 63, 63, 63, 64, 64, & 64, 64, 65, 65, 65, 65, 65, 66, 66, 66, 66, 66, & 67, 67, 67, 67, 68, 68, 68, 68, 68, 68, 68, 69, & 69, 69, 69, 70, 70, 70, 70, 71, 71, 71, 71, 72, & 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 73, & 73, 73, 73, 74, 74, 74, 74, 75, 75, 75, 75, 75, & 75, 76, 76, 76, 76, 77, 77, 77, 77, 77, 77, 77, & 78, 78, 78, 78, 78, 78, 79, 79, 79, 79, 79, 80, & 80, 80, 80, 80, 81, 81, 81, 81, 82, 82, 82, 82, & 82, 83, 83, 83, 83, 83, 84, 84, 84, 84, 84, 84, & 84, 85, 85, 85, 85, 85, 85, 85, 85, 85, 85, 86, & 86, 86, 86, 86, 86, 87, 87, 87, 87, 87, 87, 87 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_IROW_1 = (/ & 87, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, 88, & 88, 88, 88, 88, 88, 89, 89, 89, 89, 89, 89, 89, & 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, 89, & 89, 89, 89, 89, 89, 90, 90, 90, 90, 90, 90, 90, & 90, 90, 90, 90, 91, 91, 91, 91, 91, 91, 91, 91, & 91, 91, 91, 91, 91, 91, 91, 91, 91, 91, 91, 91, & 91, 91, 91, 91, 91, 91, 91, 91, 91, 91, 91, 91, & 91, 91, 91, 92, 92, 92, 92, 92, 92, 92, 92, 93, & 93, 93, 93, 93, 93, 93, 93, 93, 94, 94, 94, 94, & 94, 94, 94, 94, 94, 94, 95, 95, 95, 95, 95, 95, & 95, 95, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, & 96, 96, 97, 97, 97, 97, 97, 97, 97, 98, 98, 98, & 98, 98, 98, 98, 98, 98, 98, 99, 99, 99, 99, 99, & 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, & 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, 99, & 100,100,100,100,100,100,100,100,100,100,100,100, & 100,100,100,100,100,100,100,101,101,101,101,101, & 101,101,101,102,102,102,102,102,102,102,102,103, & 103,103,103,103,103,103,104,104,104,104,104,104, & 104,104,104,104,104,105,105,105,105,105,105,105, & 105,106,106,106,106,106,106,106,106,107,107,107, & 107,107,107,107,107,107,107,108,108,108,108,108, & 108,108,108,108,108,108,108,109,109,109,109,109, & 109,109,109,109,109,109,110,110,110,110,110,110, & 110,110,110,110,110,111,111,111,111,111,111,111, & 111,112,112,112,112,112,112,112,112,112,112,112, & 112,112,112,113,113,113,113,113,113,113,113,113, & 113,113,113,114,114,114,114,114,114,114,114,114, & 114,114,114,114,114,114,114,114,114,114,114,114, & 114,114,114,114,114,114,114,114,114,114,114,114 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_IROW_2 = (/ & 114,114,114,114,114,114,114,114,114,114,114,114, & 114,114,114,114,114,114,114,114,114,114,114,114, & 114,114,115,115,115,115,115,115,115,115,115,115, & 115,115,115,115,115,115,115,115,115,115,115,115, & 115,115,115,116,116,116,116,116,116,116,116,116, & 116,116,116,116,117,117,117,117,117,117,117,117, & 117,117,118,118,118,118,118,118,118,118,118,118, & 118,118,118,118,118,118,118,118,118,118,118,118, & 118,118,118,118,118,118,118,118,118,118,118,118, & 118,118,119,119,119,119,119,119,119,119,119,119, & 119,119,119,119,119,119,119,119,119,119,119,119, & 119,119,119,119,119,119,119,119,119,119,119,120, & 120,120,120,120,120,120,120,120,120,120,120,121, & 121,121,121,121,121,121,121,121,121,121,121,122, & 122,122,122,122,122,122,122,122,122,122,122,122, & 122,123,123,123,123,123,123,123,123,123,123,124, & 124,124,124,124,124,124,124,124,124,124,124,124, & 125,125,125,125,125,125,125,125,125,125,126,126, & 126,126,126,126,126,126,126,126,126,126,126,126, & 126,126,126,127,127,127,127,127,127,127,127,127, & 127,127,127,127,127,127,127,127,127,127,127,127, & 127,127,127,127,127,127,127,127,127,127,127,127, & 127,127,127,127,127,127,128,128,128,128,128,128, & 128,128,128,128,128,128,128,128,128,128,128,128, & 128,128,128,128,128,128,128,128,128,128,128,128, & 128,128,128,128,128,128,128,128,128,128,128,128, & 128,128,128,128,128,128,128,128,128,128,128,128, & 128,128,128,128,128,128,128,128,128,129,129,129, & 129,129,129,129,129,129,129,129,129,129,129,129, & 129,129,129,129,129,129,129,129,129,129,129,129 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_IROW_3 = (/ & 129,129,129,129,129,129,129,129,129,129,129,129, & 129,129,129,129,129,129,129,129,130,130,130,130, & 130,130,130,130,130,130,130,130,130,130,130,130, & 130,130,130,130,130,130,130,130,130,130,130,130, & 130,130,130,130,130,130,130,130,130,130,130,130, & 130,130,130,130,130,130,130,131,131,131,131,131, & 131,131,131,131,131,131,131,131,131,131,131,131, & 131,131,131,131,131,131,131,131,131,131,131,131, & 132,132,132,132,132,132,132,132,132,132,132,132, & 132,132,132,132,132,132,132,132,132,132,132,132, & 132,132,132,132,132,132,132,132,132,132,132,132, & 132,132,132,132,132,132,132,132,132,132,132,132, & 132,132,132,132,132,132,132,132,132,132,132,132, & 132,132,132,132,132,132,132,132,132,132,132,132, & 132,132,132,132,132,132,132,132,132,132,132,132, & 132,132,132,132,132,133,133,133,133,133,133,133, & 133,133,133,133,133,133,133,133,133,133,133,133, & 133,133,133,133,133,133,133,133,133,133,133,133, & 133,133,133,133,133,133,133,133,133,133,133,133, & 133,133,133,133,133,133,134,134,134,134,134,134, & 134,134,134,134,134,134,134,134,134,134,134,134, & 134,134,134,134,134,134,134,134,134,134,134,134, & 134,134,134,134,134,134,134,134,134,134,134,134, & 134,134,135,135,135,135,135,135,135,135,135,135, & 135,135,135,135,135,135,135,135,135,135,135,135, & 135,135,135,136,136,136,136,136,136,136,136,136, & 136,136,136,136,136,136,136,136,136,136,136,136, & 136,136,136,136,136,136,136,136,136,136,136,136, & 136,136,136,136,136,136,136,136,136,136,136,136, & 136,136,136,136,136,136,136,136,136,136,136,136 /) INTEGER, PARAMETER, DIMENSION(41) :: LU_IROW_4 = (/ & 136,136,136,136,136,136,136,136,136,136,136,136, & 136,136,136,136,136,136,136,136,137,137,137,137, & 137,137,137,137,137,137,137,137,137,137,137,137, & 137,137,137,137,137 /) INTEGER, PARAMETER, DIMENSION(1481) :: LU_IROW = (/& LU_IROW_0, LU_IROW_1, LU_IROW_2, LU_IROW_3, LU_IROW_4 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_ICOL_0 = (/ & 1,114, 2, 46, 3, 99, 4, 55, 5,128, 6,135, & 7, 49, 8, 79, 9, 50, 10, 45, 11,137, 12, 48, & 132, 13, 60,132, 14, 81,132, 15, 17,136, 16, 17, & 134, 17, 18,132,134,136, 18,132, 19, 21,136, 20, & 21,134, 21, 22,132,134,136, 22,132, 23, 25,136, & 24, 25,134, 25, 26,132,134,136, 26,132, 27, 29, & 136, 28, 29,134, 29, 30,132,134,136, 30,132, 31, & 93,132, 32,125, 33, 82, 34, 74, 35, 47, 36,104, & 37, 93, 38, 96, 39, 72, 40, 89, 41, 65, 42, 65, & 68, 70, 72, 75, 78, 81, 84, 87, 89, 91, 92, 95, & 97,100,102,119,121,122,123,130,132,133,134,135, & 136, 43, 45, 46, 47, 49, 50, 55, 65, 72, 74, 79, & 82, 89, 93, 96, 99,104,114,125,128,135,137, 44, & 132, 45,113,128, 46,132,136, 47, 74,132, 48, 60, & 129,132, 49,128,130, 50, 92, 95,128,134, 51,129, & 132, 52,129,132, 53,128,132,134, 54,130,132,136, & 55,128,129, 56, 94,132,136, 57, 98,132,136, 58, & 113,132,136, 59,128,132,136, 60,129,132, 61,112, & 132,136, 62,132,133,136, 63,108,132,136, 64,122, & 132,136, 65, 75, 95,128,134, 66, 87,126,132,136, & 67,110,132,136, 68,103,125,132,134,135,136, 69, & 109,132,136, 70,105,132,136, 71,117,132,136, 72, & 84, 92, 93,100,103,119,124,132,134,135,136, 73, & 101,132,136, 74,123,132,136, 65, 75, 95,128,134, & 136, 76,111,132,136, 76, 77,111,130,132,134,136, & 78,111,130,132,134,135, 79,122,128,132,135, 80, & 103,107,132,136, 81,129,132,135, 82,123,132,134, & 135, 83,123,130,134,136, 47, 74, 84,123,132,134, & 136, 57, 85, 98,126,130,132,133,134,136,137, 81, & 86,120,129,132,135, 78, 87,111,130,132,134,135 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_ICOL_1 = (/ & 136, 88, 94, 98,101,102,105,106,109,110,111,116, & 117,123,126,132,133, 89, 92, 94, 97, 98,101,102, & 105,106,108,109,110,111,112,116,117,119,123,126, & 130,132,133,134,136, 83, 84, 90,100,111,123,129, & 130,132,134,136, 70, 76, 77, 81, 84, 85, 86, 90, & 91, 96, 97, 98,100,104,105,111,114,115,118,119, & 120,121,122,123,124,126,127,129,130,132,133,134, & 135,136,137, 50, 92, 95,104,128,132,134,136, 82, & 93,101,107,123,132,134,135,136, 44, 56, 94,126, & 130,132,133,134,136,137, 79, 95,122,128,132,134, & 135,136, 67, 96, 97,103,110,125,130,132,133,134, & 135,136, 97,121,129,130,133,134,136, 44, 57, 98, & 126,130,132,133,134,136,137, 51, 52, 55, 60, 83, & 87, 90, 99,100,102,103,107,111,114,115,117,118, & 121,123,127,128,129,130,131,132,133,134,135,136, & 73, 76, 83, 84, 93,100,101,103,107,111,123,124, & 125,130,132,133,134,135,136, 73,101,120,130,132, & 133,134,136, 64,102,122,130,132,133,134,136, 80, & 103,107,125,132,134,136, 82,103,104,105,107,123, & 125,132,134,135,136, 70,105,121,130,132,133,134, & 136,106,129,130,131,132,133,134,136, 80, 82,103, & 107,123,125,132,134,135,136, 63, 85, 98,108,122, & 126,130,132,133,134,136,137, 69, 86,109,120,129, & 130,132,133,134,135,136, 67, 86,110,120,129,130, & 132,133,134,135,136, 76,111,124,130,132,133,134, & 136, 52, 61,112,113,126,127,129,130,131,132,133, & 134,136,137, 45, 58,113,127,128,129,130,131,132, & 133,134,136, 60, 62, 63, 64, 68, 69, 70, 73, 74, & 79, 81, 82, 84, 86, 90, 92, 93, 94, 95, 96, 97, & 98,100,101,102,103,104,105,106,107,108,109,110 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_ICOL_2 = (/ & 111,112,113,114,115,116,117,119,120,121,122,123, & 124,125,126,127,128,129,130,131,132,133,134,135, & 136,137, 58, 61, 69, 86,106,109,110,112,113,115, & 116,118,120,126,127,128,129,130,131,132,133,134, & 135,136,137, 80,103,107,116,123,125,130,132,133, & 134,135,136,137, 71, 81,117,129,130,132,133,134, & 135,136, 63, 73, 76, 83, 84, 92, 93, 95, 96, 97, & 101,103,104,105,107,108,110,111,118,119,120,121, & 122,123,124,125,126,128,129,130,132,133,134,135, & 136,137, 69, 70, 76, 77, 83, 84, 92, 95,102,103, & 104,105,107,108,109,111,119,120,121,122,123,124, & 125,126,128,129,130,132,133,134,135,136,137, 74, & 81,117,120,123,129,130,132,133,134,135,136, 74, & 81,117,121,123,129,130,132,133,134,135,136, 64, & 79,120,121,122,123,128,129,130,132,133,134,135, & 136, 74, 81,123,129,130,132,133,134,135,136, 74, & 82,117,123,124,125,129,130,132,133,134,135,136, & 117,123,125,129,130,132,133,134,135,136, 51, 66, & 87,111,116,123,124,125,126,129,130,132,133,134, & 135,136,137, 56, 66, 67, 71, 80, 87, 88, 94, 97, & 98,101,102,103,105,106,107,109,110,111,116,117, & 120,121,122,123,124,125,126,127,128,129,130,131, & 132,133,134,135,136,137, 45, 49, 50, 53, 55, 59, & 65, 67, 68, 71, 75, 79, 80, 82, 83, 84, 87, 92, & 93, 94, 95, 96, 97, 98, 99,100,101,102,103,104, & 105,106,107,108,109,110,111,112,113,114,115,116, & 117,118,119,120,121,122,123,124,125,126,127,128, & 129,130,131,132,133,134,135,136,137, 49, 51, 52, & 55, 59, 60, 65, 75, 81, 86, 90, 92, 93, 95, 99, & 100,101,102,103,104,105,107,111,114,115,116,117 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_ICOL_3 = (/ & 118,119,120,121,122,123,124,125,126,127,128,129, & 130,131,132,133,134,135,136,137, 49, 54, 63, 64, & 65, 73, 75, 85, 93, 94, 95, 97, 98,101,102,105, & 106,107,108,109,110,111,112,113,115,116,117,118, & 119,120,121,122,123,124,125,126,127,128,129,130, & 131,132,133,134,135,136,137, 71, 73, 77, 78,101, & 102,103,105,106,107,111,117,120,121,122,123,124, & 125,126,128,129,130,131,132,133,134,135,136,137, & 44, 46, 47, 48, 51, 52, 53, 54, 56, 57, 58, 59, & 60, 61, 62, 63, 64, 66, 67, 68, 69, 70, 71, 72, & 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, & 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, & 97, 98, 99,100,101,102,103,104,105,106,107,108, & 109,110,111,112,113,114,115,116,117,118,119,120, & 121,122,123,124,125,126,127,128,129,130,131,132, & 133,134,135,136,137, 49, 54, 60, 62, 85, 86, 89, & 92, 94, 95, 96, 97, 98,101,102,103,104,105,106, & 107,108,109,110,111,112,113,115,116,117,118,119, & 120,121,122,123,124,125,126,127,128,129,130,131, & 132,133,134,135,136,137, 53, 55, 75, 83, 84, 87, & 92, 94, 95, 97, 98,101,102,103,104,105,106,107, & 108,109,110,111,112,113,116,117,120,121,122,123, & 124,125,126,127,128,129,130,131,132,133,134,135, & 136,137, 55, 78, 79, 81, 82, 86,111,113,120,121, & 122,123,124,125,127,128,129,130,131,132,133,134, & 135,136,137, 46, 48, 56, 57, 58, 59, 60, 61, 62, & 66, 67, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, & 79, 80, 81, 82, 83, 84, 86, 87, 88, 90, 91, 92, & 93, 94, 95, 96, 97, 98,100,101,102,103,104,105, & 106,107,108,109,110,111,112,113,114,115,116,117 /) INTEGER, PARAMETER, DIMENSION(41) :: LU_ICOL_4 = (/ & 118,119,120,121,122,123,124,125,126,127,128,129, & 130,131,132,133,134,135,136,137, 71, 97,106,108, & 116,117,121,122,123,125,126,128,129,130,131,132, & 133,134,135,136,137 /) INTEGER, PARAMETER, DIMENSION(1481) :: LU_ICOL = (/& LU_ICOL_0, LU_ICOL_1, LU_ICOL_2, LU_ICOL_3, LU_ICOL_4 /) INTEGER, PARAMETER, DIMENSION(138) :: LU_CROW = (/ & 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, & 26, 29, 32, 35, 38, 43, 45, 48, 51, 56, 58, 61, & 64, 69, 71, 74, 77, 82, 84, 87, 89, 91, 93, 95, & 97, 99,101,103,105,107,134,156,158,161,164,167, & 171,174,179,182,185,189,193,196,200,204,208,212, & 215,219,223,227,231,236,241,245,252,256,260,264, & 276,280,284,290,294,301,307,312,317,321,326,331, & 338,348,354,362,378,402,413,448,456,465,475,483, & 495,502,512,541,560,568,576,583,594,602,610,620, & 632,643,654,662,676,688,747,772,785,795,831,864, & 876,888,902,912,925,935,952,991,1054,1101,1148,1177, & 1266,1315,1359,1384,1461,1482 /) INTEGER, PARAMETER, DIMENSION(138) :: LU_DIAG = (/ & 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, & 26, 29, 32, 35, 38, 43, 45, 48, 51, 56, 58, 61, & 64, 69, 71, 74, 77, 82, 84, 87, 89, 91, 93, 95, & 97, 99,101,103,105,107,134,156,158,161,164,167, & 171,174,179,182,185,189,193,196,200,204,208,212, & 215,219,223,227,231,236,241,245,252,256,260,264, & 276,280,285,290,295,301,307,312,317,321,326,333, & 339,349,355,362,378,404,421,449,457,467,476,484, & 495,504,519,546,561,569,577,585,595,602,613,623, & 634,645,655,664,678,724,756,775,787,813,847,867, & 879,892,904,916,927,943,980,1044,1092,1140,1170,1260, & 1310,1355,1381,1459,1481,1482 /) END MODULE gckpp_JacobianSP