! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! Sparse Hessian 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_HessianSP.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_HessianSP PUBLIC SAVE ! Hessian Sparse Data ! INTEGER, PARAMETER, DIMENSION(360) :: IHESS_I_0 = (/ & 12, 13, 14, 15, 16, 17, 17, 17, 18, 19, 20, 21, & 21, 21, 22, 23, 24, 25, 25, 25, 26, 27, 28, 29, & 29, 29, 30, 31, 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, 42, 42, 42, 44, 45, & 46, 46, 46, 47, 47, 48, 48, 48, 49, 50, 50, 51, & 51, 52, 52, 53, 53, 54, 54, 55, 56, 56, 57, 57, & 58, 58, 59, 59, 60, 60, 61, 61, 62, 62, 63, 63, & 64, 64, 65, 65, 66, 66, 66, 67, 67, 68, 68, 68, & 68, 69, 69, 70, 70, 71, 71, 72, 72, 72, 72, 72, & 72, 72, 72, 72, 72, 72, 73, 73, 74, 74, 75, 75, & 75, 76, 76, 77, 77, 77, 77, 78, 78, 78, 78, 79, & 79, 79, 80, 80, 80, 81, 81, 81, 82, 82, 82, 83, & 83, 83, 83, 84, 84, 84, 85, 85, 85, 85, 85, 85, & 85, 85, 86, 86, 86, 86, 87, 87, 87, 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, 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, 92, 92, 92, 92, 93, 93, 93, 93, 93, 94, 94, & 94, 94, 94, 94, 94, 94, 94, 95, 95, 95, 96, 96, & 96, 96, 96, 96, 96, 96, 96, 97, 97, 97, 97, 97, & 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,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,102,102,102,102,102,102, & 103,103,103,103,104,104,104,104,104,104,104,105 /) INTEGER, PARAMETER, DIMENSION(360) :: IHESS_I_1 = (/ & 105,105,105,105,105,106,106,106,106,106,106,107, & 107,107,107,108,108,108,108,108,108,108,109,109, & 109,109,109,109,110,110,110,110,110,110,111,111, & 111,111,111,111,112,112,112,112,112,112,112,112, & 112,112,112,112,112,112,112,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, & 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,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,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,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,119,120,120,120,120,120,120,120,120,120,120, & 121,121,121,121,121,121,121,121,121,121,121,122, & 122,122,122,122,122,122,122,123,123,123,123,123, & 123,123,124,124,124,124,124,124,124,124,124,124, & 125,125,125,125,125,125,125,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,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 /) INTEGER, PARAMETER, DIMENSION(360) :: IHESS_I_2 = (/ & 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,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,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,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,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,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 /) INTEGER, PARAMETER, DIMENSION(77) :: IHESS_I_3 = (/ & 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,136,136,136, & 136,136,136,136,136,136,136,136,136,136,137,137, & 137,137,137,137,137 /) INTEGER, PARAMETER, DIMENSION(1157) :: IHESS_I = (/& IHESS_I_0, IHESS_I_1, IHESS_I_2, IHESS_I_3 /) INTEGER, PARAMETER, DIMENSION(360) :: IHESS_J_0 = (/ & 48, 60, 81, 17, 17, 17, 17, 18, 18, 21, 21, 21, & 21, 22, 22, 25, 25, 25, 25, 26, 26, 29, 29, 29, & 29, 30, 30, 93, 68, 70, 72, 75, 75, 78, 81, 84, & 84, 87, 87, 89, 91, 92, 92, 95, 95, 97,100,102, & 102,102,102,102,119,121,122,122,123,130, 44,113, & 46,132,136, 47, 74, 48, 60, 60,128, 92, 95, 51, & 51, 52, 52, 53,132, 54,130,128, 56, 94, 57, 98, & 58,113, 59,128, 60, 60, 61,112, 62,133, 63,108, & 64,122, 75, 95, 66, 87,126, 67,110, 68,103,103, & 125, 69,109, 70,105, 71,117, 72, 84, 84, 92, 92, & 93,100,103,103,119,124, 73,101, 74,123, 75, 75, & 75, 76,111, 76, 77,111,111, 78, 78,111,111, 79, & 79,122, 80,103,107, 81, 81, 81, 82, 82,123, 83, & 83,123,123, 47, 84, 84, 57, 85, 98, 98, 98,126, & 126,126, 81, 86, 86, 86, 78, 87, 87, 88, 94, 98, & 101,102,105,106,109,110,111,116,117,123,123,126, & 89, 92, 92, 94, 97, 98,101,102,105,106,108,109, & 110,111,112,116,117,119,123,126,130,130, 83, 83, & 84, 84, 90, 90,100,111,111, 77, 81, 84, 84, 86, & 90, 90, 91, 97, 97, 97, 97,100,105,105,105,105, & 111,111,111,114,114,115,118,118,119,120,121,122, & 124, 92, 92, 92,104, 82, 93,101,107,107, 44, 56, & 94, 94, 94, 94,126,126,126, 79, 95, 95, 67, 96, & 97, 97, 97, 97,103,103,125, 97, 97, 97, 97,121, & 44, 57, 98, 98, 98, 98,126,126,126, 51, 52, 60, & 83, 87, 90, 99,102,103,107,111,114,115,117,117, & 118,121,127,128,129, 83, 83, 84, 84,100,101,101, & 101,101,103,103,107,107,111,111,111,111,124,125, & 73,101,101,101,101,120, 64,102,102,102,102,102, & 80,103,103,125, 82,103,103,104,105,107,107, 70 /) INTEGER, PARAMETER, DIMENSION(360) :: IHESS_J_1 = (/ & 105,105,105,105,121,106,106,106,106,129,131, 80, & 82,107,107, 63, 85,108,108,108,108,122, 69, 86, & 109,109,109,109, 67, 86,110,110,110,110, 76,111, & 111,111,111,124, 52, 52, 61,112,112,112,112,112, & 113,113,113,113,126,126,126, 58,113,113,113,113, & 113,127,127, 60, 60, 62, 68, 79, 81, 82, 84, 84, & 86, 92, 92, 93, 94, 95, 97, 97, 97, 97, 98,100, & 101,101,101,101,102,103,103,105,105,105,105,106, & 107,107,108,108,108,108,109,109,109,110,110,110, & 111,112,113,114,114,115,116,116,116,117,117,120, & 121,122,122,122,122,123,123,123,123,123,124,126, & 130,133,133,133, 61, 86,106,106,106,106,109,109, & 109,110,110,110,112,112,112,112,115,115,116,116, & 116,120,126,126,126, 80,116,116,116,116,132, 71, & 81,117,117,117,117,117,117, 63, 83, 83, 84, 84, & 92, 92, 93, 96, 97, 97, 97, 97,101,101,101,101, & 105,105,105,105,108,108,111,111,111,111,118,118, & 119,120,121,124, 69, 70, 77, 83, 83, 84, 84, 92, & 92, 95, 95,102,102,102,102,102,103,103,105,105, & 105,105,107,107,108,109,109,111,111,111,119,122, & 124,125, 81,117,117,120,120,123,123,123,123,123, & 81,117,117,121,121,121,123,123,123,123,123, 64, & 121,121,122,122,122,122,122, 74, 81,123,123,123, & 123,123, 74,117,117,117,123,123,123,123,124,124, & 117,117,117,117,123,125,125, 51, 51, 66,116,116, & 116,126,126,126,126, 56, 66, 87, 87, 88, 94, 94, & 94, 97,101,109,109,110,110,111,111,116,116,116, & 117,117,126,126,126,127,127, 53, 59, 68, 71, 75, & 75, 79, 83, 84, 87, 92, 92, 92, 94, 95, 95, 96, & 97, 97, 97, 97, 98,101,102,103,103,105,106,107 /) INTEGER, PARAMETER, DIMENSION(360) :: IHESS_J_2 = (/ & 107,108,109,110,110,110,111,112,113,113,116,116, & 116,117,117,117,117,117,122,122,123,126,128,128, & 128,128,128,129,129,129,129,130,133,134,134, 51, & 52, 60, 81, 86, 90, 92, 92, 93, 95, 99,103,103, & 114,115,117,118,121,127,128,128,129,129,129,129, & 129, 54, 75, 75, 94, 97, 98,101,101,101,101,102, & 105,106,106,106,106,108,108,108,108,109,110,111, & 112,113,115,115,116,117,118,118,122,122,123,126, & 128,130,130,130,130, 71, 73, 77, 78,101,101,101, & 102,103,103,105,106,106,111,111,123,126,126,126, & 129,131, 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, 78, 79, 80, 81, & 81, 82, 82, 83, 84, 84, 85, 86, 86, 87, 88, 89, & 90, 91, 92, 93, 96, 97, 97, 99,100,101,102,102, & 102,102,102,103,104,105,106,107,108,111,113,114, & 115,117,118,119,120,120,121,121,122,123,124,124, & 125,125,127,128,129,129,130,131,132,132,132,132, & 132,134,135, 60, 60, 62, 86, 89, 94, 94, 97, 97, & 98, 98,101,101,102,102,105,105,106,106,108,108, & 109,109,110,110,111,111,112,112,113,113,116,116, & 117,117,119,122,122,122,122,123,123,126,126,126, & 130,130,130,130,133,133,133,133, 75, 83, 84, 87, & 92, 94, 95, 97, 98,101,102,103,105,106,107,108, & 109,110,111,112,113,116,117,122,123,126,128,129, & 130,132,133,134,134, 78, 79, 81, 82, 86,113,120, & 121,122,124,125,128,130,132,132,133,134,135, 46, & 48, 72, 75, 77, 78, 79, 81, 83, 84, 84, 86, 87, & 87, 88, 90, 90, 91, 92, 94, 94, 94, 94, 95, 95, & 97, 97, 97, 98, 98, 98, 98,100,101,101,101,101 /) INTEGER, PARAMETER, DIMENSION(77) :: IHESS_J_3 = (/ & 102,102,103,103,105,105,105,105,106,106,107,107, & 108,108,109,109,109,109,110,110,111,111,111,111, & 112,112,112,112,112,113,113,114,114,115,116,116, & 117,117,117,117,117,117,119,120,121,122,122,123, & 123,123,123,123,126,126,126,126,128,129,129,130, & 130,132,132,133,133,133,133,134,135,136, 71, 97, & 106,108,116,126,132 /) INTEGER, PARAMETER, DIMENSION(1157) :: IHESS_J = (/& IHESS_J_0, IHESS_J_1, IHESS_J_2, IHESS_J_3 /) INTEGER, PARAMETER, DIMENSION(360) :: IHESS_K_0 = (/ & 132,132,132,136,134,134,136,132,132,136,134,134, & 136,132,132,136,134,134,136,132,132,136,134,134, & 136,132,132,132,132,132,132,134,136,135,135,134, & 136,134,136,132,132,134,136,134,136,130,132,102, & 130,133,134,136,132,135,134,136,130,134,132,128, & 132,132,136,132,132,132,129,132,130,128,134,129, & 132,129,132,132,134,132,136,129,132,136,132,136, & 132,136,132,136,129,132,132,136,132,136,132,136, & 132,136,128,134,132,136,136,132,136,132,134,136, & 135,132,136,132,136,132,136,132,134,136,134,136, & 132,132,134,136,132,135,132,136,132,136,128,134, & 136,132,136,132,132,130,134,132,135,130,134,132, & 135,128,132,136,136,129,132,135,132,135,134,134, & 136,130,134,132,134,136,132,132,130,133,134,130, & 133,134,135,129,132,135,132,134,136,132,133,133, & 133,133,133,133,133,133,133,133,133,123,133,133, & 132,134,136,130,130,130,130,130,130,130,130,130, & 130,130,130,130,130,132,130,130,133,136,134,136, & 134,136,129,132,132,130,134,132,135,134,136,135, & 129,132,132,130,133,134,136,132,130,133,134,136, & 130,133,134,129,132,132,129,132,132,135,135,134, & 135,128,134,136,132,135,132,134,134,136,132,132, & 130,133,134,136,130,133,134,132,134,136,132,132, & 130,133,134,136,134,136,135,130,133,134,136,129, & 132,132,130,133,134,136,130,133,134,129,129,129, & 134,134,129,132,134,134,134,134,129,129,130,133, & 129,129,129,132,131,134,136,134,136,132,130,133, & 134,136,134,136,134,136,130,133,134,136,135,135, & 132,130,133,134,136,132,132,102,130,133,134,136, & 132,134,136,132,135,134,136,132,134,134,136,132 /) INTEGER, PARAMETER, DIMENSION(360) :: IHESS_K_1 = (/ & 130,133,134,136,132,130,133,134,136,131,132,132, & 132,134,136,132,132,130,133,134,136,136,132,132, & 130,133,134,136,132,129,130,133,134,136,132,130, & 133,134,136,132,129,132,132,112,130,133,134,136, & 130,133,134,136,130,133,134,132,128,130,133,134, & 136,129,132,129,132,132,132,135,135,135,134,136, & 135,134,136,132,133,134,130,133,134,136,133,132, & 130,133,134,136,133,134,136,130,133,134,136,133, & 134,136,130,133,134,136,130,133,134,130,133,134, & 133,133,133,129,132,132,130,133,134,130,133,135, & 135,130,133,134,136,123,130,133,134,136,135,133, & 133,133,134,135,132,135,130,133,134,136,130,133, & 134,130,133,134,112,130,133,134,129,132,130,133, & 134,135,130,133,134,132,130,133,134,136,137,132, & 129,117,129,130,133,134,136,132,134,136,134,136, & 134,136,132,132,130,133,134,136,130,133,134,136, & 130,133,134,136,130,133,130,133,134,136,129,132, & 132,135,135,135,132,132,132,134,136,134,136,134, & 136,134,136,102,130,133,134,136,134,136,130,133, & 134,136,134,136,133,130,133,130,133,134,132,136, & 135,135,135,130,133,132,135,123,130,133,134,136, & 135,130,133,129,132,135,123,130,133,134,136,132, & 129,132,128,130,133,134,136,132,132,123,130,133, & 134,136,132,117,129,134,123,130,133,134,132,135, & 117,129,134,136,134,132,135,129,132,132,130,133, & 134,130,133,134,136,132,132,134,136,132,130,133, & 134,130,136,130,133,130,133,130,134,130,133,134, & 130,133,130,133,134,129,132,132,132,132,132,128, & 134,135,134,134,134,128,134,136,134,134,136,132, & 130,133,134,136,134,134,134,134,136,134,134,134 /) INTEGER, PARAMETER, DIMENSION(360) :: IHESS_K_2 = (/ & 136,134,134,130,133,134,134,134,128,134,130,133, & 134,117,129,130,133,134,128,134,134,134,129,130, & 132,135,136,129,132,134,136,134,134,135,136,129, & 129,129,129,129,129,134,136,132,134,132,134,136, & 129,129,129,129,129,129,129,135,129,131,132,134, & 136,132,134,136,130,130,130,130,133,134,136,130, & 130,130,133,134,136,130,133,134,136,130,130,130, & 130,130,129,132,130,130,129,132,133,134,130,130, & 130,130,133,134,136,132,132,132,135,130,133,136, & 130,134,136,130,130,133,130,133,130,130,133,134, & 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,136,132,132,132,135,132,132,132, & 135,132,135,136,134,136,132,132,135,136,132,132, & 132,132,136,132,132,134,136,132,132,136,102,130, & 133,134,136,136,132,136,136,136,136,136,136,132, & 132,136,132,132,132,135,132,135,136,136,132,135, & 132,135,132,132,132,136,136,132,132,134,135,136, & 137,136,136,129,132,132,135,132,130,133,130,133, & 130,133,130,133,130,133,130,133,130,133,130,133, & 130,133,130,133,130,133,130,133,130,133,130,133, & 130,133,132,130,133,134,136,130,133,130,133,134, & 130,133,134,136,133,134,135,136,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,129,134, & 134,134,134,135,136,135,135,135,135,135,136,135, & 135,136,135,135,135,136,132,135,135,135,136,132, & 132,132,136,132,135,135,135,134,134,136,135,134, & 136,132,129,132,132,136,130,133,134,136,134,136, & 130,133,136,130,133,134,136,132,130,133,134,136 /) INTEGER, PARAMETER, DIMENSION(77) :: IHESS_K_3 = (/ & 133,136,134,136,130,133,134,136,133,136,134,136, & 133,136,130,133,134,136,133,136,130,133,134,136, & 112,130,133,134,136,133,136,129,132,132,133,136, & 117,129,130,133,134,136,132,135,135,133,136,123, & 130,133,134,136,130,133,134,136,136,132,136,133, & 136,135,136,133,134,135,136,136,136,136,132,133, & 134,134,136,134,137 /) INTEGER, PARAMETER, DIMENSION(1157) :: IHESS_K = (/& IHESS_K_0, IHESS_K_1, IHESS_K_2, IHESS_K_3 /) END MODULE gckpp_HessianSP