! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! 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 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_JacobianSP PUBLIC SAVE ! Sparse Jacobian Data INTEGER, PARAMETER, DIMENSION(360) :: LU_IROW_0 = (/ & 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, 18, 18, 18, 19, 19, 19, 20, 20, 20, & 21, 21, 21, 21, 21, 21, 21, 21, 22, 22, 22, 22, & 22, 22, 22, 22, 22, 22, 22, 23, 23, 24, 24, 25, & 25, 26, 26, 27, 27, 27, 28, 28, 28, 29, 29, 29, & 30, 30, 30, 30, 31, 31, 31, 31, 31, 32, 32, 32, & 32, 32, 33, 33, 33, 33, 33, 34, 34, 34, 35, 35, & 35, 36, 36, 36, 37, 37, 37, 37, 38, 38, 38, 39, & 39, 39, 39, 40, 40, 40, 40, 41, 41, 41, 41, 42, & 42, 42, 42, 43, 43, 43, 43, 44, 44, 44, 45, 45, & 45, 45, 46, 46, 46, 46, 47, 47, 47, 47, 48, 48, & 48, 48, 49, 49, 49, 49, 50, 50, 50, 50, 51, 51, & 51, 51, 52, 52, 52, 52, 53, 53, 53, 53, 54, 54, & 54, 54, 55, 55, 55, 55, 56, 56, 56, 56, 56, 56, & 57, 57, 57, 57, 57, 58, 58, 58, 58, 58, 59, 59, & 59, 59, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, & 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, 60, & 60, 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, 63, 63, 63, 64, 64, 64, 64, 64, & 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, 67, 67, 67, 67, 67, 67, & 67, 67, 67, 67, 67, 67, 67, 67, 67, 68, 68, 68, & 68, 68, 68, 68, 69, 69, 69, 69, 69, 69, 69, 70, & 70, 70, 70, 70, 70, 71, 71, 71, 71, 71, 71, 71, & 71, 71, 71, 72, 72, 72, 72, 72, 72, 72, 72, 73, & 73, 73, 73, 73, 73, 73, 73, 74, 74, 74, 74, 74 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_IROW_1 = (/ & 74, 74, 74, 74, 74, 75, 75, 75, 75, 75, 75, 75, & 75, 76, 76, 76, 76, 76, 76, 76, 76, 76, 76, 76, & 77, 77, 77, 77, 77, 77, 77, 77, 77, 77, 77, 77, & 77, 77, 77, 77, 77, 77, 77, 77, 77, 77, 78, 78, & 78, 78, 78, 78, 78, 78, 78, 78, 79, 79, 79, 79, & 79, 79, 79, 79, 79, 79, 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, & 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, 82, 82, 82, 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, 85, 85, 85, & 85, 85, 85, 85, 85, 85, 85, 85, 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, 88, 88, 88, 88, 88, 88, 88, 88, 88, 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, 91, 91, 91, & 91, 91, 91, 91, 91, 91, 91, 91, 91, 91, 92, 92, & 92, 92, 92, 92, 92, 92, 92, 92, 92, 93, 93, 93, & 93, 93, 93, 93, 93, 93, 93, 93, 93, 93, 94, 94, & 94, 94, 94, 94, 94, 94, 94, 94, 94, 94, 94, 94, & 94, 94, 94, 94, 94, 94, 95, 95, 95, 95, 95, 95 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_IROW_2 = (/ & 95, 95, 95, 95, 95, 95, 95, 95, 95, 95, 95, 95, & 95, 95, 95, 95, 95, 95, 95, 95, 95, 95, 95, 95, & 95, 95, 95, 95, 95, 95, 95, 95, 95, 95, 95, 96, & 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, & 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, & 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, & 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, & 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, 96, & 96, 96, 96, 96, 96, 96, 97, 97, 97, 97, 97, 97, & 97, 97, 97, 97, 97, 97, 97, 97, 97, 98, 98, 98, & 98, 98, 98, 98, 98, 98, 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, 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,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,101,101,101,101,101,101,101,101,101,101, & 101,101,101,101,101,101,101,101,101,101,101,101, & 101,101,101,101,101,101,101,101,101,101,101,102, & 102,102,102,102,102,102,102,102,102,102,102,102, & 102,102,102,102,102,102,102,102,102,102,102,102, & 102,102,102,102,102,102,102,102,102,102,102,102, & 102,102,102,102,102,102,102,102,102,102,102,102, & 102,102,102,102,102,102,102,102,102,102,102,102, & 102,102,102,102,102,102,102,102,102,102,103,103, & 103,103,103,103,103,103,103,103,103,103,103,103, & 103,103,103,103,103,103,103,103,103,103,103,103 /) INTEGER, PARAMETER, DIMENSION(12) :: LU_IROW_3 = (/ & 103,103,103,103,103,103,103,103,103,103,103,103 /) INTEGER, PARAMETER, DIMENSION(1092) :: LU_IROW = (/& LU_IROW_0, LU_IROW_1, LU_IROW_2, LU_IROW_3 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_ICOL_0 = (/ & 1, 2, 82, 3, 29, 4, 67, 5, 38, 6, 99, 7, & 98, 8, 34, 9, 58, 10, 27, 11, 86, 12, 30,102, & 13, 44,102, 14, 59,102, 15, 31, 96, 16, 31,103, & 17, 32, 96, 18, 32,103, 19, 33, 96, 20, 33,103, & 21, 59, 60, 91, 98,101,102,103, 22, 27, 29, 34, & 38, 58, 67, 82, 86, 98, 99, 23,102, 24,102, 25, & 102, 26,102, 27, 90, 99, 28, 94, 99, 29, 96,102, & 30, 44,100,102, 25, 31, 96,102,103, 23, 32, 96, & 102,103, 24, 33, 96,102,103, 34, 99,101, 35,100, & 102, 36,100,102, 37, 99,102,103, 38, 99,100, 39, & 95, 96,102, 40, 85, 96,102, 41, 65, 96,102, 42, & 66, 96,102, 43, 90, 96,102, 44,100,102, 45, 96, & 97,102, 46, 87, 96,102, 47, 94, 96,102, 48, 96, & 101,102, 49, 96, 99,102, 50, 83, 96,102, 51, 76, & 96,102, 52, 78, 96,102, 53, 72, 96,102, 54, 79, & 96,102, 55, 75, 96,102, 56, 64, 69, 80, 96,102, & 57, 70, 92, 96,102, 58, 97, 98, 99,102, 59, 98, & 100,102, 53, 54, 59, 60, 61, 62, 63, 71, 72, 79, & 81, 82, 84, 88, 91, 93, 95, 96, 98,100,101,102, & 103, 59, 61, 93, 98,100,102, 42, 62, 66, 85, 86, & 95, 96,101,102,103, 54, 55, 63, 64, 71, 75, 79, & 80, 95, 96, 98,101,102,103, 64, 95, 96,101,103, & 26, 41, 65, 85, 86, 95, 96,101,102,103, 26, 42, & 66, 85, 86, 95, 96,101,102,103, 35, 36, 38, 44, & 64, 67, 70, 72, 75, 78, 79, 80, 81, 82, 84, 88, & 89, 91, 95, 96, 99,100,101,102,103, 68, 86, 95, & 96,101,102,103, 69, 91, 95, 96,100,101,103, 70, & 92, 95, 96,101,103, 57, 70, 71, 92, 95, 96, 98, & 101,102,103, 53, 72, 91, 95, 96,101,102,103, 73, & 89, 95, 96,100,101,102,103, 62, 66, 74, 85, 86 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_ICOL_1 = (/ & 95, 96,101,102,103, 55, 75, 93, 95, 96,101,102, & 103, 51, 61, 76, 93, 95, 96, 98,100,101,102,103, & 53, 54, 58, 71, 72, 74, 77, 79, 80, 83, 85, 86, & 91, 92, 95, 96, 97, 98, 99,101,102,103, 52, 59, & 78, 95, 96, 98,100,101,102,103, 54, 71, 79, 92, & 95, 96, 98,101,102,103, 56, 64, 69, 80, 91, 95, & 96,100,101,102,103, 40, 41, 51, 52, 55, 56, 64, & 65, 68, 69, 75, 76, 78, 80, 81, 83, 85, 86, 91, & 93, 95, 96, 98,100,101,102,103, 1, 39, 44, 45, & 47, 48, 50, 53, 55, 57, 58, 59, 61, 63, 64, 65, & 66, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, & 79, 80, 82, 83, 84, 85, 86, 87, 89, 90, 91, 92, & 93, 94, 95, 96, 97, 98, 99,100,101,102,103, 50, & 61, 83, 93, 95, 96, 98,100,101,102,103, 43, 46, & 50, 61, 68, 73, 76, 83, 84, 85, 86, 87, 88, 89, & 90, 93, 95, 96, 98,100,101,102,103, 35, 40, 68, & 85, 86, 95, 96,100,101,102,103, 68, 73, 74, 85, & 86, 89, 95, 96,100,101,102,103, 36, 46, 81, 83, & 85, 86, 87, 89, 90, 91, 93, 95, 96, 98,100,101, & 102,103, 54, 55, 64, 69, 71, 72, 73, 74, 75, 77, & 79, 80, 83, 85, 86, 88, 89, 91, 92, 93, 95, 96, & 97, 98, 99,100,101,102,103, 70, 72, 73, 74, 75, & 79, 85, 86, 89, 91, 92, 93, 95, 96, 98,100,101, & 102,103, 27, 43, 81, 83, 85, 86, 89, 90, 91, 92, & 93, 95, 96, 98, 99,100,101,102,103, 57, 59, 70, & 78, 91, 92, 95, 96, 98,100,101,102,103, 57, 59, & 70, 92, 95, 96, 98,100,101,102,103, 57, 59, 70, & 78, 92, 93, 95, 96, 98,100,101,102,103, 28, 47, & 63, 64, 71, 75, 79, 80, 91, 92, 93, 94, 95, 96, & 98, 99,100,101,102,103, 34, 39, 44, 48, 61, 62 /) INTEGER, PARAMETER, DIMENSION(360) :: LU_ICOL_2 = (/ & 64, 65, 66, 68, 69, 70, 72, 73, 74, 75, 76, 78, & 79, 80, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, & 93, 94, 95, 96, 97, 98, 99,100,101,102,103, 29, & 30, 31, 32, 33, 39, 40, 41, 42, 43, 44, 46, 47, & 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, & 61, 62, 63, 64, 65, 66, 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, 45, 58, 71, 91, 92, 93, & 95, 96, 97, 98, 99,100,101,102,103, 38, 58, 59, & 61, 71, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, & 100,101,102,103, 27, 28, 34, 37, 38, 49, 51, 52, & 56, 58, 64, 65, 66, 67, 68, 69, 70, 72, 73, 74, & 75, 76, 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, 34, 35, 36, 38, 44, 49, 59, 61, & 67, 70, 72, 75, 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, 34, 45, 48, 55, 62, 64, & 65, 66, 68, 69, 70, 72, 73, 74, 75, 76, 77, 78, & 79, 80, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, & 93, 94, 95, 96, 97, 98, 99,100,101,102,103, 23, & 24, 25, 26, 29, 30, 35, 36, 37, 39, 40, 41, 42, & 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, & 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, & 67, 69, 70, 71, 72, 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, 31, 32, & 33, 37, 38, 64, 65, 66, 68, 69, 70, 72, 73, 74, & 75, 76, 78, 79, 80, 83, 85, 86, 87, 89, 90, 91 /) INTEGER, PARAMETER, DIMENSION(12) :: LU_ICOL_3 = (/ & 92, 93, 94, 95, 96, 97, 98, 99,100,101,102,103 /) INTEGER, PARAMETER, DIMENSION(1092) :: LU_ICOL = (/& LU_ICOL_0, LU_ICOL_1, LU_ICOL_2, LU_ICOL_3 /) INTEGER, PARAMETER, DIMENSION(104) :: LU_CROW = (/ & 1, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, & 25, 28, 31, 34, 37, 40, 43, 46, 49, 57, 68, 70, & 72, 74, 76, 79, 82, 85, 89, 94, 99,104,107,110, & 113,117,120,124,128,132,136,140,143,147,151,155, & 159,163,167,171,175,179,183,187,193,198,203,207, & 230,236,246,260,265,275,285,310,317,324,330,340, & 348,356,366,374,385,407,417,427,438,465,516,527, & 550,561,573,591,620,639,658,671,682,695,715,756, & 823,838,857,905,943,984,1055,1093 /) INTEGER, PARAMETER, DIMENSION(104) :: LU_DIAG = (/ & 1, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, & 25, 28, 31, 34, 37, 40, 43, 46, 49, 57, 68, 70, & 72, 74, 76, 79, 82, 85, 90, 95,100,104,107,110, & 113,117,120,124,128,132,136,140,143,147,151,155, & 159,163,167,171,175,179,183,187,193,198,203,210, & 231,237,248,260,267,277,290,310,317,324,332,341, & 348,358,367,376,391,409,419,430,452,495,518,535, & 553,565,579,606,628,646,662,674,687,706,747,815, & 831,851,900,939,981,1053,1092,1093 /) END MODULE gckpp_JacobianSP