! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! 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 : Wed Sep 15 15:20:39 2010 ! Working directory : /mnt/lstr04/srv/home/c/ccarouge/KPP/geoschem_kppfiles/v8-03-02/standard ! 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, 11, 12, & 12, 12, 13, 13, 13, 13, 13, 13, 13, 13, 14, 14, & 14, 14, 14, 14, 14, 14, 14, 14, 14, 15, 15, 15, & 16, 16, 17, 17, 17, 18, 18, 18, 19, 19, 19, 20, & 20, 20, 20, 21, 21, 21, 22, 22, 22, 23, 23, 23, & 24, 24, 24, 24, 25, 25, 25, 26, 26, 26, 26, 27, & 27, 27, 27, 28, 28, 28, 28, 29, 29, 29, 29, 30, & 30, 30, 30, 31, 31, 31, 31, 32, 32, 32, 32, 33, & 33, 33, 33, 34, 34, 34, 35, 35, 35, 35, 36, 36, & 36, 36, 37, 37, 37, 37, 38, 38, 38, 38, 39, 39, & 39, 39, 40, 40, 40, 40, 41, 41, 41, 41, 42, 42, & 42, 42, 43, 43, 43, 43, 43, 43, 44, 44, 44, 44, & 44, 45, 45, 45, 45, 45, 46, 46, 46, 46, 47, 47, & 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, & 47, 47, 47, 47, 47, 47, 47, 47, 47, 48, 48, 48, & 48, 48, 48, 49, 49, 49, 49, 49, 49, 49, 49, 49, & 49, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, & 50, 50, 50, 51, 51, 51, 51, 51, 52, 52, 52, 52, & 52, 52, 52, 52, 52, 52, 53, 53, 53, 53, 53, 53, & 53, 53, 53, 53, 54, 54, 54, 54, 54, 54, 54, 54, & 54, 54, 54, 54, 54, 54, 54, 54, 54, 54, 54, 54, & 54, 54, 54, 54, 54, 55, 55, 55, 55, 55, 55, 55, & 56, 56, 56, 56, 56, 56, 56, 57, 57, 57, 57, 57, & 57, 58, 58, 58, 58, 58, 58, 58, 58, 58, 58, 59, & 59, 59, 59, 59, 59, 59, 59, 60, 60, 60, 60, 60, & 60, 60, 60, 61, 61, 61, 61, 61, 61, 61, 61, 61, & 61, 61, 61, 61, 61, 61, 61, 61, 61, 61, 61, 62, & 62, 62, 62, 62, 62, 62, 62, 62, 62, 63, 63, 63, & 63, 63, 63, 63, 63, 63, 63, 63, 64, 64, 64, 64 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_IROW_1 = (/ & 64, 64, 64, 64, 65, 65, 65, 65, 65, 65, 65, 65, & 65, 65, 65, 66, 66, 66, 66, 66, 66, 66, 66, 66, & 66, 67, 67, 67, 67, 67, 67, 67, 67, 67, 67, 68, & 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, & 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, & 68, 68, 69, 69, 69, 69, 69, 69, 69, 69, 69, 69, & 69, 69, 69, 69, 69, 69, 69, 69, 69, 69, 69, 69, & 69, 69, 69, 69, 69, 69, 69, 69, 69, 69, 69, 69, & 69, 69, 69, 69, 69, 69, 69, 69, 69, 69, 69, 69, & 69, 69, 69, 69, 70, 70, 70, 70, 70, 70, 70, 70, & 70, 70, 70, 71, 71, 71, 71, 71, 71, 71, 71, 71, & 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, 71, & 71, 71, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, & 72, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, 73, & 73, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, 74, & 74, 74, 74, 74, 74, 74, 74, 75, 75, 75, 75, 75, & 75, 75, 75, 75, 75, 75, 75, 75, 75, 75, 75, 75, & 75, 75, 75, 75, 75, 75, 75, 75, 75, 75, 75, 75, & 76, 76, 76, 76, 76, 76, 76, 76, 76, 76, 76, 76, & 76, 76, 76, 76, 76, 76, 76, 77, 77, 77, 77, 77, & 77, 77, 77, 77, 77, 77, 77, 77, 78, 78, 78, 78, & 78, 78, 78, 78, 78, 78, 78, 78, 78, 78, 78, 79, & 79, 79, 79, 79, 79, 79, 79, 79, 79, 79, 80, 80, & 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 80, 81, & 81, 81, 81, 81, 81, 81, 81, 81, 81, 81, 81, 81, & 81, 81, 81, 81, 81, 81, 82, 82, 82, 82, 82, 82, & 82, 82, 82, 82, 82, 82, 82, 82, 82, 82, 82, 82, & 82, 82, 82, 82, 82, 82, 82, 82, 82, 82, 82, 82, & 82, 82, 82, 82, 82, 82, 82, 82, 82, 82, 82, 82, & 82, 82, 82, 82, 82, 82, 83, 83, 83, 83, 83, 83 /) INTEGER, PARAMETER, DIMENSION(321) :: LU_IROW_2 = (/ & 83, 83, 83, 83, 83, 83, 83, 83, 83, 83, 83, 83, & 83, 83, 83, 83, 83, 83, 83, 83, 83, 83, 83, 83, & 83, 83, 83, 83, 83, 83, 83, 83, 83, 83, 83, 84, & 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, & 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, & 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, & 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, & 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, 84, & 84, 84, 84, 85, 85, 85, 85, 85, 85, 85, 85, 85, & 85, 85, 85, 85, 85, 85, 85, 85, 85, 85, 85, 85, & 85, 85, 85, 85, 85, 85, 85, 85, 85, 85, 85, 85, & 85, 85, 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, & 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, & 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, 86, & 86, 86, 86, 86, 86, 86, 86, 87, 87, 87, 87, 87, & 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, & 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, 87, & 87, 87, 87, 87, 87, 87, 87, 87, 87, 88, 88, 88, & 88, 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, & 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, & 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, & 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, & 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, & 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, 90, & 90, 90, 90, 90, 90, 90, 90, 90, 90 /) INTEGER, PARAMETER, DIMENSION(1041) :: LU_IROW = (/& LU_IROW_0, LU_IROW_1, LU_IROW_2 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_ICOL_0 = (/ & 1, 69, 2, 17, 3, 54, 4, 25, 5, 82, 6, 89, & 7, 21, 8, 44, 9, 18, 10, 73, 11, 20, 90, 12, & 34, 90, 13, 46, 47, 80, 83, 85, 89, 90, 14, 17, & 18, 21, 25, 44, 54, 69, 73, 82, 89, 15, 46, 90, & 16, 90, 17, 84, 90, 18, 81, 82, 19, 82, 88, 20, & 34, 87, 90, 21, 82, 83, 22, 87, 90, 23, 87, 90, & 24, 82, 85, 90, 25, 82, 87, 26, 78, 84, 90, 27, & 83, 84, 90, 28, 84, 86, 90, 29, 82, 84, 90, 30, & 72, 84, 90, 31, 52, 84, 90, 32, 53, 84, 90, 33, & 81, 84, 90, 34, 87, 90, 35, 74, 84, 90, 36, 84, & 88, 90, 37, 70, 84, 90, 38, 63, 84, 90, 39, 67, & 84, 90, 40, 59, 84, 90, 41, 66, 84, 90, 42, 64, & 84, 90, 43, 51, 56, 65, 84, 90, 44, 78, 82, 89, & 90, 45, 57, 79, 84, 90, 46, 87, 89, 90, 40, 41, & 46, 47, 48, 49, 50, 58, 59, 66, 68, 69, 71, 75, & 77, 80, 83, 84, 85, 86, 87, 89, 90, 46, 48, 77, & 87, 89, 90, 32, 49, 53, 72, 73, 83, 84, 85, 86, & 90, 41, 42, 50, 51, 58, 64, 65, 66, 83, 84, 85, & 86, 89, 90, 51, 83, 84, 85, 86, 16, 31, 52, 72, & 73, 83, 84, 85, 86, 90, 16, 32, 53, 72, 73, 83, & 84, 85, 86, 90, 22, 23, 25, 34, 51, 54, 57, 59, & 64, 65, 66, 67, 68, 69, 71, 75, 76, 80, 82, 83, & 84, 85, 86, 87, 90, 55, 73, 83, 84, 85, 86, 90, & 56, 80, 83, 84, 85, 86, 87, 57, 79, 83, 84, 85, & 86, 45, 57, 58, 79, 83, 84, 85, 86, 89, 90, 40, & 59, 80, 83, 84, 85, 86, 90, 60, 76, 83, 84, 85, & 86, 87, 90, 40, 41, 44, 58, 59, 61, 62, 65, 66, & 70, 78, 79, 80, 82, 83, 84, 85, 86, 89, 90, 49, & 53, 62, 72, 73, 83, 84, 85, 86, 90, 38, 48, 63, & 77, 83, 84, 85, 86, 87, 89, 90, 42, 64, 77, 83 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_ICOL_1 = (/ & 84, 85, 86, 90, 43, 51, 56, 65, 80, 83, 84, 85, & 86, 87, 90, 41, 58, 66, 79, 83, 84, 85, 86, 89, & 90, 39, 46, 67, 83, 84, 85, 86, 87, 89, 90, 30, & 31, 38, 39, 42, 43, 51, 52, 55, 56, 63, 64, 65, & 67, 68, 70, 72, 73, 77, 80, 83, 84, 85, 86, 87, & 89, 90, 26, 27, 28, 34, 36, 37, 40, 42, 44, 45, & 46, 48, 50, 51, 52, 53, 55, 56, 57, 58, 59, 60, & 61, 62, 63, 64, 65, 66, 67, 69, 70, 71, 72, 73, & 74, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, & 87, 88, 89, 90, 37, 48, 70, 77, 83, 84, 85, 86, & 87, 89, 90, 33, 35, 37, 48, 55, 60, 63, 70, 71, & 72, 73, 74, 75, 76, 77, 81, 83, 84, 85, 86, 87, & 89, 90, 22, 30, 55, 72, 73, 83, 84, 85, 86, 87, & 90, 55, 60, 62, 72, 73, 76, 83, 84, 85, 86, 87, & 90, 23, 35, 68, 70, 72, 73, 74, 76, 77, 80, 81, & 83, 84, 85, 86, 87, 89, 90, 41, 42, 51, 56, 58, & 59, 60, 61, 62, 64, 65, 66, 70, 72, 73, 75, 76, & 77, 78, 79, 80, 82, 83, 84, 85, 86, 87, 89, 90, & 57, 59, 60, 62, 64, 66, 72, 73, 76, 77, 79, 80, & 83, 84, 85, 86, 87, 89, 90, 45, 46, 57, 67, 77, & 79, 83, 84, 85, 86, 87, 89, 90, 26, 44, 58, 77, & 78, 79, 80, 82, 83, 84, 85, 86, 87, 89, 90, 45, & 46, 57, 79, 83, 84, 85, 86, 87, 89, 90, 45, 46, & 57, 67, 79, 80, 83, 84, 85, 86, 87, 89, 90, 18, & 33, 68, 70, 72, 73, 76, 77, 79, 80, 81, 82, 83, & 84, 85, 86, 87, 89, 90, 18, 19, 21, 24, 25, 29, & 38, 39, 43, 44, 51, 52, 53, 54, 55, 56, 57, 59, & 60, 62, 63, 64, 65, 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, 21, 26, 27, 42, 49, 51 /) INTEGER, PARAMETER, DIMENSION(321) :: LU_ICOL_2 = (/ & 52, 53, 55, 56, 57, 59, 60, 61, 62, 63, 64, 65, & 66, 67, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, & 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 17, & 20, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, & 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, & 51, 52, 53, 55, 56, 57, 58, 59, 60, 61, 62, 63, & 64, 65, 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, 24, 25, 51, 52, 53, 55, 56, 57, 59, & 60, 62, 63, 64, 65, 66, 67, 70, 72, 73, 74, 76, & 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, & 89, 90, 21, 27, 28, 34, 48, 49, 51, 52, 53, 55, & 56, 57, 59, 60, 62, 63, 64, 65, 66, 67, 70, 71, & 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, & 84, 85, 86, 87, 88, 89, 90, 21, 22, 23, 25, 29, & 34, 46, 48, 54, 57, 59, 64, 65, 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, 19, 36, 50, & 51, 58, 64, 65, 66, 77, 79, 80, 82, 83, 84, 85, & 86, 87, 88, 89, 90, 25, 44, 46, 48, 58, 77, 78, & 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, & 16, 17, 20, 22, 23, 24, 26, 27, 28, 29, 30, 31, & 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, & 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 56, & 57, 58, 59, 61, 62, 63, 64, 65, 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 /) INTEGER, PARAMETER, DIMENSION(1041) :: LU_ICOL = (/& LU_ICOL_0, LU_ICOL_1, LU_ICOL_2 /) INTEGER, PARAMETER, DIMENSION(91) :: LU_CROW = (/ & 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 24, & 27, 35, 46, 49, 51, 54, 57, 60, 64, 67, 70, 73, & 77, 80, 84, 88, 92, 96,100,104,108,112,115,119, & 123,127,131,135,139,143,147,153,158,163,167,190, & 196,206,220,225,235,245,270,277,284,290,300,308, & 316,336,346,357,365,376,386,396,423,473,484,507, & 518,530,548,577,596,609,624,635,648,667,715,756, & 820,855,896,934,954,973,1042 /) INTEGER, PARAMETER, DIMENSION(91) :: LU_DIAG = (/ & 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 24, & 27, 35, 46, 49, 51, 54, 57, 60, 64, 67, 70, 73, & 77, 80, 84, 88, 92, 96,100,104,108,112,115,119, & 123,127,131,135,139,143,147,153,158,163,170,191, & 197,208,220,227,237,250,270,277,284,292,301,308, & 321,338,348,358,368,378,388,410,452,475,492,510, & 522,536,563,585,600,613,627,640,658,706,748,813, & 849,891,930,951,971,1041,1042 /) END MODULE gckpp_JacobianSP