! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! The ODE Jacobian of Chemical Model 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_Jacobian.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_Jacobian USE gckpp_Parameters USE gckpp_JacobianSP IMPLICIT NONE CONTAINS ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! Jac_SP - the Jacobian of Variables in sparse matrix representation ! Arguments : ! V - Concentrations of variable species (local) ! F - Concentrations of fixed species (local) ! RCT - Rate constants (local) ! JVS - sparse Jacobian of variables ! ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ SUBROUTINE Jac_SP ( V, F, RCT, JVS ) ! V - Concentrations of variable species (local) REAL(kind=dp) :: V(NVAR) ! F - Concentrations of fixed species (local) REAL(kind=dp) :: F(NFIX) ! RCT - Rate constants (local) REAL(kind=dp) :: RCT(NREACT) ! JVS - sparse Jacobian of variables REAL(kind=dp) :: JVS(LU_NONZERO) ! Local variables ! B - Temporary array REAL(kind=dp) :: B(548) ! B(1) = dA(1)/dV(98) B(1) = RCT(1)*V(103) ! B(2) = dA(1)/dV(103) B(2) = RCT(1)*V(98) ! B(3) = dA(2)/dV(98) B(3) = RCT(2)*V(102) ! B(4) = dA(2)/dV(102) B(4) = RCT(2)*V(98) ! B(5) = dA(3)/dV(96) B(5) = RCT(3)*V(98) ! B(6) = dA(3)/dV(98) B(6) = RCT(3)*V(96) ! B(7) = dA(4)/dV(98) B(7) = RCT(4)*V(99) ! B(8) = dA(4)/dV(99) B(8) = RCT(4)*V(98) ! B(9) = dA(5)/dV(95) B(9) = RCT(5)*V(98) ! B(10) = dA(5)/dV(98) B(10) = RCT(5)*V(95) ! B(11) = dA(6)/dV(102) B(11) = RCT(6)*2*V(102) ! B(12) = dA(7)/dV(102) B(12) = RCT(7)*2*V(102) ! B(13) = dA(8)/dV(96) B(13) = RCT(8)*V(102) ! B(14) = dA(8)/dV(102) B(14) = RCT(8)*V(96) ! B(15) = dA(9)/dV(29) B(15) = RCT(9)*V(102) ! B(16) = dA(9)/dV(102) B(16) = RCT(9)*V(29) ! B(17) = dA(10)/dV(96) B(17) = RCT(10)*V(103) ! B(18) = dA(10)/dV(103) B(18) = RCT(10)*V(96) ! B(19) = dA(11)/dV(96) B(19) = RCT(11)*2*V(96) ! B(20) = dA(12)/dV(102) B(20) = RCT(12)*F(9) ! B(22) = dA(13)/dV(60) B(22) = RCT(13)*V(102) ! B(23) = dA(13)/dV(102) B(23) = RCT(13)*V(60) ! B(24) = dA(14)/dV(102) B(24) = RCT(14)*F(2) ! B(26) = dA(15)/dV(95) B(26) = RCT(15)*V(103) ! B(27) = dA(15)/dV(103) B(27) = RCT(15)*V(95) ! B(28) = dA(16)/dV(95) B(28) = RCT(16)*V(96) ! B(29) = dA(16)/dV(96) B(29) = RCT(16)*V(95) ! B(30) = dA(17)/dV(95) B(30) = RCT(17)*2*V(95) ! B(31) = dA(18)/dV(95) B(31) = RCT(18)*2*V(95) ! B(32) = dA(19)/dV(39) B(32) = RCT(19)*V(102) ! B(33) = dA(19)/dV(102) B(33) = RCT(19)*V(39) ! B(34) = dA(20)/dV(39) B(34) = RCT(20)*V(102) ! B(35) = dA(20)/dV(102) B(35) = RCT(20)*V(39) ! B(36) = dA(21)/dV(82) B(36) = RCT(21)*V(102) ! B(37) = dA(21)/dV(102) B(37) = RCT(21)*V(82) ! B(38) = dA(22)/dV(99) B(38) = RCT(22)*V(102) ! B(39) = dA(22)/dV(102) B(39) = RCT(22)*V(99) ! B(40) = dA(23)/dV(67) B(40) = RCT(23)*V(102) ! B(41) = dA(23)/dV(102) B(41) = RCT(23)*V(67) ! B(42) = dA(24)/dV(102) B(42) = RCT(24)*V(103) ! B(43) = dA(24)/dV(103) B(43) = RCT(24)*V(102) ! B(44) = dA(25)/dV(37) B(44) = RCT(25)*V(102) ! B(45) = dA(25)/dV(102) B(45) = RCT(25)*V(37) ! B(46) = dA(26)/dV(96) B(46) = RCT(26)*V(99) ! B(47) = dA(26)/dV(99) B(47) = RCT(26)*V(96) ! B(48) = dA(27)/dV(49) B(48) = RCT(27) ! B(49) = dA(28)/dV(49) B(49) = RCT(28)*V(102) ! B(50) = dA(28)/dV(102) B(50) = RCT(28)*V(49) ! B(51) = dA(29)/dV(96) B(51) = RCT(29)*V(100) ! B(52) = dA(29)/dV(100) B(52) = RCT(29)*V(96) ! B(53) = dA(30)/dV(100) B(53) = RCT(30)*V(103) ! B(54) = dA(30)/dV(103) B(54) = RCT(30)*V(100) ! B(55) = dA(31)/dV(100) B(55) = RCT(31)*V(102) ! B(56) = dA(31)/dV(102) B(56) = RCT(31)*V(100) ! B(57) = dA(32)/dV(99) B(57) = RCT(32)*V(100) ! B(58) = dA(32)/dV(100) B(58) = RCT(32)*V(99) ! B(59) = dA(33)/dV(38) B(59) = RCT(33) ! B(60) = dA(34)/dV(102) B(60) = RCT(34)*F(11) ! B(62) = dA(35)/dV(102) B(62) = RCT(35)*F(13) ! B(64) = dA(36)/dV(99) B(64) = RCT(36)*V(100) ! B(65) = dA(36)/dV(100) B(65) = RCT(36)*V(99) ! B(66) = dA(37)/dV(82) B(66) = RCT(37)*V(100) ! B(67) = dA(37)/dV(100) B(67) = RCT(37)*V(82) ! B(68) = dA(38)/dV(84) B(68) = RCT(38)*V(102) ! B(69) = dA(38)/dV(102) B(69) = RCT(38)*V(84) ! B(70) = dA(39)/dV(84) B(70) = RCT(39)*V(100) ! B(71) = dA(39)/dV(100) B(71) = RCT(39)*V(84) ! B(72) = dA(40)/dV(99) B(72) = RCT(40)*V(101) ! B(73) = dA(40)/dV(101) B(73) = RCT(40)*V(99) ! B(74) = dA(41)/dV(34) B(74) = RCT(41) ! B(75) = dA(42)/dV(101) B(75) = RCT(42)*V(103) ! B(76) = dA(42)/dV(103) B(76) = RCT(42)*V(101) ! B(77) = dA(43)/dV(36) B(77) = RCT(43)*V(102) ! B(78) = dA(43)/dV(102) B(78) = RCT(43)*V(36) ! B(79) = dA(44)/dV(87) B(79) = RCT(44)*V(103) ! B(80) = dA(44)/dV(103) B(80) = RCT(44)*V(87) ! B(81) = dA(45)/dV(26) B(81) = RCT(45)*V(102) ! B(82) = dA(45)/dV(102) B(82) = RCT(45)*V(26) ! B(83) = dA(46)/dV(26) B(83) = RCT(46)*V(102) ! B(84) = dA(46)/dV(102) B(84) = RCT(46)*V(26) ! B(85) = dA(47)/dV(65) B(85) = RCT(47)*V(103) ! B(86) = dA(47)/dV(103) B(86) = RCT(47)*V(65) ! B(87) = dA(48)/dV(83) B(87) = RCT(48)*V(103) ! B(88) = dA(48)/dV(103) B(88) = RCT(48)*V(83) ! B(89) = dA(49)/dV(35) B(89) = RCT(49)*V(102) ! B(90) = dA(49)/dV(102) B(90) = RCT(49)*V(35) ! B(91) = dA(50)/dV(85) B(91) = RCT(50)*V(103) ! B(92) = dA(50)/dV(103) B(92) = RCT(50)*V(85) ! B(93) = dA(51)/dV(85) B(93) = RCT(51)*V(103) ! B(94) = dA(51)/dV(103) B(94) = RCT(51)*V(85) ! B(95) = dA(52)/dV(68) B(95) = RCT(52)*V(103) ! B(96) = dA(52)/dV(103) B(96) = RCT(52)*V(68) ! B(97) = dA(53)/dV(74) B(97) = RCT(53)*V(103) ! B(98) = dA(53)/dV(103) B(98) = RCT(53)*V(74) ! B(99) = dA(54)/dV(73) B(99) = RCT(54)*V(103) ! B(100) = dA(54)/dV(103) B(100) = RCT(54)*V(73) ! B(101) = dA(55)/dV(92) B(101) = RCT(55)*V(103) ! B(102) = dA(55)/dV(103) B(102) = RCT(55)*V(92) ! B(103) = dA(56)/dV(70) B(103) = RCT(56)*V(103) ! B(104) = dA(56)/dV(103) B(104) = RCT(56)*V(70) ! B(105) = dA(57)/dV(70) B(105) = RCT(57)*V(103) ! B(106) = dA(57)/dV(103) B(106) = RCT(57)*V(70) ! B(107) = dA(58)/dV(79) B(107) = RCT(58)*V(103) ! B(108) = dA(58)/dV(103) B(108) = RCT(58)*V(79) ! B(109) = dA(59)/dV(80) B(109) = RCT(59)*V(103) ! B(110) = dA(59)/dV(103) B(110) = RCT(59)*V(80) ! B(111) = dA(60)/dV(75) B(111) = RCT(60)*V(103) ! B(112) = dA(60)/dV(103) B(112) = RCT(60)*V(75) ! B(113) = dA(61)/dV(75) B(113) = RCT(61)*V(103) ! B(114) = dA(61)/dV(103) B(114) = RCT(61)*V(75) ! B(115) = dA(62)/dV(72) B(115) = RCT(62)*V(103) ! B(116) = dA(62)/dV(103) B(116) = RCT(62)*V(72) ! B(117) = dA(63)/dV(72) B(117) = RCT(63)*V(103) ! B(118) = dA(63)/dV(103) B(118) = RCT(63)*V(72) ! B(119) = dA(64)/dV(64) B(119) = RCT(64)*V(103) ! B(120) = dA(64)/dV(103) B(120) = RCT(64)*V(64) ! B(121) = dA(65)/dV(69) B(121) = RCT(65)*V(103) ! B(122) = dA(65)/dV(103) B(122) = RCT(65)*V(69) ! B(123) = dA(66)/dV(66) B(123) = RCT(66)*V(103) ! B(124) = dA(66)/dV(103) B(124) = RCT(66)*V(66) ! B(125) = dA(67)/dV(78) B(125) = RCT(67)*V(103) ! B(126) = dA(67)/dV(103) B(126) = RCT(67)*V(78) ! B(127) = dA(68)/dV(76) B(127) = RCT(68)*V(103) ! B(128) = dA(68)/dV(103) B(128) = RCT(68)*V(76) ! B(129) = dA(69)/dV(35) B(129) = RCT(69)*V(100) ! B(130) = dA(69)/dV(100) B(130) = RCT(69)*V(35) ! B(131) = dA(70)/dV(86) B(131) = RCT(70)*V(102) ! B(132) = dA(70)/dV(102) B(132) = RCT(70)*V(86) ! B(133) = dA(71)/dV(102) B(133) = RCT(71)*F(1) ! B(135) = dA(72)/dV(81) B(135) = RCT(72)*V(102) ! B(136) = dA(72)/dV(102) B(136) = RCT(72)*V(81) ! B(137) = dA(73)/dV(90) B(137) = RCT(73)*V(99) ! B(138) = dA(73)/dV(99) B(138) = RCT(73)*V(90) ! B(139) = dA(74)/dV(27) B(139) = RCT(74) ! B(140) = dA(75)/dV(94) B(140) = RCT(75)*V(99) ! B(141) = dA(75)/dV(99) B(141) = RCT(75)*V(94) ! B(142) = dA(76)/dV(28) B(142) = RCT(76) ! B(143) = dA(77)/dV(97) B(143) = RCT(77)*V(99) ! B(144) = dA(77)/dV(99) B(144) = RCT(77)*V(97) ! B(145) = dA(78)/dV(58) B(145) = RCT(78) ! B(146) = dA(79)/dV(99) B(146) = RCT(79)*F(5) ! B(149) = dA(81)/dV(90) B(149) = RCT(81)*V(103) ! B(150) = dA(81)/dV(103) B(150) = RCT(81)*V(90) ! B(151) = dA(82)/dV(94) B(151) = RCT(82)*V(103) ! B(152) = dA(82)/dV(103) B(152) = RCT(82)*V(94) ! B(153) = dA(83)/dV(97) B(153) = RCT(83)*V(103) ! B(154) = dA(83)/dV(103) B(154) = RCT(83)*V(97) ! B(155) = dA(84)/dV(103) B(155) = RCT(84)*F(5) ! B(157) = dA(85)/dV(81) B(157) = RCT(85)*V(100) ! B(158) = dA(85)/dV(100) B(158) = RCT(85)*V(81) ! B(159) = dA(86)/dV(62) B(159) = RCT(86)*V(102) ! B(160) = dA(86)/dV(102) B(160) = RCT(86)*V(62) ! B(161) = dA(87)/dV(62) B(161) = RCT(87)*V(102) ! B(162) = dA(87)/dV(102) B(162) = RCT(87)*V(62) ! B(163) = dA(88)/dV(65) B(163) = RCT(88)*V(95) ! B(164) = dA(88)/dV(95) B(164) = RCT(88)*V(65) ! B(165) = dA(89)/dV(83) B(165) = RCT(89)*V(95) ! B(166) = dA(89)/dV(95) B(166) = RCT(89)*V(83) ! B(167) = dA(90)/dV(85) B(167) = RCT(90)*V(96) ! B(168) = dA(90)/dV(96) B(168) = RCT(90)*V(85) ! B(169) = dA(91)/dV(68) B(169) = RCT(91)*V(96) ! B(170) = dA(91)/dV(96) B(170) = RCT(91)*V(68) ! B(171) = dA(92)/dV(74) B(171) = RCT(92)*V(96) ! B(172) = dA(92)/dV(96) B(172) = RCT(92)*V(74) ! B(173) = dA(93)/dV(73) B(173) = RCT(93)*V(96) ! B(174) = dA(93)/dV(96) B(174) = RCT(93)*V(73) ! B(175) = dA(94)/dV(92) B(175) = RCT(94)*V(96) ! B(176) = dA(94)/dV(96) B(176) = RCT(94)*V(92) ! B(177) = dA(95)/dV(70) B(177) = RCT(95)*V(96) ! B(178) = dA(95)/dV(96) B(178) = RCT(95)*V(70) ! B(179) = dA(96)/dV(79) B(179) = RCT(96)*V(96) ! B(180) = dA(96)/dV(96) B(180) = RCT(96)*V(79) ! B(181) = dA(97)/dV(80) B(181) = RCT(97)*V(96) ! B(182) = dA(97)/dV(96) B(182) = RCT(97)*V(80) ! B(183) = dA(98)/dV(75) B(183) = RCT(98)*V(96) ! B(184) = dA(98)/dV(96) B(184) = RCT(98)*V(75) ! B(185) = dA(99)/dV(72) B(185) = RCT(99)*V(96) ! B(186) = dA(99)/dV(96) B(186) = RCT(99)*V(72) ! B(187) = dA(100)/dV(64) B(187) = RCT(100)*V(96) ! B(188) = dA(100)/dV(96) B(188) = RCT(100)*V(64) ! B(189) = dA(101)/dV(69) B(189) = RCT(101)*V(96) ! B(190) = dA(101)/dV(96) B(190) = RCT(101)*V(69) ! B(191) = dA(102)/dV(66) B(191) = RCT(102)*V(96) ! B(192) = dA(102)/dV(96) B(192) = RCT(102)*V(66) ! B(193) = dA(103)/dV(78) B(193) = RCT(103)*V(96) ! B(194) = dA(103)/dV(96) B(194) = RCT(103)*V(78) ! B(195) = dA(104)/dV(76) B(195) = RCT(104)*V(96) ! B(196) = dA(104)/dV(96) B(196) = RCT(104)*V(76) ! B(197) = dA(105)/dV(89) B(197) = RCT(105)*V(102) ! B(198) = dA(105)/dV(102) B(198) = RCT(105)*V(89) ! B(199) = dA(106)/dV(87) B(199) = RCT(106)*V(95) ! B(200) = dA(106)/dV(95) B(200) = RCT(106)*V(87) ! B(201) = dA(107)/dV(89) B(201) = RCT(107)*V(100) ! B(202) = dA(107)/dV(100) B(202) = RCT(107)*V(89) ! B(203) = dA(108)/dV(85) B(203) = RCT(108)*V(95) ! B(204) = dA(108)/dV(95) B(204) = RCT(108)*V(85) ! B(205) = dA(109)/dV(68) B(205) = RCT(109)*V(95) ! B(206) = dA(109)/dV(95) B(206) = RCT(109)*V(68) ! B(207) = dA(110)/dV(74) B(207) = RCT(110)*V(95) ! B(208) = dA(110)/dV(95) B(208) = RCT(110)*V(74) ! B(209) = dA(111)/dV(73) B(209) = RCT(111)*V(95) ! B(210) = dA(111)/dV(95) B(210) = RCT(111)*V(73) ! B(211) = dA(112)/dV(92) B(211) = RCT(112)*V(95) ! B(212) = dA(112)/dV(95) B(212) = RCT(112)*V(92) ! B(213) = dA(113)/dV(70) B(213) = RCT(113)*V(95) ! B(214) = dA(113)/dV(95) B(214) = RCT(113)*V(70) ! B(215) = dA(114)/dV(79) B(215) = RCT(114)*V(95) ! B(216) = dA(114)/dV(95) B(216) = RCT(114)*V(79) ! B(217) = dA(115)/dV(80) B(217) = RCT(115)*V(95) ! B(218) = dA(115)/dV(95) B(218) = RCT(115)*V(80) ! B(219) = dA(116)/dV(75) B(219) = RCT(116)*V(95) ! B(220) = dA(116)/dV(95) B(220) = RCT(116)*V(75) ! B(221) = dA(117)/dV(72) B(221) = RCT(117)*V(95) ! B(222) = dA(117)/dV(95) B(222) = RCT(117)*V(72) ! B(223) = dA(118)/dV(64) B(223) = RCT(118)*V(95) ! B(224) = dA(118)/dV(95) B(224) = RCT(118)*V(64) ! B(225) = dA(119)/dV(69) B(225) = RCT(119)*V(95) ! B(226) = dA(119)/dV(95) B(226) = RCT(119)*V(69) ! B(227) = dA(120)/dV(66) B(227) = RCT(120)*V(95) ! B(228) = dA(120)/dV(95) B(228) = RCT(120)*V(66) ! B(229) = dA(121)/dV(78) B(229) = RCT(121)*V(95) ! B(230) = dA(121)/dV(95) B(230) = RCT(121)*V(78) ! B(231) = dA(122)/dV(76) B(231) = RCT(122)*V(95) ! B(232) = dA(122)/dV(95) B(232) = RCT(122)*V(76) ! B(233) = dA(123)/dV(102) B(233) = RCT(123)*F(4) ! B(235) = dA(124)/dV(102) B(235) = RCT(124)*F(16) ! B(237) = dA(125)/dV(87) B(237) = RCT(125)*2*V(87) ! B(238) = dA(126)/dV(87) B(238) = RCT(126)*2*V(87) ! B(239) = dA(127)/dV(87) B(239) = RCT(127)*V(96) ! B(240) = dA(127)/dV(96) B(240) = RCT(127)*V(87) ! B(241) = dA(128)/dV(65) B(241) = RCT(128)*V(96) ! B(242) = dA(128)/dV(96) B(242) = RCT(128)*V(65) ! B(243) = dA(129)/dV(83) B(243) = RCT(129)*V(96) ! B(244) = dA(129)/dV(96) B(244) = RCT(129)*V(83) ! B(245) = dA(130)/dV(96) B(245) = RCT(130)*V(101) ! B(246) = dA(130)/dV(101) B(246) = RCT(130)*V(96) ! B(247) = dA(131)/dV(90) B(247) = RCT(131)*V(96) ! B(248) = dA(131)/dV(96) B(248) = RCT(131)*V(90) ! B(249) = dA(132)/dV(94) B(249) = RCT(132)*V(96) ! B(250) = dA(132)/dV(96) B(250) = RCT(132)*V(94) ! B(251) = dA(133)/dV(96) B(251) = RCT(133)*V(97) ! B(252) = dA(133)/dV(97) B(252) = RCT(133)*V(96) ! B(253) = dA(134)/dV(96) B(253) = RCT(134)*F(5) ! B(255) = dA(135)/dV(61) B(255) = RCT(135)*V(102) ! B(256) = dA(135)/dV(102) B(256) = RCT(135)*V(61) ! B(257) = dA(136)/dV(61) B(257) = RCT(136)*V(98) ! B(258) = dA(136)/dV(98) B(258) = RCT(136)*V(61) ! B(259) = dA(137)/dV(58) B(259) = RCT(137)*V(102) ! B(260) = dA(137)/dV(102) B(260) = RCT(137)*V(58) ! B(261) = dA(138)/dV(58) B(261) = RCT(138)*V(98) ! B(262) = dA(138)/dV(98) B(262) = RCT(138)*V(58) ! B(263) = dA(139)/dV(63) B(263) = RCT(139)*V(102) ! B(264) = dA(139)/dV(102) B(264) = RCT(139)*V(63) ! B(265) = dA(140)/dV(61) B(265) = RCT(140)*V(100) ! B(266) = dA(140)/dV(100) B(266) = RCT(140)*V(61) ! B(267) = dA(141)/dV(102) B(267) = RCT(141)*F(8) ! B(269) = dA(142)/dV(88) B(269) = RCT(142)*V(102) ! B(270) = dA(142)/dV(102) B(270) = RCT(142)*V(88) ! B(271) = dA(143)/dV(100) B(271) = RCT(143)*F(8) ! B(273) = dA(144)/dV(88) B(273) = RCT(144)*V(100) ! B(274) = dA(144)/dV(100) B(274) = RCT(144)*V(88) ! B(275) = dA(145)/dV(59) B(275) = RCT(145)*V(102) ! B(276) = dA(145)/dV(102) B(276) = RCT(145)*V(59) ! B(277) = dA(146)/dV(93) B(277) = RCT(146)*V(102) ! B(278) = dA(146)/dV(102) B(278) = RCT(146)*V(93) ! B(279) = dA(147)/dV(91) B(279) = RCT(147)*V(102) ! B(280) = dA(147)/dV(102) B(280) = RCT(147)*V(91) ! B(281) = dA(148)/dV(77) B(281) = RCT(148)*V(102) ! B(282) = dA(148)/dV(102) B(282) = RCT(148)*V(77) ! B(283) = dA(149)/dV(65) B(283) = RCT(149)*V(101) ! B(284) = dA(149)/dV(101) B(284) = RCT(149)*V(65) ! B(285) = dA(150)/dV(83) B(285) = RCT(150)*V(101) ! B(286) = dA(150)/dV(101) B(286) = RCT(150)*V(83) ! B(287) = dA(151)/dV(65) B(287) = RCT(151)*V(101) ! B(288) = dA(151)/dV(101) B(288) = RCT(151)*V(65) ! B(289) = dA(152)/dV(83) B(289) = RCT(152)*V(101) ! B(290) = dA(152)/dV(101) B(290) = RCT(152)*V(83) ! B(291) = dA(153)/dV(59) B(291) = RCT(153)*V(98) ! B(292) = dA(153)/dV(98) B(292) = RCT(153)*V(59) ! B(293) = dA(154)/dV(93) B(293) = RCT(154)*V(98) ! B(294) = dA(154)/dV(98) B(294) = RCT(154)*V(93) ! B(295) = dA(155)/dV(91) B(295) = RCT(155)*V(98) ! B(296) = dA(155)/dV(98) B(296) = RCT(155)*V(91) ! B(297) = dA(156)/dV(59) B(297) = RCT(156)*V(100) ! B(298) = dA(156)/dV(100) B(298) = RCT(156)*V(59) ! B(299) = dA(157)/dV(91) B(299) = RCT(157)*V(100) ! B(300) = dA(157)/dV(100) B(300) = RCT(157)*V(91) ! B(301) = dA(158)/dV(91) B(301) = RCT(158)*V(100) ! B(302) = dA(158)/dV(100) B(302) = RCT(158)*V(91) ! B(303) = dA(159)/dV(90) B(303) = RCT(159)*V(95) ! B(304) = dA(159)/dV(95) B(304) = RCT(159)*V(90) ! B(305) = dA(160)/dV(94) B(305) = RCT(160)*V(95) ! B(306) = dA(160)/dV(95) B(306) = RCT(160)*V(94) ! B(307) = dA(161)/dV(95) B(307) = RCT(161)*V(97) ! B(308) = dA(161)/dV(97) B(308) = RCT(161)*V(95) ! B(309) = dA(162)/dV(95) B(309) = RCT(162)*F(5) ! B(311) = dA(163)/dV(90) B(311) = RCT(163)*V(95) ! B(312) = dA(163)/dV(95) B(312) = RCT(163)*V(90) ! B(313) = dA(164)/dV(94) B(313) = RCT(164)*V(95) ! B(314) = dA(164)/dV(95) B(314) = RCT(164)*V(94) ! B(315) = dA(165)/dV(95) B(315) = RCT(165)*V(97) ! B(316) = dA(165)/dV(97) B(316) = RCT(165)*V(95) ! B(317) = dA(166)/dV(95) B(317) = RCT(166)*F(5) ! B(319) = dA(167)/dV(52) B(319) = RCT(167)*V(102) ! B(320) = dA(167)/dV(102) B(320) = RCT(167)*V(52) ! B(321) = dA(168)/dV(51) B(321) = RCT(168)*V(102) ! B(322) = dA(168)/dV(102) B(322) = RCT(168)*V(51) ! B(323) = dA(169)/dV(46) B(323) = RCT(169)*V(102) ! B(324) = dA(169)/dV(102) B(324) = RCT(169)*V(46) ! B(325) = dA(170)/dV(41) B(325) = RCT(170)*V(102) ! B(326) = dA(170)/dV(102) B(326) = RCT(170)*V(41) ! B(327) = dA(171)/dV(42) B(327) = RCT(171)*V(102) ! B(328) = dA(171)/dV(102) B(328) = RCT(171)*V(42) ! B(329) = dA(172)/dV(40) B(329) = RCT(172)*V(102) ! B(330) = dA(172)/dV(102) B(330) = RCT(172)*V(40) ! B(331) = dA(173)/dV(43) B(331) = RCT(173)*V(102) ! B(332) = dA(173)/dV(102) B(332) = RCT(173)*V(43) ! B(333) = dA(174)/dV(50) B(333) = RCT(174)*V(102) ! B(334) = dA(174)/dV(102) B(334) = RCT(174)*V(50) ! B(335) = dA(175)/dV(47) B(335) = RCT(175)*V(102) ! B(336) = dA(175)/dV(102) B(336) = RCT(175)*V(47) ! B(337) = dA(176)/dV(102) B(337) = RCT(176)*F(6) ! B(339) = dA(177)/dV(57) B(339) = RCT(177)*V(102) ! B(340) = dA(177)/dV(102) B(340) = RCT(177)*V(57) ! B(341) = dA(178)/dV(54) B(341) = RCT(178)*V(102) ! B(342) = dA(178)/dV(102) B(342) = RCT(178)*V(54) ! B(343) = dA(179)/dV(56) B(343) = RCT(179)*V(102) ! B(344) = dA(179)/dV(102) B(344) = RCT(179)*V(56) ! B(345) = dA(180)/dV(55) B(345) = RCT(180)*V(102) ! B(346) = dA(180)/dV(102) B(346) = RCT(180)*V(55) ! B(347) = dA(181)/dV(53) B(347) = RCT(181)*V(102) ! B(348) = dA(181)/dV(102) B(348) = RCT(181)*V(53) ! B(349) = dA(182)/dV(45) B(349) = RCT(182)*V(102) ! B(350) = dA(182)/dV(102) B(350) = RCT(182)*V(45) ! B(351) = dA(183)/dV(48) B(351) = RCT(183)*V(102) ! B(352) = dA(183)/dV(102) B(352) = RCT(183)*V(48) ! B(353) = dA(184)/dV(36) B(353) = RCT(184)*V(100) ! B(354) = dA(184)/dV(100) B(354) = RCT(184)*V(36) ! B(355) = dA(185)/dV(102) B(355) = RCT(185)*F(12) ! B(357) = dA(186)/dV(71) B(357) = RCT(186)*V(102) ! B(358) = dA(186)/dV(102) B(358) = RCT(186)*V(71) ! B(359) = dA(187)/dV(71) B(359) = RCT(187)*V(98) ! B(360) = dA(187)/dV(98) B(360) = RCT(187)*V(71) ! B(361) = dA(188)/dV(101) B(361) = RCT(188)*2*V(101) ! B(362) = dA(189)/dV(95) B(362) = RCT(189)*V(101) ! B(363) = dA(189)/dV(101) B(363) = RCT(189)*V(95) ! B(364) = dA(190)/dV(95) B(364) = RCT(190)*V(101) ! B(365) = dA(190)/dV(101) B(365) = RCT(190)*V(95) ! B(366) = dA(191)/dV(85) B(366) = RCT(191)*V(101) ! B(367) = dA(191)/dV(101) B(367) = RCT(191)*V(85) ! B(368) = dA(192)/dV(74) B(368) = RCT(192)*V(101) ! B(369) = dA(192)/dV(101) B(369) = RCT(192)*V(74) ! B(370) = dA(193)/dV(73) B(370) = RCT(193)*V(101) ! B(371) = dA(193)/dV(101) B(371) = RCT(193)*V(73) ! B(372) = dA(194)/dV(92) B(372) = RCT(194)*V(101) ! B(373) = dA(194)/dV(101) B(373) = RCT(194)*V(92) ! B(374) = dA(195)/dV(70) B(374) = RCT(195)*V(101) ! B(375) = dA(195)/dV(101) B(375) = RCT(195)*V(70) ! B(376) = dA(196)/dV(79) B(376) = RCT(196)*V(101) ! B(377) = dA(196)/dV(101) B(377) = RCT(196)*V(79) ! B(378) = dA(197)/dV(80) B(378) = RCT(197)*V(101) ! B(379) = dA(197)/dV(101) B(379) = RCT(197)*V(80) ! B(380) = dA(198)/dV(75) B(380) = RCT(198)*V(101) ! B(381) = dA(198)/dV(101) B(381) = RCT(198)*V(75) ! B(382) = dA(199)/dV(72) B(382) = RCT(199)*V(101) ! B(383) = dA(199)/dV(101) B(383) = RCT(199)*V(72) ! B(384) = dA(200)/dV(66) B(384) = RCT(200)*V(101) ! B(385) = dA(200)/dV(101) B(385) = RCT(200)*V(66) ! B(386) = dA(201)/dV(68) B(386) = RCT(201)*V(101) ! B(387) = dA(201)/dV(101) B(387) = RCT(201)*V(68) ! B(388) = dA(202)/dV(64) B(388) = RCT(202)*V(101) ! B(389) = dA(202)/dV(101) B(389) = RCT(202)*V(64) ! B(390) = dA(203)/dV(69) B(390) = RCT(203)*V(101) ! B(391) = dA(203)/dV(101) B(391) = RCT(203)*V(69) ! B(392) = dA(204)/dV(78) B(392) = RCT(204)*V(101) ! B(393) = dA(204)/dV(101) B(393) = RCT(204)*V(78) ! B(394) = dA(205)/dV(76) B(394) = RCT(205)*V(101) ! B(395) = dA(205)/dV(101) B(395) = RCT(205)*V(76) ! B(396) = dA(206)/dV(85) B(396) = RCT(206)*V(101) ! B(397) = dA(206)/dV(101) B(397) = RCT(206)*V(85) ! B(398) = dA(207)/dV(74) B(398) = RCT(207)*V(101) ! B(399) = dA(207)/dV(101) B(399) = RCT(207)*V(74) ! B(400) = dA(208)/dV(73) B(400) = RCT(208)*V(101) ! B(401) = dA(208)/dV(101) B(401) = RCT(208)*V(73) ! B(402) = dA(209)/dV(92) B(402) = RCT(209)*V(101) ! B(403) = dA(209)/dV(101) B(403) = RCT(209)*V(92) ! B(404) = dA(210)/dV(70) B(404) = RCT(210)*V(101) ! B(405) = dA(210)/dV(101) B(405) = RCT(210)*V(70) ! B(406) = dA(211)/dV(79) B(406) = RCT(211)*V(101) ! B(407) = dA(211)/dV(101) B(407) = RCT(211)*V(79) ! B(408) = dA(212)/dV(75) B(408) = RCT(212)*V(101) ! B(409) = dA(212)/dV(101) B(409) = RCT(212)*V(75) ! B(410) = dA(213)/dV(72) B(410) = RCT(213)*V(101) ! B(411) = dA(213)/dV(101) B(411) = RCT(213)*V(72) ! B(412) = dA(214)/dV(68) B(412) = RCT(214)*V(101) ! B(413) = dA(214)/dV(101) B(413) = RCT(214)*V(68) ! B(414) = dA(215)/dV(80) B(414) = RCT(215)*V(101) ! B(415) = dA(215)/dV(101) B(415) = RCT(215)*V(80) ! B(416) = dA(216)/dV(64) B(416) = RCT(216)*V(101) ! B(417) = dA(216)/dV(101) B(417) = RCT(216)*V(64) ! B(418) = dA(217)/dV(69) B(418) = RCT(217)*V(101) ! B(419) = dA(217)/dV(101) B(419) = RCT(217)*V(69) ! B(420) = dA(218)/dV(78) B(420) = RCT(218)*V(101) ! B(421) = dA(218)/dV(101) B(421) = RCT(218)*V(78) ! B(422) = dA(219)/dV(76) B(422) = RCT(219)*V(101) ! B(423) = dA(219)/dV(101) B(423) = RCT(219)*V(76) ! B(424) = dA(220)/dV(66) B(424) = RCT(220)*V(101) ! B(425) = dA(220)/dV(101) B(425) = RCT(220)*V(66) ! B(426) = dA(221)/dV(87) B(426) = RCT(221)*V(101) ! B(427) = dA(221)/dV(101) B(427) = RCT(221)*V(87) ! B(428) = dA(222)/dV(87) B(428) = RCT(222)*V(101) ! B(429) = dA(222)/dV(101) B(429) = RCT(222)*V(87) ! B(430) = dA(223)/dV(90) B(430) = RCT(223)*V(101) ! B(431) = dA(223)/dV(101) B(431) = RCT(223)*V(90) ! B(432) = dA(224)/dV(94) B(432) = RCT(224)*V(101) ! B(433) = dA(224)/dV(101) B(433) = RCT(224)*V(94) ! B(434) = dA(225)/dV(97) B(434) = RCT(225)*V(101) ! B(435) = dA(225)/dV(101) B(435) = RCT(225)*V(97) ! B(436) = dA(226)/dV(101) B(436) = RCT(226)*F(5) ! B(438) = dA(227)/dV(100) B(438) = RCT(227)*2*V(100) ! B(457) = dA(246)/dV(99) B(457) = RCT(246) ! B(458) = dA(247)/dV(98) B(458) = RCT(247) ! B(459) = dA(248)/dV(34) B(459) = RCT(248) ! B(460) = dA(249)/dV(67) B(460) = RCT(249) ! B(461) = dA(250)/dV(82) B(461) = RCT(250) ! B(462) = dA(251)/dV(38) B(462) = RCT(251) ! B(463) = dA(252)/dV(29) B(463) = RCT(252) ! B(464) = dA(253)/dV(58) B(464) = RCT(253) ! B(465) = dA(254)/dV(27) B(465) = RCT(254) ! B(466) = dA(255)/dV(86) B(466) = RCT(255) ! B(467) = dA(256)/dV(96) B(467) = RCT(256) ! B(468) = dA(257)/dV(99) B(468) = RCT(257) ! B(469) = dA(258)/dV(100) B(469) = RCT(258) ! B(470) = dA(259)/dV(38) B(470) = RCT(259) ! B(471) = dA(260)/dV(44) B(471) = RCT(260)*V(102) ! B(472) = dA(260)/dV(102) B(472) = RCT(260)*V(44) ! B(473) = dA(261)/dV(44) B(473) = RCT(261)*V(102) ! B(474) = dA(261)/dV(102) B(474) = RCT(261)*V(44) ! B(475) = dA(262)/dV(44) B(475) = RCT(262)*V(100) ! B(476) = dA(262)/dV(100) B(476) = RCT(262)*V(44) ! B(477) = dA(263)/dV(30) B(477) = RCT(263)*V(102) ! B(478) = dA(263)/dV(102) B(478) = RCT(263)*V(30) ! B(479) = dA(264)/dV(25) B(479) = RCT(264)*V(102) ! B(480) = dA(264)/dV(102) B(480) = RCT(264)*V(25) ! B(481) = dA(265)/dV(23) B(481) = RCT(265)*V(102) ! B(482) = dA(265)/dV(102) B(482) = RCT(265)*V(23) ! B(483) = dA(266)/dV(24) B(483) = RCT(266)*V(102) ! B(484) = dA(266)/dV(102) B(484) = RCT(266)*V(24) ! B(485) = dA(267)/dV(1) B(485) = RCT(267) ! B(486) = dA(268)/dV(31) B(486) = RCT(268)*V(96) ! B(487) = dA(268)/dV(96) B(487) = RCT(268)*V(31) ! B(488) = dA(269)/dV(31) B(488) = RCT(269)*V(103) ! B(489) = dA(269)/dV(103) B(489) = RCT(269)*V(31) ! B(490) = dA(270)/dV(32) B(490) = RCT(270)*V(96) ! B(491) = dA(270)/dV(96) B(491) = RCT(270)*V(32) ! B(492) = dA(271)/dV(32) B(492) = RCT(271)*V(103) ! B(493) = dA(271)/dV(103) B(493) = RCT(271)*V(32) ! B(494) = dA(272)/dV(33) B(494) = RCT(272)*V(96) ! B(495) = dA(272)/dV(96) B(495) = RCT(272)*V(33) ! B(496) = dA(273)/dV(33) B(496) = RCT(273)*V(103) ! B(497) = dA(273)/dV(103) B(497) = RCT(273)*V(33) ! B(498) = dA(274)/dV(98) B(498) = RCT(274) ! B(499) = dA(275)/dV(99) B(499) = RCT(275) ! B(500) = dA(276)/dV(29) B(500) = RCT(276) ! B(501) = dA(277)/dV(39) B(501) = RCT(277) ! B(502) = dA(278)/dV(82) B(502) = RCT(278) ! B(503) = dA(279)/dV(82) B(503) = RCT(279) ! B(504) = dA(280)/dV(67) B(504) = RCT(280) ! B(505) = dA(281)/dV(37) B(505) = RCT(281) ! B(506) = dA(282)/dV(49) B(506) = RCT(282) ! B(507) = dA(283)/dV(100) B(507) = RCT(283) ! B(508) = dA(284)/dV(100) B(508) = RCT(284) ! B(509) = dA(285)/dV(38) B(509) = RCT(285) ! B(510) = dA(286)/dV(38) B(510) = RCT(286) ! B(511) = dA(287)/dV(49) B(511) = RCT(287) ! B(512) = dA(288)/dV(84) B(512) = RCT(288) ! B(513) = dA(289)/dV(84) B(513) = RCT(289) ! B(514) = dA(290)/dV(34) B(514) = RCT(290) ! B(515) = dA(291)/dV(81) B(515) = RCT(291) ! B(516) = dA(292)/dV(62) B(516) = RCT(292) ! B(517) = dA(293)/dV(62) B(517) = RCT(293) ! B(518) = dA(294)/dV(89) B(518) = RCT(294) ! B(520) = dA(296)/dV(63) B(520) = RCT(296) ! B(523) = dA(299)/dV(88) B(523) = RCT(299) ! B(524) = dA(300)/dV(88) B(524) = RCT(300) ! B(525) = dA(301)/dV(93) B(525) = RCT(301) ! B(526) = dA(302)/dV(93) B(526) = RCT(302) ! B(527) = dA(303)/dV(93) B(527) = RCT(303) ! B(528) = dA(304)/dV(91) B(528) = RCT(304) ! B(529) = dA(305)/dV(91) B(529) = RCT(305) ! B(530) = dA(306)/dV(77) B(530) = RCT(306) ! B(531) = dA(307)/dV(52) B(531) = RCT(307) ! B(532) = dA(308)/dV(51) B(532) = RCT(308) ! B(533) = dA(309)/dV(46) B(533) = RCT(309) ! B(534) = dA(310)/dV(41) B(534) = RCT(310) ! B(535) = dA(311)/dV(42) B(535) = RCT(311) ! B(536) = dA(312)/dV(40) B(536) = RCT(312) ! B(537) = dA(313)/dV(50) B(537) = RCT(313) ! B(538) = dA(314)/dV(43) B(538) = RCT(314) ! B(539) = dA(315)/dV(47) B(539) = RCT(315) ! B(541) = dA(317)/dV(57) B(541) = RCT(317) ! B(542) = dA(318)/dV(54) B(542) = RCT(318) ! B(543) = dA(319)/dV(56) B(543) = RCT(319) ! B(544) = dA(320)/dV(55) B(544) = RCT(320) ! B(545) = dA(321)/dV(53) B(545) = RCT(321) ! B(546) = dA(322)/dV(45) B(546) = RCT(322) ! B(547) = dA(323)/dV(86) B(547) = RCT(323) ! B(548) = dA(324)/dV(48) B(548) = RCT(324) ! Construct the Jacobian terms from B's ! JVS(1) = Jac_FULL(1,1) JVS(1) = -B(485) ! JVS(2) = Jac_FULL(2,2) JVS(2) = 0 ! JVS(3) = Jac_FULL(2,82) JVS(3) = B(461) ! JVS(4) = Jac_FULL(3,3) JVS(4) = 0 ! JVS(5) = Jac_FULL(3,29) JVS(5) = B(463) ! JVS(6) = Jac_FULL(4,4) JVS(6) = 0 ! JVS(7) = Jac_FULL(4,67) JVS(7) = B(460) ! JVS(8) = Jac_FULL(5,5) JVS(8) = 0 ! JVS(9) = Jac_FULL(5,38) JVS(9) = B(462) ! JVS(10) = Jac_FULL(6,6) JVS(10) = 0 ! JVS(11) = Jac_FULL(6,99) JVS(11) = B(457) ! JVS(12) = Jac_FULL(7,7) JVS(12) = 0 ! JVS(13) = Jac_FULL(7,98) JVS(13) = B(458) ! JVS(14) = Jac_FULL(8,8) JVS(14) = 0 ! JVS(15) = Jac_FULL(8,34) JVS(15) = B(459) ! JVS(16) = Jac_FULL(9,9) JVS(16) = 0 ! JVS(17) = Jac_FULL(9,58) JVS(17) = B(464) ! JVS(18) = Jac_FULL(10,10) JVS(18) = 0 ! JVS(19) = Jac_FULL(10,27) JVS(19) = B(465) ! JVS(20) = Jac_FULL(11,11) JVS(20) = 0 ! JVS(21) = Jac_FULL(11,86) JVS(21) = B(466) ! JVS(22) = Jac_FULL(12,12) JVS(22) = 0 ! JVS(23) = Jac_FULL(12,30) JVS(23) = B(477) ! JVS(24) = Jac_FULL(12,102) JVS(24) = B(478) ! JVS(25) = Jac_FULL(13,13) JVS(25) = 0 ! JVS(26) = Jac_FULL(13,44) JVS(26) = 0.25*B(473) ! JVS(27) = Jac_FULL(13,102) JVS(27) = 0.25*B(474) ! JVS(28) = Jac_FULL(14,14) JVS(28) = 0 ! JVS(29) = Jac_FULL(14,59) JVS(29) = B(275) ! JVS(30) = Jac_FULL(14,102) JVS(30) = B(276) ! JVS(31) = Jac_FULL(15,15) JVS(31) = 0 ! JVS(32) = Jac_FULL(15,31) JVS(32) = B(486) ! JVS(33) = Jac_FULL(15,96) JVS(33) = B(487) ! JVS(34) = Jac_FULL(16,16) JVS(34) = 0 ! JVS(35) = Jac_FULL(16,31) JVS(35) = B(488) ! JVS(36) = Jac_FULL(16,103) JVS(36) = B(489) ! JVS(37) = Jac_FULL(17,17) JVS(37) = 0 ! JVS(38) = Jac_FULL(17,32) JVS(38) = B(490) ! JVS(39) = Jac_FULL(17,96) JVS(39) = B(491) ! JVS(40) = Jac_FULL(18,18) JVS(40) = 0 ! JVS(41) = Jac_FULL(18,32) JVS(41) = B(492) ! JVS(42) = Jac_FULL(18,103) JVS(42) = B(493) ! JVS(43) = Jac_FULL(19,19) JVS(43) = 0 ! JVS(44) = Jac_FULL(19,33) JVS(44) = B(494) ! JVS(45) = Jac_FULL(19,96) JVS(45) = B(495) ! JVS(46) = Jac_FULL(20,20) JVS(46) = 0 ! JVS(47) = Jac_FULL(20,33) JVS(47) = B(496) ! JVS(48) = Jac_FULL(20,103) JVS(48) = B(497) ! JVS(49) = Jac_FULL(21,21) JVS(49) = 0 ! JVS(50) = Jac_FULL(21,59) JVS(50) = 0.15*B(291) ! JVS(51) = Jac_FULL(21,60) JVS(51) = B(22) ! JVS(52) = Jac_FULL(21,91) JVS(52) = 0.16*B(295) ! JVS(53) = Jac_FULL(21,98) JVS(53) = 0.15*B(292)+0.16*B(296) ! JVS(54) = Jac_FULL(21,101) JVS(54) = B(75) ! JVS(55) = Jac_FULL(21,102) JVS(55) = B(23)+B(60)+B(133) ! JVS(56) = Jac_FULL(21,103) JVS(56) = B(76) ! JVS(57) = Jac_FULL(22,22) JVS(57) = 0 ! JVS(58) = Jac_FULL(22,27) JVS(58) = B(465) ! JVS(59) = Jac_FULL(22,29) JVS(59) = B(463) ! JVS(60) = Jac_FULL(22,34) JVS(60) = B(459) ! JVS(61) = Jac_FULL(22,38) JVS(61) = B(462) ! JVS(62) = Jac_FULL(22,58) JVS(62) = B(464) ! JVS(63) = Jac_FULL(22,67) JVS(63) = B(460) ! JVS(64) = Jac_FULL(22,82) JVS(64) = B(461) ! JVS(65) = Jac_FULL(22,86) JVS(65) = B(466) ! JVS(66) = Jac_FULL(22,98) JVS(66) = B(458) ! JVS(67) = Jac_FULL(22,99) JVS(67) = B(457) ! JVS(68) = Jac_FULL(23,23) JVS(68) = -B(481) ! JVS(69) = Jac_FULL(23,102) JVS(69) = -B(482) ! JVS(70) = Jac_FULL(24,24) JVS(70) = -B(483) ! JVS(71) = Jac_FULL(24,102) JVS(71) = -B(484) ! JVS(72) = Jac_FULL(25,25) JVS(72) = -B(479) ! JVS(73) = Jac_FULL(25,102) JVS(73) = -B(480) ! JVS(74) = Jac_FULL(26,26) JVS(74) = -B(81)-B(83) ! JVS(75) = Jac_FULL(26,102) JVS(75) = -B(82)-B(84) ! JVS(76) = Jac_FULL(27,27) JVS(76) = -B(139)-B(465) ! JVS(77) = Jac_FULL(27,90) JVS(77) = B(137) ! JVS(78) = Jac_FULL(27,99) JVS(78) = B(138) ! JVS(79) = Jac_FULL(28,28) JVS(79) = -B(142) ! JVS(80) = Jac_FULL(28,94) JVS(80) = B(140) ! JVS(81) = Jac_FULL(28,99) JVS(81) = B(141) ! JVS(82) = Jac_FULL(29,29) JVS(82) = -B(15)-B(463)-B(500) ! JVS(83) = Jac_FULL(29,96) JVS(83) = B(19)+0.5*B(467) ! JVS(84) = Jac_FULL(29,102) JVS(84) = B(12)-B(16) ! JVS(85) = Jac_FULL(30,30) JVS(85) = -B(477) ! JVS(86) = Jac_FULL(30,44) JVS(86) = B(471)+0.75*B(473)+B(475) ! JVS(87) = Jac_FULL(30,100) JVS(87) = B(476) ! JVS(88) = Jac_FULL(30,102) JVS(88) = B(472)+0.75*B(474)-B(478) ! JVS(89) = Jac_FULL(31,25) JVS(89) = B(479) ! JVS(90) = Jac_FULL(31,31) JVS(90) = -B(486)-B(488) ! JVS(91) = Jac_FULL(31,96) JVS(91) = -B(487) ! JVS(92) = Jac_FULL(31,102) JVS(92) = B(480) ! JVS(93) = Jac_FULL(31,103) JVS(93) = -B(489) ! JVS(94) = Jac_FULL(32,23) JVS(94) = B(481) ! JVS(95) = Jac_FULL(32,32) JVS(95) = -B(490)-B(492) ! JVS(96) = Jac_FULL(32,96) JVS(96) = -B(491) ! JVS(97) = Jac_FULL(32,102) JVS(97) = B(482) ! JVS(98) = Jac_FULL(32,103) JVS(98) = -B(493) ! JVS(99) = Jac_FULL(33,24) JVS(99) = B(483) ! JVS(100) = Jac_FULL(33,33) JVS(100) = -B(494)-B(496) ! JVS(101) = Jac_FULL(33,96) JVS(101) = -B(495) ! JVS(102) = Jac_FULL(33,102) JVS(102) = B(484) ! JVS(103) = Jac_FULL(33,103) JVS(103) = -B(497) ! JVS(104) = Jac_FULL(34,34) JVS(104) = -B(74)-B(459)-B(514) ! JVS(105) = Jac_FULL(34,99) JVS(105) = B(72) ! JVS(106) = Jac_FULL(34,101) JVS(106) = B(73) ! JVS(107) = Jac_FULL(35,35) JVS(107) = -B(89)-B(129) ! JVS(108) = Jac_FULL(35,100) JVS(108) = -B(130) ! JVS(109) = Jac_FULL(35,102) JVS(109) = -B(90) ! JVS(110) = Jac_FULL(36,36) JVS(110) = -B(77)-B(353) ! JVS(111) = Jac_FULL(36,100) JVS(111) = -B(354) ! JVS(112) = Jac_FULL(36,102) JVS(112) = -B(78) ! JVS(113) = Jac_FULL(37,37) JVS(113) = -B(44)-B(505) ! JVS(114) = Jac_FULL(37,99) JVS(114) = 0.5*B(468) ! JVS(115) = Jac_FULL(37,102) JVS(115) = B(42)-B(45) ! JVS(116) = Jac_FULL(37,103) JVS(116) = B(43) ! JVS(117) = Jac_FULL(38,38) JVS(117) = -B(59)-B(462)-B(470)-B(509)-B(510) ! JVS(118) = Jac_FULL(38,99) JVS(118) = B(57) ! JVS(119) = Jac_FULL(38,100) JVS(119) = B(58) ! JVS(120) = Jac_FULL(39,39) JVS(120) = -B(32)-B(34)-B(501) ! JVS(121) = Jac_FULL(39,95) JVS(121) = B(28) ! JVS(122) = Jac_FULL(39,96) JVS(122) = B(29) ! JVS(123) = Jac_FULL(39,102) JVS(123) = -B(33)-B(35) ! JVS(124) = Jac_FULL(40,40) JVS(124) = -B(329)-B(536) ! JVS(125) = Jac_FULL(40,85) JVS(125) = B(167) ! JVS(126) = Jac_FULL(40,96) JVS(126) = B(168) ! JVS(127) = Jac_FULL(40,102) JVS(127) = -B(330) ! JVS(128) = Jac_FULL(41,41) JVS(128) = -B(325)-B(534) ! JVS(129) = Jac_FULL(41,65) JVS(129) = B(241) ! JVS(130) = Jac_FULL(41,96) JVS(130) = B(242) ! JVS(131) = Jac_FULL(41,102) JVS(131) = -B(326) ! JVS(132) = Jac_FULL(42,42) JVS(132) = -B(327)-B(535) ! JVS(133) = Jac_FULL(42,66) JVS(133) = B(191) ! JVS(134) = Jac_FULL(42,96) JVS(134) = B(192) ! JVS(135) = Jac_FULL(42,102) JVS(135) = -B(328) ! JVS(136) = Jac_FULL(43,43) JVS(136) = -B(331)-B(538) ! JVS(137) = Jac_FULL(43,90) JVS(137) = 0.7*B(247) ! JVS(138) = Jac_FULL(43,96) JVS(138) = 0.7*B(248) ! JVS(139) = Jac_FULL(43,102) JVS(139) = -B(332) ! JVS(140) = Jac_FULL(44,44) JVS(140) = -B(471)-B(473)-B(475) ! JVS(141) = Jac_FULL(44,100) JVS(141) = -B(476) ! JVS(142) = Jac_FULL(44,102) JVS(142) = -B(472)-B(474) ! JVS(143) = Jac_FULL(45,45) JVS(143) = -B(349)-B(546) ! JVS(144) = Jac_FULL(45,96) JVS(144) = 0.7*B(251) ! JVS(145) = Jac_FULL(45,97) JVS(145) = 0.7*B(252) ! JVS(146) = Jac_FULL(45,102) JVS(146) = -B(350) ! JVS(147) = Jac_FULL(46,46) JVS(147) = -B(323)-B(533) ! JVS(148) = Jac_FULL(46,87) JVS(148) = B(239) ! JVS(149) = Jac_FULL(46,96) JVS(149) = B(240) ! JVS(150) = Jac_FULL(46,102) JVS(150) = -B(324) ! JVS(151) = Jac_FULL(47,47) JVS(151) = -B(335)-B(539) ! JVS(152) = Jac_FULL(47,94) JVS(152) = 0.71*B(249) ! JVS(153) = Jac_FULL(47,96) JVS(153) = 0.71*B(250) ! JVS(154) = Jac_FULL(47,102) JVS(154) = -B(336) ! JVS(155) = Jac_FULL(48,48) JVS(155) = -B(351)-B(548) ! JVS(156) = Jac_FULL(48,96) JVS(156) = 0.41*B(245) ! JVS(157) = Jac_FULL(48,101) JVS(157) = 0.41*B(246) ! JVS(158) = Jac_FULL(48,102) JVS(158) = -B(352) ! JVS(159) = Jac_FULL(49,49) JVS(159) = -B(48)-B(49)-B(506)-B(511) ! JVS(160) = Jac_FULL(49,96) JVS(160) = B(46) ! JVS(161) = Jac_FULL(49,99) JVS(161) = B(47) ! JVS(162) = Jac_FULL(49,102) JVS(162) = -B(50) ! JVS(163) = Jac_FULL(50,50) JVS(163) = -B(333)-B(537) ! JVS(164) = Jac_FULL(50,83) JVS(164) = B(243) ! JVS(165) = Jac_FULL(50,96) JVS(165) = B(244) ! JVS(166) = Jac_FULL(50,102) JVS(166) = -B(334) ! JVS(167) = Jac_FULL(51,51) JVS(167) = -B(321)-B(532) ! JVS(168) = Jac_FULL(51,76) JVS(168) = B(195) ! JVS(169) = Jac_FULL(51,96) JVS(169) = B(196) ! JVS(170) = Jac_FULL(51,102) JVS(170) = -B(322) ! JVS(171) = Jac_FULL(52,52) JVS(171) = -B(319)-B(531) ! JVS(172) = Jac_FULL(52,78) JVS(172) = B(193) ! JVS(173) = Jac_FULL(52,96) JVS(173) = B(194) ! JVS(174) = Jac_FULL(52,102) JVS(174) = -B(320) ! JVS(175) = Jac_FULL(53,53) JVS(175) = -B(347)-B(545) ! JVS(176) = Jac_FULL(53,72) JVS(176) = B(185) ! JVS(177) = Jac_FULL(53,96) JVS(177) = B(186) ! JVS(178) = Jac_FULL(53,102) JVS(178) = -B(348) ! JVS(179) = Jac_FULL(54,54) JVS(179) = -B(341)-B(542) ! JVS(180) = Jac_FULL(54,79) JVS(180) = B(179) ! JVS(181) = Jac_FULL(54,96) JVS(181) = B(180) ! JVS(182) = Jac_FULL(54,102) JVS(182) = -B(342) ! JVS(183) = Jac_FULL(55,55) JVS(183) = -B(345)-B(544) ! JVS(184) = Jac_FULL(55,75) JVS(184) = B(183) ! JVS(185) = Jac_FULL(55,96) JVS(185) = B(184) ! JVS(186) = Jac_FULL(55,102) JVS(186) = -B(346) ! JVS(187) = Jac_FULL(56,56) JVS(187) = -B(343)-B(543) ! JVS(188) = Jac_FULL(56,64) JVS(188) = B(187) ! JVS(189) = Jac_FULL(56,69) JVS(189) = B(189) ! JVS(190) = Jac_FULL(56,80) JVS(190) = B(181) ! JVS(191) = Jac_FULL(56,96) JVS(191) = B(182)+B(188)+B(190) ! JVS(192) = Jac_FULL(56,102) JVS(192) = -B(344) ! JVS(193) = Jac_FULL(57,57) JVS(193) = -B(339)-B(541) ! JVS(194) = Jac_FULL(57,70) JVS(194) = B(177) ! JVS(195) = Jac_FULL(57,92) JVS(195) = B(175) ! JVS(196) = Jac_FULL(57,96) JVS(196) = B(176)+B(178) ! JVS(197) = Jac_FULL(57,102) JVS(197) = -B(340) ! JVS(198) = Jac_FULL(58,58) JVS(198) = -B(145)-B(259)-B(261)-B(464) ! JVS(199) = Jac_FULL(58,97) JVS(199) = B(143) ! JVS(200) = Jac_FULL(58,98) JVS(200) = -B(262) ! JVS(201) = Jac_FULL(58,99) JVS(201) = B(144) ! JVS(202) = Jac_FULL(58,102) JVS(202) = -B(260) ! JVS(203) = Jac_FULL(59,59) JVS(203) = -B(275)-B(291)-B(297) ! JVS(204) = Jac_FULL(59,98) JVS(204) = -B(292) ! JVS(205) = Jac_FULL(59,100) JVS(205) = -B(298) ! JVS(206) = Jac_FULL(59,102) JVS(206) = -B(276) ! JVS(207) = Jac_FULL(60,53) JVS(207) = 0.5*B(545) ! JVS(208) = Jac_FULL(60,54) JVS(208) = 0.67*B(542) ! JVS(209) = Jac_FULL(60,59) JVS(209) = 0.05*B(291) ! JVS(210) = Jac_FULL(60,60) JVS(210) = -B(22) ! JVS(211) = Jac_FULL(60,61) JVS(211) = 0.42*B(257) ! JVS(212) = Jac_FULL(60,62) JVS(212) = B(517) ! JVS(213) = Jac_FULL(60,63) JVS(213) = 0.4*B(263)+B(520) ! JVS(214) = Jac_FULL(60,71) JVS(214) = 0.4*B(359) ! JVS(215) = Jac_FULL(60,72) JVS(215) = 0.15*B(221)+0.83*B(382) ! JVS(216) = Jac_FULL(60,79) JVS(216) = 0.61*B(107)+0.33*B(215)+0.65*B(376) ! JVS(217) = Jac_FULL(60,81) JVS(217) = B(515) ! JVS(218) = Jac_FULL(60,82) JVS(218) = B(36)+B(66)+B(502)+B(503) ! JVS(219) = Jac_FULL(60,84) JVS(219) = 0.05*B(68)+B(512)+B(513) ! JVS(220) = Jac_FULL(60,88) JVS(220) = B(269)+B(273)+B(523)+B(524) ! JVS(221) = Jac_FULL(60,91) JVS(221) = 0.2*B(295)+B(529) ! JVS(222) = Jac_FULL(60,93) JVS(222) = 0.05*B(293)+B(525)+B(526) ! JVS(223) = Jac_FULL(60,95) JVS(223) = 0.33*B(216)+0.15*B(222)+B(309) ! JVS(224) = Jac_FULL(60,96) JVS(224) = 0 ! JVS(225) = Jac_FULL(60,98) JVS(225) = 0.42*B(258)+0.05*B(292)+0.05*B(294)+0.2*B(296)+0.4*B(360) ! JVS(226) = Jac_FULL(60,100) JVS(226) = B(67)+2*B(271)+B(274) ! JVS(227) = Jac_FULL(60,101) JVS(227) = 0.65*B(377)+0.83*B(383)+B(436) ! JVS(228) = Jac_FULL(60,102) JVS(228) = -B(23)+B(37)+0.05*B(69)+0.4*B(264)+2*B(267)+B(270)+0.5*B(337) ! JVS(229) = Jac_FULL(60,103) JVS(229) = 0.61*B(108)+B(155) ! JVS(230) = Jac_FULL(61,59) JVS(230) = 0.07*B(291) ! JVS(231) = Jac_FULL(61,61) JVS(231) = -B(255)-B(257)-B(265) ! JVS(232) = Jac_FULL(61,93) JVS(232) = B(525) ! JVS(233) = Jac_FULL(61,98) JVS(233) = -B(258)+0.07*B(292) ! JVS(234) = Jac_FULL(61,100) JVS(234) = -B(266) ! JVS(235) = Jac_FULL(61,102) JVS(235) = -B(256) ! JVS(236) = Jac_FULL(62,42) JVS(236) = 0.5*B(327)+B(535) ! JVS(237) = Jac_FULL(62,62) JVS(237) = -B(159)-B(161)-B(516)-B(517) ! JVS(238) = Jac_FULL(62,66) JVS(238) = B(123)+0.75*B(227)+B(384)+B(424) ! JVS(239) = Jac_FULL(62,85) JVS(239) = 0.32*B(91)+0.16*B(203)+0.32*B(366) ! JVS(240) = Jac_FULL(62,86) JVS(240) = 0.32*B(547) ! JVS(241) = Jac_FULL(62,95) JVS(241) = 0.16*B(204)+0.75*B(228) ! JVS(242) = Jac_FULL(62,96) JVS(242) = 0 ! JVS(243) = Jac_FULL(62,101) JVS(243) = 0.32*B(367)+B(385)+B(425) ! JVS(244) = Jac_FULL(62,102) JVS(244) = -B(160)-B(162)+0.5*B(328) ! JVS(245) = Jac_FULL(62,103) JVS(245) = 0.32*B(92)+B(124) ! JVS(246) = Jac_FULL(63,54) JVS(246) = 0.26*B(542) ! JVS(247) = Jac_FULL(63,55) JVS(247) = 0.7*B(544) ! JVS(248) = Jac_FULL(63,63) JVS(248) = -B(263)-B(520) ! JVS(249) = Jac_FULL(63,64) JVS(249) = 0.6*B(119) ! JVS(250) = Jac_FULL(63,71) JVS(250) = 0.28*B(359) ! JVS(251) = Jac_FULL(63,75) JVS(251) = 0.72*B(111)+0.36*B(219)+0.72*B(380) ! JVS(252) = Jac_FULL(63,79) JVS(252) = 0.24*B(107)+0.13*B(215)+0.26*B(376) ! JVS(253) = Jac_FULL(63,80) JVS(253) = 0.95*B(109)+0.5*B(217)+B(378) ! JVS(254) = Jac_FULL(63,95) JVS(254) = 0.13*B(216)+0.5*B(218)+0.36*B(220) ! JVS(255) = Jac_FULL(63,96) JVS(255) = 0 ! JVS(256) = Jac_FULL(63,98) JVS(256) = 0.28*B(360) ! JVS(257) = Jac_FULL(63,101) JVS(257) = 0.26*B(377)+B(379)+0.72*B(381) ! JVS(258) = Jac_FULL(63,102) JVS(258) = -B(264) ! JVS(259) = Jac_FULL(63,103) JVS(259) = 0.24*B(108)+0.95*B(110)+0.72*B(112)+0.6*B(120) ! JVS(260) = Jac_FULL(64,64) JVS(260) = -B(119)-B(187)-B(223)-B(388)-B(416) ! JVS(261) = Jac_FULL(64,95) JVS(261) = -B(224) ! JVS(262) = Jac_FULL(64,96) JVS(262) = -B(188) ! JVS(263) = Jac_FULL(64,101) JVS(263) = -B(389)-B(417) ! JVS(264) = Jac_FULL(64,103) JVS(264) = -B(120) ! JVS(265) = Jac_FULL(65,26) JVS(265) = B(83) ! JVS(266) = Jac_FULL(65,41) JVS(266) = 0.5*B(325) ! JVS(267) = Jac_FULL(65,65) JVS(267) = -B(85)-B(163)-B(241)-B(283)-B(287) ! JVS(268) = Jac_FULL(65,85) JVS(268) = 0.05*B(91)+0.03*B(203)+0.05*B(366) ! JVS(269) = Jac_FULL(65,86) JVS(269) = 0.05*B(547) ! JVS(270) = Jac_FULL(65,95) JVS(270) = -B(164)+0.03*B(204) ! JVS(271) = Jac_FULL(65,96) JVS(271) = -B(242) ! JVS(272) = Jac_FULL(65,101) JVS(272) = -B(284)-B(288)+0.05*B(367) ! JVS(273) = Jac_FULL(65,102) JVS(273) = B(84)+0.5*B(326) ! JVS(274) = Jac_FULL(65,103) JVS(274) = -B(86)+0.05*B(92) ! JVS(275) = Jac_FULL(66,26) JVS(275) = B(81) ! JVS(276) = Jac_FULL(66,42) JVS(276) = 0.5*B(327) ! JVS(277) = Jac_FULL(66,66) JVS(277) = -B(123)-B(191)-B(227)-B(384)-B(424) ! JVS(278) = Jac_FULL(66,85) JVS(278) = 0.18*B(91)+0.09*B(203)+0.18*B(366) ! JVS(279) = Jac_FULL(66,86) JVS(279) = 0.18*B(547) ! JVS(280) = Jac_FULL(66,95) JVS(280) = 0.09*B(204)-B(228) ! JVS(281) = Jac_FULL(66,96) JVS(281) = -B(192) ! JVS(282) = Jac_FULL(66,101) JVS(282) = 0.18*B(367)-B(385)-B(425) ! JVS(283) = Jac_FULL(66,102) JVS(283) = B(82)+0.5*B(328) ! JVS(284) = Jac_FULL(66,103) JVS(284) = 0.18*B(92)-B(124) ! JVS(285) = Jac_FULL(67,35) JVS(285) = B(129) ! JVS(286) = Jac_FULL(67,36) JVS(286) = B(353) ! JVS(287) = Jac_FULL(67,38) JVS(287) = 2*B(470) ! JVS(288) = Jac_FULL(67,44) JVS(288) = B(475) ! JVS(289) = Jac_FULL(67,64) JVS(289) = 0.1*B(119) ! JVS(290) = Jac_FULL(67,67) JVS(290) = -B(40)-B(460)-B(504) ! JVS(291) = Jac_FULL(67,70) JVS(291) = B(105) ! JVS(292) = Jac_FULL(67,72) JVS(292) = B(117) ! JVS(293) = Jac_FULL(67,75) JVS(293) = B(113) ! JVS(294) = Jac_FULL(67,78) JVS(294) = 0.85*B(125)+0.425*B(229)+0.85*B(392) ! JVS(295) = Jac_FULL(67,79) JVS(295) = 0.08*B(107) ! JVS(296) = Jac_FULL(67,80) JVS(296) = 0.05*B(109) ! JVS(297) = Jac_FULL(67,81) JVS(297) = B(157) ! JVS(298) = Jac_FULL(67,82) JVS(298) = B(66) ! JVS(299) = Jac_FULL(67,84) JVS(299) = B(70) ! JVS(300) = Jac_FULL(67,88) JVS(300) = B(273) ! JVS(301) = Jac_FULL(67,89) JVS(301) = B(201) ! JVS(302) = Jac_FULL(67,91) JVS(302) = B(301) ! JVS(303) = Jac_FULL(67,95) JVS(303) = 0.425*B(230) ! JVS(304) = Jac_FULL(67,96) JVS(304) = 0 ! JVS(305) = Jac_FULL(67,99) JVS(305) = B(38)+0.5*B(468) ! JVS(306) = Jac_FULL(67,100) JVS(306) = B(67)+B(71)+B(130)+B(158)+B(202)+B(271)+B(274)+B(302)+B(354)+B(469)+B(476) ! JVS(307) = Jac_FULL(67,101) JVS(307) = 0.85*B(393) ! JVS(308) = Jac_FULL(67,102) JVS(308) = B(39)-B(41) ! JVS(309) = Jac_FULL(67,103) JVS(309) = B(106)+0.08*B(108)+0.05*B(110)+B(114)+B(118)+0.1*B(120)+0.85*B(126) ! JVS(310) = Jac_FULL(68,68) JVS(310) = -B(95)-B(169)-B(205)-B(386)-B(412) ! JVS(311) = Jac_FULL(68,86) JVS(311) = B(131) ! JVS(312) = Jac_FULL(68,95) JVS(312) = -B(206) ! JVS(313) = Jac_FULL(68,96) JVS(313) = -B(170) ! JVS(314) = Jac_FULL(68,101) JVS(314) = -B(387)-B(413) ! JVS(315) = Jac_FULL(68,102) JVS(315) = B(132) ! JVS(316) = Jac_FULL(68,103) JVS(316) = -B(96) ! JVS(317) = Jac_FULL(69,69) JVS(317) = -B(121)-B(189)-B(225)-B(390)-B(418) ! JVS(318) = Jac_FULL(69,91) JVS(318) = B(299) ! JVS(319) = Jac_FULL(69,95) JVS(319) = -B(226) ! JVS(320) = Jac_FULL(69,96) JVS(320) = -B(190) ! JVS(321) = Jac_FULL(69,100) JVS(321) = B(300) ! JVS(322) = Jac_FULL(69,101) JVS(322) = -B(391)-B(419) ! JVS(323) = Jac_FULL(69,103) JVS(323) = -B(122) ! JVS(324) = Jac_FULL(70,70) JVS(324) = -B(103)-B(105)-B(177)-B(213)-B(374)-B(404) ! JVS(325) = Jac_FULL(70,92) JVS(325) = 0.07*B(211)+0.136*B(372) ! JVS(326) = Jac_FULL(70,95) JVS(326) = 0.07*B(212)-B(214) ! JVS(327) = Jac_FULL(70,96) JVS(327) = -B(178) ! JVS(328) = Jac_FULL(70,101) JVS(328) = 0.136*B(373)-B(375)-B(405) ! JVS(329) = Jac_FULL(70,103) JVS(329) = -B(104)-B(106) ! JVS(330) = Jac_FULL(71,57) JVS(330) = 0.509*B(339)+0.373*B(541) ! JVS(331) = Jac_FULL(71,70) JVS(331) = B(103)+0.5*B(213)+B(374) ! JVS(332) = Jac_FULL(71,71) JVS(332) = -B(357)-B(359) ! JVS(333) = Jac_FULL(71,92) JVS(333) = 0.34*B(101)+0.06*B(211)+0.127*B(372) ! JVS(334) = Jac_FULL(71,95) JVS(334) = 0.06*B(212)+0.5*B(214) ! JVS(335) = Jac_FULL(71,96) JVS(335) = 0 ! JVS(336) = Jac_FULL(71,98) JVS(336) = -B(360) ! JVS(337) = Jac_FULL(71,101) JVS(337) = 0.127*B(373)+B(375) ! JVS(338) = Jac_FULL(71,102) JVS(338) = 0.509*B(340)-B(358) ! JVS(339) = Jac_FULL(71,103) JVS(339) = 0.34*B(102)+B(104) ! JVS(340) = Jac_FULL(72,53) JVS(340) = B(347) ! JVS(341) = Jac_FULL(72,72) JVS(341) = -B(115)-B(117)-B(185)-B(221)-B(382)-B(410) ! JVS(342) = Jac_FULL(72,91) JVS(342) = 0.43*B(279) ! JVS(343) = Jac_FULL(72,95) JVS(343) = -B(222) ! JVS(344) = Jac_FULL(72,96) JVS(344) = -B(186) ! JVS(345) = Jac_FULL(72,101) JVS(345) = -B(383)-B(411) ! JVS(346) = Jac_FULL(72,102) JVS(346) = 0.43*B(280)+B(348) ! JVS(347) = Jac_FULL(72,103) JVS(347) = -B(116)-B(118) ! JVS(348) = Jac_FULL(73,73) JVS(348) = -B(99)-B(173)-B(209)-B(370)-B(400) ! JVS(349) = Jac_FULL(73,89) JVS(349) = B(197)+B(201) ! JVS(350) = Jac_FULL(73,95) JVS(350) = -B(210) ! JVS(351) = Jac_FULL(73,96) JVS(351) = -B(174) ! JVS(352) = Jac_FULL(73,100) JVS(352) = B(202) ! JVS(353) = Jac_FULL(73,101) JVS(353) = -B(371)-B(401) ! JVS(354) = Jac_FULL(73,102) JVS(354) = B(198) ! JVS(355) = Jac_FULL(73,103) JVS(355) = -B(100) ! JVS(356) = Jac_FULL(74,62) JVS(356) = B(159)+B(161) ! JVS(357) = Jac_FULL(74,66) JVS(357) = 0 ! JVS(358) = Jac_FULL(74,74) JVS(358) = -B(97)-B(171)-B(207)-B(368)-B(398) ! JVS(359) = Jac_FULL(74,85) JVS(359) = 0 ! JVS(360) = Jac_FULL(74,86) JVS(360) = 0 ! JVS(361) = Jac_FULL(74,95) JVS(361) = -B(208) ! JVS(362) = Jac_FULL(74,96) JVS(362) = -B(172) ! JVS(363) = Jac_FULL(74,101) JVS(363) = -B(369)-B(399) ! JVS(364) = Jac_FULL(74,102) JVS(364) = B(160)+B(162) ! JVS(365) = Jac_FULL(74,103) JVS(365) = -B(98) ! JVS(366) = Jac_FULL(75,55) JVS(366) = 0.5*B(345) ! JVS(367) = Jac_FULL(75,75) JVS(367) = -B(111)-B(113)-B(183)-B(219)-B(380)-B(408) ! JVS(368) = Jac_FULL(75,93) JVS(368) = B(277) ! JVS(369) = Jac_FULL(75,95) JVS(369) = -B(220) ! JVS(370) = Jac_FULL(75,96) JVS(370) = -B(184) ! JVS(371) = Jac_FULL(75,101) JVS(371) = -B(381)-B(409) ! JVS(372) = Jac_FULL(75,102) JVS(372) = B(278)+0.5*B(346) ! JVS(373) = Jac_FULL(75,103) JVS(373) = -B(112)-B(114) ! JVS(374) = Jac_FULL(76,51) JVS(374) = B(321) ! JVS(375) = Jac_FULL(76,61) JVS(375) = B(265) ! JVS(376) = Jac_FULL(76,76) JVS(376) = -B(127)-B(195)-B(231)-B(394)-B(422) ! JVS(377) = Jac_FULL(76,93) JVS(377) = 0 ! JVS(378) = Jac_FULL(76,95) JVS(378) = -B(232) ! JVS(379) = Jac_FULL(76,96) JVS(379) = -B(196) ! JVS(380) = Jac_FULL(76,98) JVS(380) = 0 ! JVS(381) = Jac_FULL(76,100) JVS(381) = B(266) ! JVS(382) = Jac_FULL(76,101) JVS(382) = -B(395)-B(423) ! JVS(383) = Jac_FULL(76,102) JVS(383) = B(322) ! JVS(384) = Jac_FULL(76,103) JVS(384) = -B(128) ! JVS(385) = Jac_FULL(77,53) JVS(385) = B(545) ! JVS(386) = Jac_FULL(77,54) JVS(386) = 0.36*B(542) ! JVS(387) = Jac_FULL(77,58) JVS(387) = 0.59*B(259) ! JVS(388) = Jac_FULL(77,71) JVS(388) = 0.2*B(359) ! JVS(389) = Jac_FULL(77,72) JVS(389) = B(115)+B(221)+0.83*B(382) ! JVS(390) = Jac_FULL(77,74) JVS(390) = 0.2*B(207) ! JVS(391) = Jac_FULL(77,77) JVS(391) = -B(281)-B(530) ! JVS(392) = Jac_FULL(77,79) JVS(392) = 0.33*B(107)+0.18*B(215)+0.36*B(376) ! JVS(393) = Jac_FULL(77,80) JVS(393) = 0.95*B(109)+0.5*B(217)+B(378) ! JVS(394) = Jac_FULL(77,83) JVS(394) = 0.16*B(165)+0.65*B(289) ! JVS(395) = Jac_FULL(77,85) JVS(395) = 0 ! JVS(396) = Jac_FULL(77,86) JVS(396) = 0 ! JVS(397) = Jac_FULL(77,91) JVS(397) = 0 ! JVS(398) = Jac_FULL(77,92) JVS(398) = 0 ! JVS(399) = Jac_FULL(77,95) JVS(399) = 0.16*B(166)+0.2*B(208)+0.18*B(216)+0.5*B(218)+B(222) ! JVS(400) = Jac_FULL(77,96) JVS(400) = 0 ! JVS(401) = Jac_FULL(77,97) JVS(401) = 0 ! JVS(402) = Jac_FULL(77,98) JVS(402) = 0.2*B(360) ! JVS(403) = Jac_FULL(77,99) JVS(403) = 0 ! JVS(404) = Jac_FULL(77,101) JVS(404) = 0.65*B(290)+0.36*B(377)+B(379)+0.83*B(383) ! JVS(405) = Jac_FULL(77,102) JVS(405) = 0.59*B(260)-B(282) ! JVS(406) = Jac_FULL(77,103) JVS(406) = 0.33*B(108)+0.95*B(110)+B(116) ! JVS(407) = Jac_FULL(78,52) JVS(407) = B(319) ! JVS(408) = Jac_FULL(78,59) JVS(408) = B(297) ! JVS(409) = Jac_FULL(78,78) JVS(409) = -B(125)-B(193)-B(229)-B(392)-B(420) ! JVS(410) = Jac_FULL(78,95) JVS(410) = -B(230) ! JVS(411) = Jac_FULL(78,96) JVS(411) = -B(194) ! JVS(412) = Jac_FULL(78,98) JVS(412) = 0 ! JVS(413) = Jac_FULL(78,100) JVS(413) = B(298) ! JVS(414) = Jac_FULL(78,101) JVS(414) = -B(393)-B(421) ! JVS(415) = Jac_FULL(78,102) JVS(415) = B(320) ! JVS(416) = Jac_FULL(78,103) JVS(416) = -B(126) ! JVS(417) = Jac_FULL(79,54) JVS(417) = B(341) ! JVS(418) = Jac_FULL(79,71) JVS(418) = 0.44*B(357) ! JVS(419) = Jac_FULL(79,79) JVS(419) = -B(107)-B(179)-B(215)-B(376)-B(406) ! JVS(420) = Jac_FULL(79,92) JVS(420) = 0 ! JVS(421) = Jac_FULL(79,95) JVS(421) = -B(216) ! JVS(422) = Jac_FULL(79,96) JVS(422) = -B(180) ! JVS(423) = Jac_FULL(79,98) JVS(423) = 0 ! JVS(424) = Jac_FULL(79,101) JVS(424) = -B(377)-B(407) ! JVS(425) = Jac_FULL(79,102) JVS(425) = B(342)+0.44*B(358) ! JVS(426) = Jac_FULL(79,103) JVS(426) = -B(108) ! JVS(427) = Jac_FULL(80,56) JVS(427) = 0.5*B(343) ! JVS(428) = Jac_FULL(80,64) JVS(428) = 0 ! JVS(429) = Jac_FULL(80,69) JVS(429) = 0 ! JVS(430) = Jac_FULL(80,80) JVS(430) = -B(109)-B(181)-B(217)-B(378)-B(414) ! JVS(431) = Jac_FULL(80,91) JVS(431) = 0 ! JVS(432) = Jac_FULL(80,95) JVS(432) = -B(218) ! JVS(433) = Jac_FULL(80,96) JVS(433) = -B(182) ! JVS(434) = Jac_FULL(80,100) JVS(434) = 0 ! JVS(435) = Jac_FULL(80,101) JVS(435) = -B(379)-B(415) ! JVS(436) = Jac_FULL(80,102) JVS(436) = 0.5*B(344) ! JVS(437) = Jac_FULL(80,103) JVS(437) = -B(110) ! JVS(438) = Jac_FULL(81,40) JVS(438) = 0.5*B(329)+B(536) ! JVS(439) = Jac_FULL(81,41) JVS(439) = 0.5*B(325)+B(534) ! JVS(440) = Jac_FULL(81,51) JVS(440) = B(532) ! JVS(441) = Jac_FULL(81,52) JVS(441) = B(531) ! JVS(442) = Jac_FULL(81,55) JVS(442) = 0.5*B(345) ! JVS(443) = Jac_FULL(81,56) JVS(443) = 0.5*B(343)+B(543) ! JVS(444) = Jac_FULL(81,64) JVS(444) = 0.25*B(223)+B(416) ! JVS(445) = Jac_FULL(81,65) JVS(445) = B(85)+0.75*B(163)+B(283)+B(287) ! JVS(446) = Jac_FULL(81,68) JVS(446) = 0.57*B(95)+0.54*B(205)+0.57*B(386)+B(412) ! JVS(447) = Jac_FULL(81,69) JVS(447) = 0.25*B(225)+B(418) ! JVS(448) = Jac_FULL(81,75) JVS(448) = 0 ! JVS(449) = Jac_FULL(81,76) JVS(449) = 0.25*B(231)+B(422) ! JVS(450) = Jac_FULL(81,78) JVS(450) = 0.25*B(229)+B(420) ! JVS(451) = Jac_FULL(81,80) JVS(451) = 0.25*B(217)+B(414) ! JVS(452) = Jac_FULL(81,81) JVS(452) = -B(135)-B(157)-B(515) ! JVS(453) = Jac_FULL(81,83) JVS(453) = 0.09*B(165)+0.35*B(289) ! JVS(454) = Jac_FULL(81,85) JVS(454) = 0.13*B(91)+0.07*B(203)+0.13*B(366) ! JVS(455) = Jac_FULL(81,86) JVS(455) = 0.13*B(547) ! JVS(456) = Jac_FULL(81,91) JVS(456) = 0 ! JVS(457) = Jac_FULL(81,93) JVS(457) = 0 ! JVS(458) = Jac_FULL(81,95) JVS(458) = 0.75*B(164)+0.09*B(166)+0.07*B(204)+0.54*B(206)+0.25*B(218)+0.25*B(224)+0.25*B(226)+0.25*B(230)+0.25*B(232) ! JVS(459) = Jac_FULL(81,96) JVS(459) = 0 ! JVS(460) = Jac_FULL(81,98) JVS(460) = 0 ! JVS(461) = Jac_FULL(81,100) JVS(461) = -B(158) ! JVS(462) = Jac_FULL(81,101) JVS(462) = B(284)+B(288)+0.35*B(290)+0.13*B(367)+0.57*B(387)+B(413)+B(415)+B(417)+B(419)+B(421)+B(423) ! JVS(463) = Jac_FULL(81,102) JVS(463) = -B(136)+B(235)+0.5*B(326)+0.5*B(330)+0.5*B(344)+0.5*B(346) ! JVS(464) = Jac_FULL(81,103) JVS(464) = B(86)+0.13*B(92)+0.57*B(96) ! JVS(465) = Jac_FULL(82,1) JVS(465) = 3*B(485) ! JVS(466) = Jac_FULL(82,39) JVS(466) = B(34)+B(501) ! JVS(467) = Jac_FULL(82,44) JVS(467) = B(471)+B(475) ! JVS(468) = Jac_FULL(82,45) JVS(468) = B(546) ! JVS(469) = Jac_FULL(82,47) JVS(469) = B(539) ! JVS(470) = Jac_FULL(82,48) JVS(470) = 0.5*B(351) ! JVS(471) = Jac_FULL(82,50) JVS(471) = B(537) ! JVS(472) = Jac_FULL(82,53) JVS(472) = 0.5*B(545) ! JVS(473) = Jac_FULL(82,55) JVS(473) = 0.3*B(544) ! JVS(474) = Jac_FULL(82,57) JVS(474) = 0.627*B(541) ! JVS(475) = Jac_FULL(82,58) JVS(475) = 2.23*B(259)+0.6*B(261) ! JVS(476) = Jac_FULL(82,59) JVS(476) = 0.9*B(291) ! JVS(477) = Jac_FULL(82,61) JVS(477) = 0.535*B(257) ! JVS(478) = Jac_FULL(82,63) JVS(478) = B(520) ! JVS(479) = Jac_FULL(82,64) JVS(479) = 0.3*B(119)+1.25*B(223)+B(388) ! JVS(480) = Jac_FULL(82,65) JVS(480) = 0.75*B(163) ! JVS(481) = Jac_FULL(82,66) JVS(481) = 0.75*B(227) ! JVS(482) = Jac_FULL(82,68) JVS(482) = 0.39*B(95)+0.95*B(205)+0.39*B(386) ! JVS(483) = Jac_FULL(82,69) JVS(483) = B(121)+1.25*B(225)+B(390) ! JVS(484) = Jac_FULL(82,70) JVS(484) = 0.75*B(103)+1.13*B(213)+0.75*B(374) ! JVS(485) = Jac_FULL(82,71) JVS(485) = 0.12*B(359) ! JVS(486) = Jac_FULL(82,72) JVS(486) = B(115)+0.85*B(221)+0.17*B(382) ! JVS(487) = Jac_FULL(82,73) JVS(487) = 0.75*B(209) ! JVS(488) = Jac_FULL(82,74) JVS(488) = 0.96*B(97)+0.5*B(207)+0.2*B(368) ! JVS(489) = Jac_FULL(82,75) JVS(489) = 0.28*B(111)+0.89*B(219)+0.28*B(380) ! JVS(490) = Jac_FULL(82,76) JVS(490) = B(127)+1.25*B(231)+B(394) ! JVS(491) = Jac_FULL(82,77) JVS(491) = B(530) ! JVS(492) = Jac_FULL(82,78) JVS(492) = 0.15*B(125)+0.83*B(229)+0.15*B(392) ! JVS(493) = Jac_FULL(82,79) JVS(493) = 0.35*B(107)+0.95*B(215)+0.4*B(376) ! JVS(494) = Jac_FULL(82,80) JVS(494) = 0.75*B(217) ! JVS(495) = Jac_FULL(82,82) JVS(495) = -B(36)-B(66)-B(461)-B(502)-B(503) ! JVS(496) = Jac_FULL(82,83) JVS(496) = B(87)+1.25*B(165)+B(285) ! JVS(497) = Jac_FULL(82,84) JVS(497) = 0.05*B(68) ! JVS(498) = Jac_FULL(82,85) JVS(498) = 0.75*B(203) ! JVS(499) = Jac_FULL(82,86) JVS(499) = 0 ! JVS(500) = Jac_FULL(82,87) JVS(500) = 0.75*B(199) ! JVS(501) = Jac_FULL(82,89) JVS(501) = 0 ! JVS(502) = Jac_FULL(82,90) JVS(502) = B(303)+B(311) ! JVS(503) = Jac_FULL(82,91) JVS(503) = 0.7*B(295)+B(529) ! JVS(504) = Jac_FULL(82,92) JVS(504) = 0.56*B(101)+1.1*B(211)+0.69*B(372) ! JVS(505) = Jac_FULL(82,93) JVS(505) = 0.8*B(293)+B(526) ! JVS(506) = Jac_FULL(82,94) JVS(506) = B(151)+0.29*B(249)+2*B(305)+B(313)+B(432) ! JVS(507) = Jac_FULL(82,95) JVS(507) = B(9)+B(26)+B(30)+2*B(31)+0.75*B(164)+1.25*B(166)+0.75*B(200)+0.75*B(204)+0.95*B(206)+0.5*B(208)+0.75*B(210)& &+1.1*B(212)+1.13*B(214)+0.95*B(216)+0.75*B(218)+0.89*B(220)+0.85*B(222)+1.25*B(224)+1.25*B(226)+0.75*B(228)& &+0.83*B(230)+1.25*B(232)+B(304)+2*B(306)+2*B(307)+B(309)+B(312)+B(314)+B(315)+B(317)+B(362)+B(364) ! JVS(508) = Jac_FULL(82,96) JVS(508) = 0.29*B(250) ! JVS(509) = Jac_FULL(82,97) JVS(509) = B(153)+2*B(308)+B(316)+B(434) ! JVS(510) = Jac_FULL(82,98) JVS(510) = B(10)+0.535*B(258)+0.6*B(262)+0.9*B(292)+0.8*B(294)+0.7*B(296)+0.12*B(360) ! JVS(511) = Jac_FULL(82,99) JVS(511) = 0 ! JVS(512) = Jac_FULL(82,100) JVS(512) = -B(67)+B(476) ! JVS(513) = Jac_FULL(82,101) JVS(513) = B(286)+B(363)+B(365)+0.2*B(369)+0.69*B(373)+0.75*B(375)+0.4*B(377)+0.28*B(381)+0.17*B(383)+0.39*B(387)& &+B(389)+B(391)+0.15*B(393)+B(395)+B(433)+B(435) ! JVS(514) = Jac_FULL(82,102) JVS(514) = B(35)-B(37)+B(62)+0.05*B(69)+2.23*B(260)+0.5*B(352)+B(355)+B(472) ! JVS(515) = Jac_FULL(82,103) JVS(515) = B(27)+B(88)+0.39*B(96)+0.96*B(98)+0.56*B(102)+0.75*B(104)+0.35*B(108)+0.28*B(112)+B(116)+0.3*B(120)+B(122)& &+0.15*B(126)+B(128)+B(152)+B(154) ! JVS(516) = Jac_FULL(83,50) JVS(516) = B(333) ! JVS(517) = Jac_FULL(83,61) JVS(517) = B(255) ! JVS(518) = Jac_FULL(83,83) JVS(518) = -B(87)-B(165)-B(243)-B(285)-B(289) ! JVS(519) = Jac_FULL(83,93) JVS(519) = 0 ! JVS(520) = Jac_FULL(83,95) JVS(520) = -B(166) ! JVS(521) = Jac_FULL(83,96) JVS(521) = -B(244) ! JVS(522) = Jac_FULL(83,98) JVS(522) = 0 ! JVS(523) = Jac_FULL(83,100) JVS(523) = 0 ! JVS(524) = Jac_FULL(83,101) JVS(524) = -B(286)-B(290) ! JVS(525) = Jac_FULL(83,102) JVS(525) = B(256)+B(334) ! JVS(526) = Jac_FULL(83,103) JVS(526) = -B(88) ! JVS(527) = Jac_FULL(84,43) JVS(527) = 0.5*B(331)+B(538) ! JVS(528) = Jac_FULL(84,46) JVS(528) = 0.5*B(323)+B(533) ! JVS(529) = Jac_FULL(84,50) JVS(529) = B(537) ! JVS(530) = Jac_FULL(84,61) JVS(530) = 0.5*B(257) ! JVS(531) = Jac_FULL(84,68) JVS(531) = 0.75*B(95)+0.38*B(205)+0.75*B(386) ! JVS(532) = Jac_FULL(84,73) JVS(532) = 0.93*B(99)+0.5*B(209)+B(370) ! JVS(533) = Jac_FULL(84,76) JVS(533) = B(127)+0.5*B(231)+B(394) ! JVS(534) = Jac_FULL(84,83) JVS(534) = B(87)+0.5*B(165)+B(285) ! JVS(535) = Jac_FULL(84,84) JVS(535) = -B(68)-B(70)-B(512)-B(513) ! JVS(536) = Jac_FULL(84,85) JVS(536) = 0.32*B(91)+0.16*B(203)+0.32*B(366) ! JVS(537) = Jac_FULL(84,86) JVS(537) = 0.32*B(547) ! JVS(538) = Jac_FULL(84,87) JVS(538) = B(79)+0.75*B(199)+2*B(237)+B(238)+B(426)+B(428) ! JVS(539) = Jac_FULL(84,88) JVS(539) = B(524) ! JVS(540) = Jac_FULL(84,89) JVS(540) = 0 ! JVS(541) = Jac_FULL(84,90) JVS(541) = 0 ! JVS(542) = Jac_FULL(84,93) JVS(542) = 0.04*B(293) ! JVS(543) = Jac_FULL(84,95) JVS(543) = 0.5*B(166)+0.75*B(200)+0.16*B(204)+0.38*B(206)+0.5*B(210)+0.5*B(232) ! JVS(544) = Jac_FULL(84,96) JVS(544) = 0 ! JVS(545) = Jac_FULL(84,98) JVS(545) = 0.5*B(258)+0.04*B(294) ! JVS(546) = Jac_FULL(84,100) JVS(546) = -B(71) ! JVS(547) = Jac_FULL(84,101) JVS(547) = B(286)+0.32*B(367)+B(371)+0.75*B(387)+B(395)+B(427)+B(429) ! JVS(548) = Jac_FULL(84,102) JVS(548) = -B(69)+B(233)+0.5*B(324)+0.5*B(332) ! JVS(549) = Jac_FULL(84,103) JVS(549) = B(80)+B(88)+0.32*B(92)+0.75*B(96)+0.93*B(100)+B(128) ! JVS(550) = Jac_FULL(85,35) JVS(550) = B(89)+B(129) ! JVS(551) = Jac_FULL(85,40) JVS(551) = 0.5*B(329) ! JVS(552) = Jac_FULL(85,68) JVS(552) = 0.3*B(95)+0.15*B(205)+0.3*B(386) ! JVS(553) = Jac_FULL(85,85) JVS(553) = -B(91)-B(93)-B(167)-B(203)-B(366)-B(396) ! JVS(554) = Jac_FULL(85,86) JVS(554) = 0 ! JVS(555) = Jac_FULL(85,95) JVS(555) = -B(204)+0.15*B(206) ! JVS(556) = Jac_FULL(85,96) JVS(556) = -B(168) ! JVS(557) = Jac_FULL(85,100) JVS(557) = B(130) ! JVS(558) = Jac_FULL(85,101) JVS(558) = -B(367)+0.3*B(387)-B(397) ! JVS(559) = Jac_FULL(85,102) JVS(559) = B(90)+0.5*B(330) ! JVS(560) = Jac_FULL(85,103) JVS(560) = -B(92)-B(94)+0.3*B(96) ! JVS(561) = Jac_FULL(86,68) JVS(561) = B(169) ! JVS(562) = Jac_FULL(86,73) JVS(562) = 0.07*B(99) ! JVS(563) = Jac_FULL(86,74) JVS(563) = 0.04*B(97) ! JVS(564) = Jac_FULL(86,85) JVS(564) = B(93) ! JVS(565) = Jac_FULL(86,86) JVS(565) = -B(131)-B(466)-B(547) ! JVS(566) = Jac_FULL(86,89) JVS(566) = 0 ! JVS(567) = Jac_FULL(86,95) JVS(567) = 0 ! JVS(568) = Jac_FULL(86,96) JVS(568) = B(170) ! JVS(569) = Jac_FULL(86,100) JVS(569) = 0 ! JVS(570) = Jac_FULL(86,101) JVS(570) = 0 ! JVS(571) = Jac_FULL(86,102) JVS(571) = -B(132) ! JVS(572) = Jac_FULL(86,103) JVS(572) = B(94)+0.04*B(98)+0.07*B(100) ! JVS(573) = Jac_FULL(87,36) JVS(573) = B(77)+B(353) ! JVS(574) = Jac_FULL(87,46) JVS(574) = 0.5*B(323) ! JVS(575) = Jac_FULL(87,81) JVS(575) = B(515) ! JVS(576) = Jac_FULL(87,83) JVS(576) = 0 ! JVS(577) = Jac_FULL(87,85) JVS(577) = 0.32*B(91)+0.16*B(203)+0.32*B(366) ! JVS(578) = Jac_FULL(87,86) JVS(578) = 0.32*B(547) ! JVS(579) = Jac_FULL(87,87) JVS(579) = -B(79)-B(199)-2*B(237)-2*B(238)-B(239)-B(426)-B(428) ! JVS(580) = Jac_FULL(87,89) JVS(580) = 0.85*B(518) ! JVS(581) = Jac_FULL(87,90) JVS(581) = B(149)+B(303)+B(430) ! JVS(582) = Jac_FULL(87,91) JVS(582) = 0 ! JVS(583) = Jac_FULL(87,93) JVS(583) = 0 ! JVS(584) = Jac_FULL(87,95) JVS(584) = -B(200)+0.16*B(204)+B(304) ! JVS(585) = Jac_FULL(87,96) JVS(585) = -B(240) ! JVS(586) = Jac_FULL(87,98) JVS(586) = 0 ! JVS(587) = Jac_FULL(87,100) JVS(587) = B(354) ! JVS(588) = Jac_FULL(87,101) JVS(588) = 0.32*B(367)-B(427)-B(429)+B(431) ! JVS(589) = Jac_FULL(87,102) JVS(589) = B(78)+0.5*B(324) ! JVS(590) = Jac_FULL(87,103) JVS(590) = -B(80)+0.32*B(92)+B(150) ! JVS(591) = Jac_FULL(88,54) JVS(591) = 0.58*B(542) ! JVS(592) = Jac_FULL(88,55) JVS(592) = 0.3*B(544) ! JVS(593) = Jac_FULL(88,64) JVS(593) = 0.3*B(119)+0.25*B(223)+0.5*B(388) ! JVS(594) = Jac_FULL(88,69) JVS(594) = B(121)+0.5*B(225)+B(390) ! JVS(595) = Jac_FULL(88,71) JVS(595) = 0.6*B(359) ! JVS(596) = Jac_FULL(88,72) JVS(596) = 0.17*B(382) ! JVS(597) = Jac_FULL(88,73) JVS(597) = B(173) ! JVS(598) = Jac_FULL(88,74) JVS(598) = 0.5*B(207)+0.8*B(368) ! JVS(599) = Jac_FULL(88,75) JVS(599) = 0.28*B(111)+0.14*B(219)+0.28*B(380) ! JVS(600) = Jac_FULL(88,77) JVS(600) = B(281) ! JVS(601) = Jac_FULL(88,79) JVS(601) = 0.53*B(107)+0.29*B(215)+0.58*B(376) ! JVS(602) = Jac_FULL(88,80) JVS(602) = 0 ! JVS(603) = Jac_FULL(88,83) JVS(603) = 0 ! JVS(604) = Jac_FULL(88,85) JVS(604) = 0 ! JVS(605) = Jac_FULL(88,86) JVS(605) = 0 ! JVS(606) = Jac_FULL(88,88) JVS(606) = -B(269)-B(273)-B(523)-B(524) ! JVS(607) = Jac_FULL(88,89) JVS(607) = 0 ! JVS(608) = Jac_FULL(88,91) JVS(608) = 0.8*B(295) ! JVS(609) = Jac_FULL(88,92) JVS(609) = 0 ! JVS(610) = Jac_FULL(88,93) JVS(610) = 0.82*B(293) ! JVS(611) = Jac_FULL(88,95) JVS(611) = 0.5*B(208)+0.29*B(216)+0.14*B(220)+0.25*B(224)+0.5*B(226) ! JVS(612) = Jac_FULL(88,96) JVS(612) = B(174) ! JVS(613) = Jac_FULL(88,97) JVS(613) = 0 ! JVS(614) = Jac_FULL(88,98) JVS(614) = 0.82*B(294)+0.8*B(296)+0.6*B(360) ! JVS(615) = Jac_FULL(88,99) JVS(615) = 0 ! JVS(616) = Jac_FULL(88,100) JVS(616) = -B(274) ! JVS(617) = Jac_FULL(88,101) JVS(617) = 0.8*B(369)+0.58*B(377)+0.28*B(381)+0.17*B(383)+0.5*B(389)+B(391) ! JVS(618) = Jac_FULL(88,102) JVS(618) = -B(270)+B(282) ! JVS(619) = Jac_FULL(88,103) JVS(619) = 0.53*B(108)+0.28*B(112)+0.3*B(120)+B(122) ! JVS(620) = Jac_FULL(89,70) JVS(620) = 0.25*B(213)+B(404) ! JVS(621) = Jac_FULL(89,72) JVS(621) = B(410) ! JVS(622) = Jac_FULL(89,73) JVS(622) = 0.25*B(209)+B(400) ! JVS(623) = Jac_FULL(89,74) JVS(623) = B(398) ! JVS(624) = Jac_FULL(89,75) JVS(624) = 0.25*B(219)+B(408) ! JVS(625) = Jac_FULL(89,79) JVS(625) = 0.25*B(215)+B(406) ! JVS(626) = Jac_FULL(89,85) JVS(626) = 0.19*B(91)+0.35*B(203)+0.19*B(366)+B(396) ! JVS(627) = Jac_FULL(89,86) JVS(627) = 0.19*B(547) ! JVS(628) = Jac_FULL(89,89) JVS(628) = -B(197)-B(201)-B(518) ! JVS(629) = Jac_FULL(89,91) JVS(629) = 0 ! JVS(630) = Jac_FULL(89,92) JVS(630) = 0.25*B(211)+B(402) ! JVS(631) = Jac_FULL(89,93) JVS(631) = 0 ! JVS(632) = Jac_FULL(89,95) JVS(632) = 0.35*B(204)+0.25*B(210)+0.25*B(212)+0.25*B(214)+0.25*B(216)+0.25*B(220) ! JVS(633) = Jac_FULL(89,96) JVS(633) = 0 ! JVS(634) = Jac_FULL(89,98) JVS(634) = 0 ! JVS(635) = Jac_FULL(89,100) JVS(635) = -B(202) ! JVS(636) = Jac_FULL(89,101) JVS(636) = 0.19*B(367)+B(397)+B(399)+B(401)+B(403)+B(405)+B(407)+B(409)+B(411) ! JVS(637) = Jac_FULL(89,102) JVS(637) = -B(198) ! JVS(638) = Jac_FULL(89,103) JVS(638) = 0.19*B(92) ! JVS(639) = Jac_FULL(90,27) JVS(639) = B(139) ! JVS(640) = Jac_FULL(90,43) JVS(640) = 0.5*B(331) ! JVS(641) = Jac_FULL(90,81) JVS(641) = B(135)+B(157) ! JVS(642) = Jac_FULL(90,83) JVS(642) = 0 ! JVS(643) = Jac_FULL(90,85) JVS(643) = 0 ! JVS(644) = Jac_FULL(90,86) JVS(644) = 0 ! JVS(645) = Jac_FULL(90,89) JVS(645) = 0.15*B(518) ! JVS(646) = Jac_FULL(90,90) JVS(646) = -B(137)-B(149)-B(247)-B(303)-B(311)-B(430) ! JVS(647) = Jac_FULL(90,91) JVS(647) = 0 ! JVS(648) = Jac_FULL(90,92) JVS(648) = 0 ! JVS(649) = Jac_FULL(90,93) JVS(649) = 0 ! JVS(650) = Jac_FULL(90,95) JVS(650) = -B(304)-B(312) ! JVS(651) = Jac_FULL(90,96) JVS(651) = -B(248) ! JVS(652) = Jac_FULL(90,98) JVS(652) = 0 ! JVS(653) = Jac_FULL(90,99) JVS(653) = -B(138) ! JVS(654) = Jac_FULL(90,100) JVS(654) = B(158) ! JVS(655) = Jac_FULL(90,101) JVS(655) = -B(431) ! JVS(656) = Jac_FULL(90,102) JVS(656) = B(136)+0.5*B(332) ! JVS(657) = Jac_FULL(90,103) JVS(657) = -B(150) ! JVS(658) = Jac_FULL(91,57) JVS(658) = 0.259*B(541) ! JVS(659) = Jac_FULL(91,59) JVS(659) = 0.387*B(291) ! JVS(660) = Jac_FULL(91,70) JVS(660) = 0 ! JVS(661) = Jac_FULL(91,78) JVS(661) = 0.1*B(125)+0.05*B(229)+0.1*B(392) ! JVS(662) = Jac_FULL(91,91) JVS(662) = -B(279)-B(295)-B(299)-B(301)-B(528)-B(529) ! JVS(663) = Jac_FULL(91,92) JVS(663) = 0.22*B(101)+0.14*B(211)+0.288*B(372) ! JVS(664) = Jac_FULL(91,95) JVS(664) = 0.14*B(212)+0.05*B(230) ! JVS(665) = Jac_FULL(91,96) JVS(665) = 0 ! JVS(666) = Jac_FULL(91,98) JVS(666) = 0.387*B(292)-B(296) ! JVS(667) = Jac_FULL(91,100) JVS(667) = -B(300)-B(302) ! JVS(668) = Jac_FULL(91,101) JVS(668) = 0.288*B(373)+0.1*B(393) ! JVS(669) = Jac_FULL(91,102) JVS(669) = -B(280) ! JVS(670) = Jac_FULL(91,103) JVS(670) = 0.22*B(102)+0.1*B(126) ! JVS(671) = Jac_FULL(92,57) JVS(671) = 0.491*B(339) ! JVS(672) = Jac_FULL(92,59) JVS(672) = B(275) ! JVS(673) = Jac_FULL(92,70) JVS(673) = 0 ! JVS(674) = Jac_FULL(92,92) JVS(674) = -B(101)-B(175)-B(211)-B(372)-B(402) ! JVS(675) = Jac_FULL(92,95) JVS(675) = -B(212) ! JVS(676) = Jac_FULL(92,96) JVS(676) = -B(176) ! JVS(677) = Jac_FULL(92,98) JVS(677) = 0 ! JVS(678) = Jac_FULL(92,100) JVS(678) = 0 ! JVS(679) = Jac_FULL(92,101) JVS(679) = -B(373)-B(403) ! JVS(680) = Jac_FULL(92,102) JVS(680) = B(276)+0.491*B(340) ! JVS(681) = Jac_FULL(92,103) JVS(681) = -B(102) ! JVS(682) = Jac_FULL(93,57) JVS(682) = 0.368*B(541) ! JVS(683) = Jac_FULL(93,59) JVS(683) = 0.159*B(291) ! JVS(684) = Jac_FULL(93,70) JVS(684) = 0 ! JVS(685) = Jac_FULL(93,78) JVS(685) = 0.05*B(125)+0.03*B(229)+0.05*B(392) ! JVS(686) = Jac_FULL(93,92) JVS(686) = 0.34*B(101)+0.2*B(211)+0.402*B(372) ! JVS(687) = Jac_FULL(93,93) JVS(687) = -B(277)-B(293)-B(525)-B(526)-B(527) ! JVS(688) = Jac_FULL(93,95) JVS(688) = 0.2*B(212)+0.03*B(230) ! JVS(689) = Jac_FULL(93,96) JVS(689) = 0 ! JVS(690) = Jac_FULL(93,98) JVS(690) = 0.159*B(292)-B(294) ! JVS(691) = Jac_FULL(93,100) JVS(691) = 0 ! JVS(692) = Jac_FULL(93,101) JVS(692) = 0.402*B(373)+0.05*B(393) ! JVS(693) = Jac_FULL(93,102) JVS(693) = -B(278) ! JVS(694) = Jac_FULL(93,103) JVS(694) = 0.34*B(102)+0.05*B(126) ! JVS(695) = Jac_FULL(94,28) JVS(695) = B(142) ! JVS(696) = Jac_FULL(94,47) JVS(696) = B(335) ! JVS(697) = Jac_FULL(94,63) JVS(697) = 0.8*B(263) ! JVS(698) = Jac_FULL(94,64) JVS(698) = 0 ! JVS(699) = Jac_FULL(94,71) JVS(699) = 0 ! JVS(700) = Jac_FULL(94,75) JVS(700) = 0 ! JVS(701) = Jac_FULL(94,79) JVS(701) = 0 ! JVS(702) = Jac_FULL(94,80) JVS(702) = 0 ! JVS(703) = Jac_FULL(94,91) JVS(703) = 0 ! JVS(704) = Jac_FULL(94,92) JVS(704) = 0 ! JVS(705) = Jac_FULL(94,93) JVS(705) = 0 ! JVS(706) = Jac_FULL(94,94) JVS(706) = -B(140)-B(151)-B(249)-B(305)-B(313)-B(432) ! JVS(707) = Jac_FULL(94,95) JVS(707) = -B(306)-B(314) ! JVS(708) = Jac_FULL(94,96) JVS(708) = -B(250) ! JVS(709) = Jac_FULL(94,98) JVS(709) = 0 ! JVS(710) = Jac_FULL(94,99) JVS(710) = -B(141) ! JVS(711) = Jac_FULL(94,100) JVS(711) = 0 ! JVS(712) = Jac_FULL(94,101) JVS(712) = -B(433) ! JVS(713) = Jac_FULL(94,102) JVS(713) = 0.8*B(264)+B(336) ! JVS(714) = Jac_FULL(94,103) JVS(714) = -B(152) ! JVS(715) = Jac_FULL(95,34) JVS(715) = 0.4*B(514) ! JVS(716) = Jac_FULL(95,39) JVS(716) = B(32) ! JVS(717) = Jac_FULL(95,44) JVS(717) = B(471)+B(473)+B(475) ! JVS(718) = Jac_FULL(95,48) JVS(718) = B(548) ! JVS(719) = Jac_FULL(95,61) JVS(719) = 0.305*B(257) ! JVS(720) = Jac_FULL(95,62) JVS(720) = B(516)+2*B(517) ! JVS(721) = Jac_FULL(95,64) JVS(721) = -B(223)+B(388) ! JVS(722) = Jac_FULL(95,65) JVS(722) = -B(163)+B(283) ! JVS(723) = Jac_FULL(95,66) JVS(723) = -B(227)+B(384) ! JVS(724) = Jac_FULL(95,68) JVS(724) = -B(205)+B(386) ! JVS(725) = Jac_FULL(95,69) JVS(725) = -B(225)+B(390) ! JVS(726) = Jac_FULL(95,70) JVS(726) = -B(213)+B(374) ! JVS(727) = Jac_FULL(95,72) JVS(727) = -B(221)+B(382) ! JVS(728) = Jac_FULL(95,73) JVS(728) = B(173)-B(209)+B(370) ! JVS(729) = Jac_FULL(95,74) JVS(729) = B(171)-B(207)+B(368) ! JVS(730) = Jac_FULL(95,75) JVS(730) = -B(219)+B(380) ! JVS(731) = Jac_FULL(95,76) JVS(731) = -B(231)+B(394) ! JVS(732) = Jac_FULL(95,78) JVS(732) = -B(229)+B(392) ! JVS(733) = Jac_FULL(95,79) JVS(733) = -B(215)+B(376) ! JVS(734) = Jac_FULL(95,80) JVS(734) = -B(217)+B(378) ! JVS(735) = Jac_FULL(95,83) JVS(735) = -B(165)+B(285) ! JVS(736) = Jac_FULL(95,84) JVS(736) = B(512) ! JVS(737) = Jac_FULL(95,85) JVS(737) = 0.18*B(91)-0.91*B(203)+1.18*B(366) ! JVS(738) = Jac_FULL(95,86) JVS(738) = 0.18*B(547) ! JVS(739) = Jac_FULL(95,87) JVS(739) = -B(199)+B(426) ! JVS(740) = Jac_FULL(95,88) JVS(740) = 0 ! JVS(741) = Jac_FULL(95,89) JVS(741) = 0.15*B(518) ! JVS(742) = Jac_FULL(95,90) JVS(742) = -B(303)-B(311)+B(430) ! JVS(743) = Jac_FULL(95,91) JVS(743) = 0 ! JVS(744) = Jac_FULL(95,92) JVS(744) = -B(211)+B(372) ! JVS(745) = Jac_FULL(95,93) JVS(745) = B(527) ! JVS(746) = Jac_FULL(95,94) JVS(746) = -B(305)-B(313)+B(432) ! JVS(747) = Jac_FULL(95,95) JVS(747) = -B(9)-B(26)-B(28)-2*B(30)-2*B(31)-B(164)-B(166)-B(200)-0.91*B(204)-B(206)-B(208)-B(210)-B(212)-B(214)& &-B(216)-B(218)-B(220)-B(222)-B(224)-B(226)-B(228)-B(230)-B(232)-B(304)-B(306)-B(307)-B(309)-B(312)-B(314)& &-B(315)-B(317)-B(364) ! JVS(748) = Jac_FULL(95,96) JVS(748) = -B(29)+B(172)+B(174)+0.44*B(245) ! JVS(749) = Jac_FULL(95,97) JVS(749) = -B(308)-B(316)+B(434) ! JVS(750) = Jac_FULL(95,98) JVS(750) = -B(10)+0.305*B(258) ! JVS(751) = Jac_FULL(95,99) JVS(751) = 0 ! JVS(752) = Jac_FULL(95,100) JVS(752) = B(476) ! JVS(753) = Jac_FULL(95,101) JVS(753) = B(75)+0.44*B(246)+B(284)+B(286)+2*B(361)-B(365)+1.18*B(367)+B(369)+B(371)+B(373)+B(375)+B(377)+B(379)& &+B(381)+B(383)+B(385)+B(387)+B(389)+B(391)+B(393)+B(395)+B(427)+B(431)+B(433)+B(435)+B(436) ! JVS(754) = Jac_FULL(95,102) JVS(754) = B(24)+B(33)+B(133)+B(472)+B(474) ! JVS(755) = Jac_FULL(95,103) JVS(755) = -B(27)+B(76)+0.18*B(92) ! JVS(756) = Jac_FULL(96,29) JVS(756) = B(15) ! JVS(757) = Jac_FULL(96,30) JVS(757) = B(477) ! JVS(758) = Jac_FULL(96,31) JVS(758) = -B(486) ! JVS(759) = Jac_FULL(96,32) JVS(759) = -B(490) ! JVS(760) = Jac_FULL(96,33) JVS(760) = -B(494) ! JVS(761) = Jac_FULL(96,39) JVS(761) = B(501) ! JVS(762) = Jac_FULL(96,40) JVS(762) = B(536) ! JVS(763) = Jac_FULL(96,41) JVS(763) = B(534) ! JVS(764) = Jac_FULL(96,42) JVS(764) = B(535) ! JVS(765) = Jac_FULL(96,43) JVS(765) = B(538) ! JVS(766) = Jac_FULL(96,44) JVS(766) = 0 ! JVS(767) = Jac_FULL(96,46) JVS(767) = B(533) ! JVS(768) = Jac_FULL(96,47) JVS(768) = B(539) ! JVS(769) = Jac_FULL(96,49) JVS(769) = B(48)+B(511) ! JVS(770) = Jac_FULL(96,50) JVS(770) = B(537) ! JVS(771) = Jac_FULL(96,51) JVS(771) = B(532) ! JVS(772) = Jac_FULL(96,52) JVS(772) = B(531) ! JVS(773) = Jac_FULL(96,53) JVS(773) = B(545) ! JVS(774) = Jac_FULL(96,54) JVS(774) = B(542) ! JVS(775) = Jac_FULL(96,55) JVS(775) = 0.3*B(544) ! JVS(776) = Jac_FULL(96,56) JVS(776) = B(543) ! JVS(777) = Jac_FULL(96,57) JVS(777) = B(541) ! JVS(778) = Jac_FULL(96,58) JVS(778) = 2*B(259)+B(261) ! JVS(779) = Jac_FULL(96,59) JVS(779) = 0.06*B(291) ! JVS(780) = Jac_FULL(96,60) JVS(780) = B(22) ! JVS(781) = Jac_FULL(96,61) JVS(781) = 0.3*B(257) ! JVS(782) = Jac_FULL(96,62) JVS(782) = 0 ! JVS(783) = Jac_FULL(96,63) JVS(783) = 0.2*B(263)+2*B(520) ! JVS(784) = Jac_FULL(96,64) JVS(784) = 0.3*B(119)-B(187)+0.75*B(223)+0.5*B(388) ! JVS(785) = Jac_FULL(96,65) JVS(785) = B(85)+B(163)-B(241)+B(283) ! JVS(786) = Jac_FULL(96,66) JVS(786) = B(123)-B(191)+B(227)+B(384) ! JVS(787) = Jac_FULL(96,68) JVS(787) = -B(169)+0.5*B(205) ! JVS(788) = Jac_FULL(96,69) JVS(788) = -B(189)+0.5*B(225) ! JVS(789) = Jac_FULL(96,70) JVS(789) = B(103)-B(177)+B(213)+B(374) ! JVS(790) = Jac_FULL(96,71) JVS(790) = 0.15*B(357) ! JVS(791) = Jac_FULL(96,72) JVS(791) = B(115)-B(185)+1.15*B(221)+B(382) ! JVS(792) = Jac_FULL(96,73) JVS(792) = -B(173)+0.5*B(209) ! JVS(793) = Jac_FULL(96,74) JVS(793) = -B(171)+0.3*B(207)+0.8*B(368) ! JVS(794) = Jac_FULL(96,75) JVS(794) = 0.28*B(111)-B(183)+0.64*B(219)+0.28*B(380) ! JVS(795) = Jac_FULL(96,76) JVS(795) = -B(195)+0.5*B(231) ! JVS(796) = Jac_FULL(96,77) JVS(796) = B(281)+B(530) ! JVS(797) = Jac_FULL(96,78) JVS(797) = 0.8*B(125)-B(193)+0.45*B(229)+0.8*B(392) ! JVS(798) = Jac_FULL(96,79) JVS(798) = 0.92*B(107)-B(179)+B(215)+B(376) ! JVS(799) = Jac_FULL(96,80) JVS(799) = 0.05*B(109)-B(181)+0.5*B(217) ! JVS(800) = Jac_FULL(96,81) JVS(800) = B(515) ! JVS(801) = Jac_FULL(96,82) JVS(801) = B(36)+B(66)+2*B(502) ! JVS(802) = Jac_FULL(96,83) JVS(802) = B(87)+B(165)-B(243)+B(285) ! JVS(803) = Jac_FULL(96,84) JVS(803) = 0.05*B(68)+B(512) ! JVS(804) = Jac_FULL(96,85) JVS(804) = 0.27*B(91)-B(167)+0.64*B(203)+0.27*B(366) ! JVS(805) = Jac_FULL(96,86) JVS(805) = 0.27*B(547) ! JVS(806) = Jac_FULL(96,87) JVS(806) = B(79)+B(199)+2*B(237)-B(239)+B(426) ! JVS(807) = Jac_FULL(96,88) JVS(807) = B(523) ! JVS(808) = Jac_FULL(96,89) JVS(808) = 0 ! JVS(809) = Jac_FULL(96,90) JVS(809) = -B(247)+B(303) ! JVS(810) = Jac_FULL(96,91) JVS(810) = 0.275*B(295)+B(528)+B(529) ! JVS(811) = Jac_FULL(96,92) JVS(811) = 0.9*B(101)-B(175)+0.92*B(211)+0.864*B(372) ! JVS(812) = Jac_FULL(96,93) JVS(812) = 0.06*B(293)+B(526) ! JVS(813) = Jac_FULL(96,94) JVS(813) = B(151)-B(249)+2*B(305)+B(432) ! JVS(814) = Jac_FULL(96,95) JVS(814) = B(9)+B(26)-B(28)+2*B(31)+B(164)+B(166)+B(200)+0.64*B(204)+0.5*B(206)+0.3*B(208)+0.5*B(210)+0.92*B(212)& &+B(214)+B(216)+0.5*B(218)+0.64*B(220)+1.15*B(222)+0.75*B(224)+0.5*B(226)+B(228)+0.45*B(230)+0.5*B(232)& &+B(304)+2*B(306)+B(307)+2*B(309)+B(362) ! JVS(815) = Jac_FULL(96,96) JVS(815) = -B(5)-B(13)-B(17)-2*B(19)-B(29)-B(46)-B(51)-B(168)-B(170)-B(172)-B(174)-B(176)-B(178)-B(180)-B(182)-B(184)& &-B(186)-B(188)-B(190)-B(192)-B(194)-B(196)-B(240)-B(242)-B(244)-B(245)-B(248)-B(250)-B(251)-B(253)-B(467)& &-B(487)-B(491)-B(495) ! JVS(816) = Jac_FULL(96,97) JVS(816) = -B(252)+B(308) ! JVS(817) = Jac_FULL(96,98) JVS(817) = B(3)-B(6)+B(10)+0.3*B(258)+B(262)+0.06*B(292)+0.06*B(294)+0.275*B(296) ! JVS(818) = Jac_FULL(96,99) JVS(818) = -B(47) ! JVS(819) = Jac_FULL(96,100) JVS(819) = -B(52)+B(55)+B(67)+B(271) ! JVS(820) = Jac_FULL(96,101) JVS(820) = -B(246)+B(284)+B(286)+B(363)+0.27*B(367)+0.8*B(369)+0.864*B(373)+B(375)+B(377)+0.28*B(381)+B(383)+B(385)& &+0.5*B(389)+0.8*B(393)+B(427)+B(433)+B(436) ! JVS(821) = Jac_FULL(96,102) JVS(821) = B(4)-B(14)+B(16)+B(20)+B(23)+B(37)+B(56)+B(60)+B(62)+0.05*B(69)+B(233)+B(235)+2*B(260)+0.2*B(264)+B(267)& &+B(282)+0.15*B(358)+B(478) ! JVS(822) = Jac_FULL(96,103) JVS(822) = -B(18)+B(27)+B(80)+B(86)+B(88)+0.27*B(92)+0.9*B(102)+B(104)+0.92*B(108)+0.05*B(110)+0.28*B(112)+B(116)+0.3& &*B(120)+B(124)+0.8*B(126)+B(152)+B(155) ! JVS(823) = Jac_FULL(97,45) JVS(823) = B(349) ! JVS(824) = Jac_FULL(97,58) JVS(824) = B(145) ! JVS(825) = Jac_FULL(97,71) JVS(825) = 0.41*B(357) ! JVS(826) = Jac_FULL(97,91) JVS(826) = 0.57*B(279)+B(301)+B(528) ! JVS(827) = Jac_FULL(97,92) JVS(827) = 0 ! JVS(828) = Jac_FULL(97,93) JVS(828) = B(527) ! JVS(829) = Jac_FULL(97,95) JVS(829) = -B(307)-B(315) ! JVS(830) = Jac_FULL(97,96) JVS(830) = -B(251) ! JVS(831) = Jac_FULL(97,97) JVS(831) = -B(143)-B(153)-B(252)-B(308)-B(316)-B(434) ! JVS(832) = Jac_FULL(97,98) JVS(832) = 0 ! JVS(833) = Jac_FULL(97,99) JVS(833) = -B(144) ! JVS(834) = Jac_FULL(97,100) JVS(834) = B(302) ! JVS(835) = Jac_FULL(97,101) JVS(835) = -B(435) ! JVS(836) = Jac_FULL(97,102) JVS(836) = 0.57*B(280)+B(350)+0.41*B(358) ! JVS(837) = Jac_FULL(97,103) JVS(837) = -B(154) ! JVS(838) = Jac_FULL(98,38) JVS(838) = B(510) ! JVS(839) = Jac_FULL(98,58) JVS(839) = -B(261) ! JVS(840) = Jac_FULL(98,59) JVS(840) = -0.9*B(291) ! JVS(841) = Jac_FULL(98,61) JVS(841) = -B(257) ! JVS(842) = Jac_FULL(98,71) JVS(842) = -0.7*B(359) ! JVS(843) = Jac_FULL(98,90) JVS(843) = 0.3*B(247) ! JVS(844) = Jac_FULL(98,91) JVS(844) = -0.8*B(295) ! JVS(845) = Jac_FULL(98,92) JVS(845) = 0 ! JVS(846) = Jac_FULL(98,93) JVS(846) = -0.8*B(293) ! JVS(847) = Jac_FULL(98,94) JVS(847) = 0.29*B(249) ! JVS(848) = Jac_FULL(98,95) JVS(848) = -B(9) ! JVS(849) = Jac_FULL(98,96) JVS(849) = -B(5)+0.15*B(245)+0.3*B(248)+0.29*B(250)+0.3*B(251)+0.3*B(253) ! JVS(850) = Jac_FULL(98,97) JVS(850) = 0.3*B(252) ! JVS(851) = Jac_FULL(98,98) JVS(851) = -B(1)-B(3)-B(6)-B(7)-B(10)-B(258)-B(262)-0.9*B(292)-0.8*B(294)-0.8*B(296)-0.7*B(360)-B(458)-B(498) ! JVS(852) = Jac_FULL(98,99) JVS(852) = -B(8)+B(499) ! JVS(853) = Jac_FULL(98,100) JVS(853) = B(507) ! JVS(854) = Jac_FULL(98,101) JVS(854) = 0.15*B(246) ! JVS(855) = Jac_FULL(98,102) JVS(855) = -B(4)+B(11) ! JVS(856) = Jac_FULL(98,103) JVS(856) = -B(2) ! JVS(857) = Jac_FULL(99,27) JVS(857) = B(139) ! JVS(858) = Jac_FULL(99,28) JVS(858) = B(142) ! JVS(859) = Jac_FULL(99,34) JVS(859) = B(74)+0.6*B(514) ! JVS(860) = Jac_FULL(99,37) JVS(860) = B(44) ! JVS(861) = Jac_FULL(99,38) JVS(861) = B(59)+B(509) ! JVS(862) = Jac_FULL(99,49) JVS(862) = B(48)+B(49)+B(511) ! JVS(863) = Jac_FULL(99,51) JVS(863) = B(532) ! JVS(864) = Jac_FULL(99,52) JVS(864) = B(531) ! JVS(865) = Jac_FULL(99,56) JVS(865) = 0.5*B(343)+B(543) ! JVS(866) = Jac_FULL(99,58) JVS(866) = B(145)+B(259)+B(261) ! JVS(867) = Jac_FULL(99,64) JVS(867) = 1.9*B(119)+B(223)+B(388)+B(416) ! JVS(868) = Jac_FULL(99,65) JVS(868) = B(85) ! JVS(869) = Jac_FULL(99,66) JVS(869) = B(123) ! JVS(870) = Jac_FULL(99,67) JVS(870) = B(504) ! JVS(871) = Jac_FULL(99,68) JVS(871) = 2*B(95)+B(205)+B(386)+B(412) ! JVS(872) = Jac_FULL(99,69) JVS(872) = 2*B(121)+B(225)+B(390)+B(418) ! JVS(873) = Jac_FULL(99,70) JVS(873) = B(103) ! JVS(874) = Jac_FULL(99,72) JVS(874) = B(115) ! JVS(875) = Jac_FULL(99,73) JVS(875) = 0.93*B(99) ! JVS(876) = Jac_FULL(99,74) JVS(876) = 0.96*B(97) ! JVS(877) = Jac_FULL(99,75) JVS(877) = B(111) ! JVS(878) = Jac_FULL(99,76) JVS(878) = 2*B(127)+B(231)+B(394)+B(422) ! JVS(879) = Jac_FULL(99,78) JVS(879) = 1.15*B(125)+0.575*B(229)+0.15*B(392)+B(420) ! JVS(880) = Jac_FULL(99,79) JVS(880) = 0.92*B(107) ! JVS(881) = Jac_FULL(99,80) JVS(881) = 1.95*B(109)+B(217)+B(378)+B(414) ! JVS(882) = Jac_FULL(99,81) JVS(882) = 0 ! JVS(883) = Jac_FULL(99,82) JVS(883) = 0 ! JVS(884) = Jac_FULL(99,83) JVS(884) = B(87) ! JVS(885) = Jac_FULL(99,84) JVS(885) = 0 ! JVS(886) = Jac_FULL(99,85) JVS(886) = B(91) ! JVS(887) = Jac_FULL(99,86) JVS(887) = B(547) ! JVS(888) = Jac_FULL(99,87) JVS(888) = B(79) ! JVS(889) = Jac_FULL(99,88) JVS(889) = 0 ! JVS(890) = Jac_FULL(99,89) JVS(890) = 0 ! JVS(891) = Jac_FULL(99,90) JVS(891) = -B(137)+B(149) ! JVS(892) = Jac_FULL(99,91) JVS(892) = 0 ! JVS(893) = Jac_FULL(99,92) JVS(893) = 0.9*B(101) ! JVS(894) = Jac_FULL(99,93) JVS(894) = 0 ! JVS(895) = Jac_FULL(99,94) JVS(895) = -B(140)+B(151) ! JVS(896) = Jac_FULL(99,95) JVS(896) = B(26)+B(206)+B(218)+B(224)+B(226)+0.575*B(230)+B(232) ! JVS(897) = Jac_FULL(99,96) JVS(897) = B(17)-B(46)+B(51) ! JVS(898) = Jac_FULL(99,97) JVS(898) = -B(143)+B(153) ! JVS(899) = Jac_FULL(99,98) JVS(899) = B(1)-B(7)+B(262) ! JVS(900) = Jac_FULL(99,99) JVS(900) = -B(8)-B(38)-B(47)-B(57)-B(72)-B(138)-B(141)-B(144)-B(146)-B(457)-B(468)-B(499) ! JVS(901) = Jac_FULL(99,100) JVS(901) = B(52)+2*B(53)+B(55)-B(58)+2*B(438)+B(507) ! JVS(902) = Jac_FULL(99,101) JVS(902) = -B(73)+B(75)+B(379)+B(387)+B(389)+B(391)+0.15*B(393)+B(395)+B(413)+B(415)+B(417)+B(419)+B(421)+B(423) ! JVS(903) = Jac_FULL(99,102) JVS(903) = -B(39)+B(45)+B(50)+B(56)+B(260)+0.5*B(344)+B(355) ! JVS(904) = Jac_FULL(99,103) JVS(904) = B(2)+B(18)+B(27)+2*B(54)+B(76)+B(80)+B(86)+B(88)+B(92)+2*B(96)+0.96*B(98)+0.93*B(100)+0.9*B(102)+B(104)& &+0.92*B(108)+1.95*B(110)+B(112)+B(116)+1.9*B(120)+2*B(122)+B(124)+1.15*B(126)+2*B(128)+B(150)+B(152)+B(154)& &+B(155) ! JVS(905) = Jac_FULL(100,34) JVS(905) = 0.4*B(514) ! JVS(906) = Jac_FULL(100,35) JVS(906) = -B(129) ! JVS(907) = Jac_FULL(100,36) JVS(907) = -B(353) ! JVS(908) = Jac_FULL(100,38) JVS(908) = B(59)+B(509)+B(510) ! JVS(909) = Jac_FULL(100,44) JVS(909) = -B(475) ! JVS(910) = Jac_FULL(100,49) JVS(910) = B(506) ! JVS(911) = Jac_FULL(100,59) JVS(911) = -B(297) ! JVS(912) = Jac_FULL(100,61) JVS(912) = -B(265) ! JVS(913) = Jac_FULL(100,67) JVS(913) = B(40) ! JVS(914) = Jac_FULL(100,70) JVS(914) = 0 ! JVS(915) = Jac_FULL(100,72) JVS(915) = 0 ! JVS(916) = Jac_FULL(100,75) JVS(916) = 0 ! JVS(917) = Jac_FULL(100,78) JVS(917) = 0 ! JVS(918) = Jac_FULL(100,79) JVS(918) = 0 ! JVS(919) = Jac_FULL(100,80) JVS(919) = 0 ! JVS(920) = Jac_FULL(100,81) JVS(920) = -B(157) ! JVS(921) = Jac_FULL(100,82) JVS(921) = -B(66) ! JVS(922) = Jac_FULL(100,83) JVS(922) = 0 ! JVS(923) = Jac_FULL(100,84) JVS(923) = -B(70) ! JVS(924) = Jac_FULL(100,85) JVS(924) = 0 ! JVS(925) = Jac_FULL(100,86) JVS(925) = 0 ! JVS(926) = Jac_FULL(100,87) JVS(926) = 0 ! JVS(927) = Jac_FULL(100,88) JVS(927) = -B(273) ! JVS(928) = Jac_FULL(100,89) JVS(928) = -B(201) ! JVS(929) = Jac_FULL(100,90) JVS(929) = 0 ! JVS(930) = Jac_FULL(100,91) JVS(930) = -B(299)-B(301) ! JVS(931) = Jac_FULL(100,92) JVS(931) = 0 ! JVS(932) = Jac_FULL(100,93) JVS(932) = 0 ! JVS(933) = Jac_FULL(100,94) JVS(933) = 0 ! JVS(934) = Jac_FULL(100,95) JVS(934) = 0 ! JVS(935) = Jac_FULL(100,96) JVS(935) = -B(51) ! JVS(936) = Jac_FULL(100,97) JVS(936) = 0 ! JVS(937) = Jac_FULL(100,98) JVS(937) = B(7) ! JVS(938) = Jac_FULL(100,99) JVS(938) = B(8)-B(57)-B(64) ! JVS(939) = Jac_FULL(100,100) JVS(939) = -B(52)-B(53)-B(55)-B(58)-B(65)-B(67)-B(71)-B(130)-B(158)-B(202)-B(266)-B(271)-B(274)-B(298)-B(300)-B(302)& &-B(354)-2*B(438)-B(469)-B(476)-B(507)-B(508) ! JVS(940) = Jac_FULL(100,101) JVS(940) = 0 ! JVS(941) = Jac_FULL(100,102) JVS(941) = B(41)-B(56) ! JVS(942) = Jac_FULL(100,103) JVS(942) = -B(54) ! JVS(943) = Jac_FULL(101,34) JVS(943) = B(74)+0.6*B(514) ! JVS(944) = Jac_FULL(101,45) JVS(944) = B(546) ! JVS(945) = Jac_FULL(101,48) JVS(945) = 0.5*B(351) ! JVS(946) = Jac_FULL(101,55) JVS(946) = 0.7*B(544) ! JVS(947) = Jac_FULL(101,62) JVS(947) = B(516) ! JVS(948) = Jac_FULL(101,64) JVS(948) = 0.6*B(119)+0.25*B(223)-0.5*B(388)-B(416) ! JVS(949) = Jac_FULL(101,65) JVS(949) = -B(283)-B(287) ! JVS(950) = Jac_FULL(101,66) JVS(950) = -B(384)-B(424) ! JVS(951) = Jac_FULL(101,68) JVS(951) = -B(386)-B(412) ! JVS(952) = Jac_FULL(101,69) JVS(952) = -B(390)-B(418) ! JVS(953) = Jac_FULL(101,70) JVS(953) = -B(374)-B(404) ! JVS(954) = Jac_FULL(101,72) JVS(954) = -B(382)-B(410) ! JVS(955) = Jac_FULL(101,73) JVS(955) = 0.93*B(99)+0.5*B(209)-B(400) ! JVS(956) = Jac_FULL(101,74) JVS(956) = 0.96*B(97)+B(171)+0.3*B(207)-0.8*B(368)-B(398) ! JVS(957) = Jac_FULL(101,75) JVS(957) = 0.72*B(111)+0.36*B(219)-0.28*B(380)-B(408) ! JVS(958) = Jac_FULL(101,76) JVS(958) = -B(394)-B(422) ! JVS(959) = Jac_FULL(101,77) JVS(959) = B(530) ! JVS(960) = Jac_FULL(101,78) JVS(960) = -B(392)-B(420) ! JVS(961) = Jac_FULL(101,79) JVS(961) = -B(376)-B(406) ! JVS(962) = Jac_FULL(101,80) JVS(962) = -B(378)-B(414) ! JVS(963) = Jac_FULL(101,83) JVS(963) = -B(285)-B(289) ! JVS(964) = Jac_FULL(101,84) JVS(964) = 0.95*B(68)+B(70) ! JVS(965) = Jac_FULL(101,85) JVS(965) = -B(366)-B(396) ! JVS(966) = Jac_FULL(101,86) JVS(966) = 0 ! JVS(967) = Jac_FULL(101,87) JVS(967) = -B(426)-B(428) ! JVS(968) = Jac_FULL(101,88) JVS(968) = B(269)+B(273)+B(523) ! JVS(969) = Jac_FULL(101,89) JVS(969) = 0.85*B(518) ! JVS(970) = Jac_FULL(101,90) JVS(970) = -B(430) ! JVS(971) = Jac_FULL(101,91) JVS(971) = B(529) ! JVS(972) = Jac_FULL(101,92) JVS(972) = -B(372)-B(402) ! JVS(973) = Jac_FULL(101,93) JVS(973) = B(526) ! JVS(974) = Jac_FULL(101,94) JVS(974) = -B(432) ! JVS(975) = Jac_FULL(101,95) JVS(975) = 0.3*B(208)+0.5*B(210)+0.36*B(220)+0.25*B(224)+B(307)-B(362)-B(364) ! JVS(976) = Jac_FULL(101,96) JVS(976) = B(172)-B(245) ! JVS(977) = Jac_FULL(101,97) JVS(977) = B(153)+B(308) ! JVS(978) = Jac_FULL(101,98) JVS(978) = 0 ! JVS(979) = Jac_FULL(101,99) JVS(979) = -B(72) ! JVS(980) = Jac_FULL(101,100) JVS(980) = B(71)+B(274) ! JVS(981) = Jac_FULL(101,101) JVS(981) = -B(73)-B(75)-B(246)-B(284)-B(286)-B(288)-B(290)-2*B(361)-B(363)-B(365)-B(367)-0.8*B(369)-B(373)-B(375)& &-B(377)-B(379)-0.28*B(381)-B(383)-B(385)-B(387)-0.5*B(389)-B(391)-B(393)-B(395)-B(397)-B(399)-B(401)-B(403)& &-B(405)-B(407)-B(409)-B(411)-B(413)-B(415)-B(417)-B(419)-B(421)-B(423)-B(425)-B(427)-B(429)-B(431)-B(433)& &-B(436) ! JVS(982) = Jac_FULL(101,102) JVS(982) = 0.95*B(69)+B(270)+0.5*B(352) ! JVS(983) = Jac_FULL(101,103) JVS(983) = -B(76)+0.96*B(98)+0.93*B(100)+0.72*B(112)+0.6*B(120)+B(154) ! JVS(984) = Jac_FULL(102,23) JVS(984) = -B(481) ! JVS(985) = Jac_FULL(102,24) JVS(985) = -B(483) ! JVS(986) = Jac_FULL(102,25) JVS(986) = -B(479) ! JVS(987) = Jac_FULL(102,26) JVS(987) = -B(81)-B(83) ! JVS(988) = Jac_FULL(102,29) JVS(988) = -B(15)+2*B(500) ! JVS(989) = Jac_FULL(102,30) JVS(989) = -B(477) ! JVS(990) = Jac_FULL(102,35) JVS(990) = -B(89) ! JVS(991) = Jac_FULL(102,36) JVS(991) = -B(77) ! JVS(992) = Jac_FULL(102,37) JVS(992) = -B(44)+B(505) ! JVS(993) = Jac_FULL(102,39) JVS(993) = -B(32)+B(501) ! JVS(994) = Jac_FULL(102,40) JVS(994) = -0.5*B(329)+B(536) ! JVS(995) = Jac_FULL(102,41) JVS(995) = -0.5*B(325)+B(534) ! JVS(996) = Jac_FULL(102,42) JVS(996) = -0.5*B(327)+B(535) ! JVS(997) = Jac_FULL(102,43) JVS(997) = -0.5*B(331)+B(538) ! JVS(998) = Jac_FULL(102,44) JVS(998) = -B(471)-B(473) ! JVS(999) = Jac_FULL(102,45) JVS(999) = -B(349)+B(546) ! JVS(1000) = Jac_FULL(102,46) JVS(1000) = -0.5*B(323)+B(533) ! JVS(1001) = Jac_FULL(102,47) JVS(1001) = -B(335)+B(539) ! JVS(1002) = Jac_FULL(102,48) JVS(1002) = -0.5*B(351)+B(548) ! JVS(1003) = Jac_FULL(102,49) JVS(1003) = -B(49)+B(506) ! JVS(1004) = Jac_FULL(102,50) JVS(1004) = -B(333)+B(537) ! JVS(1005) = Jac_FULL(102,51) JVS(1005) = -B(321)+B(532) ! JVS(1006) = Jac_FULL(102,52) JVS(1006) = -B(319)+B(531) ! JVS(1007) = Jac_FULL(102,53) JVS(1007) = -B(347)+B(545) ! JVS(1008) = Jac_FULL(102,54) JVS(1008) = -B(341)+B(542) ! JVS(1009) = Jac_FULL(102,55) JVS(1009) = -0.5*B(345)+B(544) ! JVS(1010) = Jac_FULL(102,56) JVS(1010) = -0.5*B(343)+B(543) ! JVS(1011) = Jac_FULL(102,57) JVS(1011) = -0.491*B(339)+B(541) ! JVS(1012) = Jac_FULL(102,58) JVS(1012) = -B(259) ! JVS(1013) = Jac_FULL(102,59) JVS(1013) = -B(275)+0.27*B(291) ! JVS(1014) = Jac_FULL(102,60) JVS(1014) = -B(22) ! JVS(1015) = Jac_FULL(102,61) JVS(1015) = -B(255)+0.135*B(257) ! JVS(1016) = Jac_FULL(102,62) JVS(1016) = -B(159)-B(161) ! JVS(1017) = Jac_FULL(102,63) JVS(1017) = -B(263) ! JVS(1018) = Jac_FULL(102,64) JVS(1018) = 0 ! JVS(1019) = Jac_FULL(102,65) JVS(1019) = 0 ! JVS(1020) = Jac_FULL(102,66) JVS(1020) = 0 ! JVS(1021) = Jac_FULL(102,67) JVS(1021) = -B(40)+B(504) ! JVS(1022) = Jac_FULL(102,69) JVS(1022) = 0 ! JVS(1023) = Jac_FULL(102,70) JVS(1023) = 0 ! JVS(1024) = Jac_FULL(102,71) JVS(1024) = -B(357)+0.1*B(359) ! JVS(1025) = Jac_FULL(102,72) JVS(1025) = 0 ! JVS(1026) = Jac_FULL(102,75) JVS(1026) = 0 ! JVS(1027) = Jac_FULL(102,76) JVS(1027) = 0 ! JVS(1028) = Jac_FULL(102,77) JVS(1028) = -B(281) ! JVS(1029) = Jac_FULL(102,78) JVS(1029) = 0 ! JVS(1030) = Jac_FULL(102,79) JVS(1030) = 0 ! JVS(1031) = Jac_FULL(102,80) JVS(1031) = 0 ! JVS(1032) = Jac_FULL(102,81) JVS(1032) = -B(135) ! JVS(1033) = Jac_FULL(102,82) JVS(1033) = -B(36) ! JVS(1034) = Jac_FULL(102,83) JVS(1034) = 0 ! JVS(1035) = Jac_FULL(102,84) JVS(1035) = -B(68) ! JVS(1036) = Jac_FULL(102,85) JVS(1036) = 0 ! JVS(1037) = Jac_FULL(102,86) JVS(1037) = -B(131) ! JVS(1038) = Jac_FULL(102,87) JVS(1038) = 0 ! JVS(1039) = Jac_FULL(102,88) JVS(1039) = -B(269) ! JVS(1040) = Jac_FULL(102,89) JVS(1040) = -B(197) ! JVS(1041) = Jac_FULL(102,90) JVS(1041) = 0 ! JVS(1042) = Jac_FULL(102,91) JVS(1042) = -B(279)+0.215*B(295) ! JVS(1043) = Jac_FULL(102,92) JVS(1043) = 0 ! JVS(1044) = Jac_FULL(102,93) JVS(1044) = -B(277)+0.08*B(293) ! JVS(1045) = Jac_FULL(102,94) JVS(1045) = 0 ! JVS(1046) = Jac_FULL(102,95) JVS(1046) = 0 ! JVS(1047) = Jac_FULL(102,96) JVS(1047) = B(5)-B(13)+B(17)+B(51)+0.44*B(245) ! JVS(1048) = Jac_FULL(102,97) JVS(1048) = 0 ! JVS(1049) = Jac_FULL(102,98) JVS(1049) = -B(3)+B(6)+0.135*B(258)+0.27*B(292)+0.08*B(294)+0.215*B(296)+0.1*B(360)+2*B(498) ! JVS(1050) = Jac_FULL(102,99) JVS(1050) = -B(38) ! JVS(1051) = Jac_FULL(102,100) JVS(1051) = B(52)-B(55) ! JVS(1052) = Jac_FULL(102,101) JVS(1052) = 0.44*B(246) ! JVS(1053) = Jac_FULL(102,102) JVS(1053) = -B(4)-2*B(11)-2*B(12)-B(14)-B(16)-B(20)-B(23)-B(24)-B(33)-B(37)-B(39)-B(41)-B(42)-B(45)-B(50)-B(56)-B(60)& &-B(62)-B(69)-B(78)-B(82)-B(84)-B(90)-B(132)-B(133)-B(136)-B(160)-B(162)-B(198)-B(233)-B(235)-B(256)-B(260)& &-B(264)-B(267)-B(270)-B(276)-B(278)-B(280)-B(282)-B(320)-B(322)-0.5*B(324)-0.5*B(326)-0.5*B(328)-0.5*B(330)& &-0.5*B(332)-B(334)-B(336)-0.5*B(337)-0.491*B(340)-B(342)-0.5*B(344)-0.5*B(346)-B(348)-B(350)-0.5*B(352)& &-B(355)-B(358)-B(472)-B(474)-B(478)-B(480)-B(482)-B(484) ! JVS(1054) = Jac_FULL(102,103) JVS(1054) = B(18)-B(43) ! JVS(1055) = Jac_FULL(103,31) JVS(1055) = -B(488) ! JVS(1056) = Jac_FULL(103,32) JVS(1056) = -B(492) ! JVS(1057) = Jac_FULL(103,33) JVS(1057) = -B(496) ! JVS(1058) = Jac_FULL(103,37) JVS(1058) = B(505) ! JVS(1059) = Jac_FULL(103,38) JVS(1059) = B(510) ! JVS(1060) = Jac_FULL(103,64) JVS(1060) = -B(119) ! JVS(1061) = Jac_FULL(103,65) JVS(1061) = -B(85) ! JVS(1062) = Jac_FULL(103,66) JVS(1062) = -B(123) ! JVS(1063) = Jac_FULL(103,68) JVS(1063) = -B(95) ! JVS(1064) = Jac_FULL(103,69) JVS(1064) = -B(121) ! JVS(1065) = Jac_FULL(103,70) JVS(1065) = -B(103)-B(105) ! JVS(1066) = Jac_FULL(103,72) JVS(1066) = -B(115)-B(117) ! JVS(1067) = Jac_FULL(103,73) JVS(1067) = -B(99) ! JVS(1068) = Jac_FULL(103,74) JVS(1068) = -B(97) ! JVS(1069) = Jac_FULL(103,75) JVS(1069) = -B(111)-B(113) ! JVS(1070) = Jac_FULL(103,76) JVS(1070) = -B(127) ! JVS(1071) = Jac_FULL(103,78) JVS(1071) = -B(125) ! JVS(1072) = Jac_FULL(103,79) JVS(1072) = -B(107) ! JVS(1073) = Jac_FULL(103,80) JVS(1073) = -B(109) ! JVS(1074) = Jac_FULL(103,83) JVS(1074) = -B(87) ! JVS(1075) = Jac_FULL(103,85) JVS(1075) = -B(91)-B(93) ! JVS(1076) = Jac_FULL(103,86) JVS(1076) = 0 ! JVS(1077) = Jac_FULL(103,87) JVS(1077) = -B(79) ! JVS(1078) = Jac_FULL(103,89) JVS(1078) = 0 ! JVS(1079) = Jac_FULL(103,90) JVS(1079) = -B(149) ! JVS(1080) = Jac_FULL(103,91) JVS(1080) = 0 ! JVS(1081) = Jac_FULL(103,92) JVS(1081) = -B(101) ! JVS(1082) = Jac_FULL(103,93) JVS(1082) = 0 ! JVS(1083) = Jac_FULL(103,94) JVS(1083) = -B(151) ! JVS(1084) = Jac_FULL(103,95) JVS(1084) = -B(26) ! JVS(1085) = Jac_FULL(103,96) JVS(1085) = -B(17) ! JVS(1086) = Jac_FULL(103,97) JVS(1086) = -B(153) ! JVS(1087) = Jac_FULL(103,98) JVS(1087) = -B(1) ! JVS(1088) = Jac_FULL(103,99) JVS(1088) = B(64)+B(499) ! JVS(1089) = Jac_FULL(103,100) JVS(1089) = -B(53)+B(65)+B(508) ! JVS(1090) = Jac_FULL(103,101) JVS(1090) = -B(75) ! JVS(1091) = Jac_FULL(103,102) JVS(1091) = -B(42) ! JVS(1092) = Jac_FULL(103,103) JVS(1092) = -B(2)-B(18)-B(27)-B(43)-B(54)-B(76)-B(80)-B(86)-B(88)-B(92)-B(94)-B(96)-B(98)-B(100)-B(102)-B(104)-B(106)& &-B(108)-B(110)-B(112)-B(114)-B(116)-B(118)-B(120)-B(122)-B(124)-B(126)-B(128)-B(150)-B(152)-B(154)-B(155)& &-B(489)-B(493)-B(497) END SUBROUTINE Jac_SP ! End of Jac_SP function ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! Jac_SP_Vec - function for sparse multiplication: sparse Jacobian times vector ! Arguments : ! JVS - sparse Jacobian of variables ! UV - User vector for variables ! JUV - Jacobian times user vector ! ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ SUBROUTINE Jac_SP_Vec ( JVS, UV, JUV ) ! JVS - sparse Jacobian of variables REAL(kind=dp) :: JVS(LU_NONZERO) ! UV - User vector for variables REAL(kind=dp) :: UV(NVAR) ! JUV - Jacobian times user vector REAL(kind=dp) :: JUV(NVAR) JUV(1) = JVS(1)*UV(1) JUV(2) = JVS(2)*UV(2)+JVS(3)*UV(82) JUV(3) = JVS(4)*UV(3)+JVS(5)*UV(29) JUV(4) = JVS(6)*UV(4)+JVS(7)*UV(67) JUV(5) = JVS(8)*UV(5)+JVS(9)*UV(38) JUV(6) = JVS(10)*UV(6)+JVS(11)*UV(99) JUV(7) = JVS(12)*UV(7)+JVS(13)*UV(98) JUV(8) = JVS(14)*UV(8)+JVS(15)*UV(34) JUV(9) = JVS(16)*UV(9)+JVS(17)*UV(58) JUV(10) = JVS(18)*UV(10)+JVS(19)*UV(27) JUV(11) = JVS(20)*UV(11)+JVS(21)*UV(86) JUV(12) = JVS(22)*UV(12)+JVS(23)*UV(30)+JVS(24)*UV(102) JUV(13) = JVS(25)*UV(13)+JVS(26)*UV(44)+JVS(27)*UV(102) JUV(14) = JVS(28)*UV(14)+JVS(29)*UV(59)+JVS(30)*UV(102) JUV(15) = JVS(31)*UV(15)+JVS(32)*UV(31)+JVS(33)*UV(96) JUV(16) = JVS(34)*UV(16)+JVS(35)*UV(31)+JVS(36)*UV(103) JUV(17) = JVS(37)*UV(17)+JVS(38)*UV(32)+JVS(39)*UV(96) JUV(18) = JVS(40)*UV(18)+JVS(41)*UV(32)+JVS(42)*UV(103) JUV(19) = JVS(43)*UV(19)+JVS(44)*UV(33)+JVS(45)*UV(96) JUV(20) = JVS(46)*UV(20)+JVS(47)*UV(33)+JVS(48)*UV(103) JUV(21) = JVS(49)*UV(21)+JVS(50)*UV(59)+JVS(51)*UV(60)+JVS(52)*UV(91)+JVS(53)*UV(98)+JVS(54)*UV(101)+JVS(55)*UV(102)& &+JVS(56)*UV(103) JUV(22) = JVS(57)*UV(22)+JVS(58)*UV(27)+JVS(59)*UV(29)+JVS(60)*UV(34)+JVS(61)*UV(38)+JVS(62)*UV(58)+JVS(63)*UV(67)& &+JVS(64)*UV(82)+JVS(65)*UV(86)+JVS(66)*UV(98)+JVS(67)*UV(99) JUV(23) = JVS(68)*UV(23)+JVS(69)*UV(102) JUV(24) = JVS(70)*UV(24)+JVS(71)*UV(102) JUV(25) = JVS(72)*UV(25)+JVS(73)*UV(102) JUV(26) = JVS(74)*UV(26)+JVS(75)*UV(102) JUV(27) = JVS(76)*UV(27)+JVS(77)*UV(90)+JVS(78)*UV(99) JUV(28) = JVS(79)*UV(28)+JVS(80)*UV(94)+JVS(81)*UV(99) JUV(29) = JVS(82)*UV(29)+JVS(83)*UV(96)+JVS(84)*UV(102) JUV(30) = JVS(85)*UV(30)+JVS(86)*UV(44)+JVS(87)*UV(100)+JVS(88)*UV(102) JUV(31) = JVS(89)*UV(25)+JVS(90)*UV(31)+JVS(91)*UV(96)+JVS(92)*UV(102)+JVS(93)*UV(103) JUV(32) = JVS(94)*UV(23)+JVS(95)*UV(32)+JVS(96)*UV(96)+JVS(97)*UV(102)+JVS(98)*UV(103) JUV(33) = JVS(99)*UV(24)+JVS(100)*UV(33)+JVS(101)*UV(96)+JVS(102)*UV(102)+JVS(103)*UV(103) JUV(34) = JVS(104)*UV(34)+JVS(105)*UV(99)+JVS(106)*UV(101) JUV(35) = JVS(107)*UV(35)+JVS(108)*UV(100)+JVS(109)*UV(102) JUV(36) = JVS(110)*UV(36)+JVS(111)*UV(100)+JVS(112)*UV(102) JUV(37) = JVS(113)*UV(37)+JVS(114)*UV(99)+JVS(115)*UV(102)+JVS(116)*UV(103) JUV(38) = JVS(117)*UV(38)+JVS(118)*UV(99)+JVS(119)*UV(100) JUV(39) = JVS(120)*UV(39)+JVS(121)*UV(95)+JVS(122)*UV(96)+JVS(123)*UV(102) JUV(40) = JVS(124)*UV(40)+JVS(125)*UV(85)+JVS(126)*UV(96)+JVS(127)*UV(102) JUV(41) = JVS(128)*UV(41)+JVS(129)*UV(65)+JVS(130)*UV(96)+JVS(131)*UV(102) JUV(42) = JVS(132)*UV(42)+JVS(133)*UV(66)+JVS(134)*UV(96)+JVS(135)*UV(102) JUV(43) = JVS(136)*UV(43)+JVS(137)*UV(90)+JVS(138)*UV(96)+JVS(139)*UV(102) JUV(44) = JVS(140)*UV(44)+JVS(141)*UV(100)+JVS(142)*UV(102) JUV(45) = JVS(143)*UV(45)+JVS(144)*UV(96)+JVS(145)*UV(97)+JVS(146)*UV(102) JUV(46) = JVS(147)*UV(46)+JVS(148)*UV(87)+JVS(149)*UV(96)+JVS(150)*UV(102) JUV(47) = JVS(151)*UV(47)+JVS(152)*UV(94)+JVS(153)*UV(96)+JVS(154)*UV(102) JUV(48) = JVS(155)*UV(48)+JVS(156)*UV(96)+JVS(157)*UV(101)+JVS(158)*UV(102) JUV(49) = JVS(159)*UV(49)+JVS(160)*UV(96)+JVS(161)*UV(99)+JVS(162)*UV(102) JUV(50) = JVS(163)*UV(50)+JVS(164)*UV(83)+JVS(165)*UV(96)+JVS(166)*UV(102) JUV(51) = JVS(167)*UV(51)+JVS(168)*UV(76)+JVS(169)*UV(96)+JVS(170)*UV(102) JUV(52) = JVS(171)*UV(52)+JVS(172)*UV(78)+JVS(173)*UV(96)+JVS(174)*UV(102) JUV(53) = JVS(175)*UV(53)+JVS(176)*UV(72)+JVS(177)*UV(96)+JVS(178)*UV(102) JUV(54) = JVS(179)*UV(54)+JVS(180)*UV(79)+JVS(181)*UV(96)+JVS(182)*UV(102) JUV(55) = JVS(183)*UV(55)+JVS(184)*UV(75)+JVS(185)*UV(96)+JVS(186)*UV(102) JUV(56) = JVS(187)*UV(56)+JVS(188)*UV(64)+JVS(189)*UV(69)+JVS(190)*UV(80)+JVS(191)*UV(96)+JVS(192)*UV(102) JUV(57) = JVS(193)*UV(57)+JVS(194)*UV(70)+JVS(195)*UV(92)+JVS(196)*UV(96)+JVS(197)*UV(102) JUV(58) = JVS(198)*UV(58)+JVS(199)*UV(97)+JVS(200)*UV(98)+JVS(201)*UV(99)+JVS(202)*UV(102) JUV(59) = JVS(203)*UV(59)+JVS(204)*UV(98)+JVS(205)*UV(100)+JVS(206)*UV(102) JUV(60) = JVS(207)*UV(53)+JVS(208)*UV(54)+JVS(209)*UV(59)+JVS(210)*UV(60)+JVS(211)*UV(61)+JVS(212)*UV(62)+JVS(213)& &*UV(63)+JVS(214)*UV(71)+JVS(215)*UV(72)+JVS(216)*UV(79)+JVS(217)*UV(81)+JVS(218)*UV(82)+JVS(219)*UV(84)& &+JVS(220)*UV(88)+JVS(221)*UV(91)+JVS(222)*UV(93)+JVS(223)*UV(95)+JVS(225)*UV(98)+JVS(226)*UV(100)+JVS(227)& &*UV(101)+JVS(228)*UV(102)+JVS(229)*UV(103) JUV(61) = JVS(230)*UV(59)+JVS(231)*UV(61)+JVS(232)*UV(93)+JVS(233)*UV(98)+JVS(234)*UV(100)+JVS(235)*UV(102) JUV(62) = JVS(236)*UV(42)+JVS(237)*UV(62)+JVS(238)*UV(66)+JVS(239)*UV(85)+JVS(240)*UV(86)+JVS(241)*UV(95)+JVS(243)& &*UV(101)+JVS(244)*UV(102)+JVS(245)*UV(103) JUV(63) = JVS(246)*UV(54)+JVS(247)*UV(55)+JVS(248)*UV(63)+JVS(249)*UV(64)+JVS(250)*UV(71)+JVS(251)*UV(75)+JVS(252)& &*UV(79)+JVS(253)*UV(80)+JVS(254)*UV(95)+JVS(256)*UV(98)+JVS(257)*UV(101)+JVS(258)*UV(102)+JVS(259)*UV(103) JUV(64) = JVS(260)*UV(64)+JVS(261)*UV(95)+JVS(262)*UV(96)+JVS(263)*UV(101)+JVS(264)*UV(103) JUV(65) = JVS(265)*UV(26)+JVS(266)*UV(41)+JVS(267)*UV(65)+JVS(268)*UV(85)+JVS(269)*UV(86)+JVS(270)*UV(95)+JVS(271)& &*UV(96)+JVS(272)*UV(101)+JVS(273)*UV(102)+JVS(274)*UV(103) JUV(66) = JVS(275)*UV(26)+JVS(276)*UV(42)+JVS(277)*UV(66)+JVS(278)*UV(85)+JVS(279)*UV(86)+JVS(280)*UV(95)+JVS(281)& &*UV(96)+JVS(282)*UV(101)+JVS(283)*UV(102)+JVS(284)*UV(103) JUV(67) = JVS(285)*UV(35)+JVS(286)*UV(36)+JVS(287)*UV(38)+JVS(288)*UV(44)+JVS(289)*UV(64)+JVS(290)*UV(67)+JVS(291)& &*UV(70)+JVS(292)*UV(72)+JVS(293)*UV(75)+JVS(294)*UV(78)+JVS(295)*UV(79)+JVS(296)*UV(80)+JVS(297)*UV(81)& &+JVS(298)*UV(82)+JVS(299)*UV(84)+JVS(300)*UV(88)+JVS(301)*UV(89)+JVS(302)*UV(91)+JVS(303)*UV(95)+JVS(305)& &*UV(99)+JVS(306)*UV(100)+JVS(307)*UV(101)+JVS(308)*UV(102)+JVS(309)*UV(103) JUV(68) = JVS(310)*UV(68)+JVS(311)*UV(86)+JVS(312)*UV(95)+JVS(313)*UV(96)+JVS(314)*UV(101)+JVS(315)*UV(102)+JVS(316)& &*UV(103) JUV(69) = JVS(317)*UV(69)+JVS(318)*UV(91)+JVS(319)*UV(95)+JVS(320)*UV(96)+JVS(321)*UV(100)+JVS(322)*UV(101)+JVS(323)& &*UV(103) JUV(70) = JVS(324)*UV(70)+JVS(325)*UV(92)+JVS(326)*UV(95)+JVS(327)*UV(96)+JVS(328)*UV(101)+JVS(329)*UV(103) JUV(71) = JVS(330)*UV(57)+JVS(331)*UV(70)+JVS(332)*UV(71)+JVS(333)*UV(92)+JVS(334)*UV(95)+JVS(336)*UV(98)+JVS(337)& &*UV(101)+JVS(338)*UV(102)+JVS(339)*UV(103) JUV(72) = JVS(340)*UV(53)+JVS(341)*UV(72)+JVS(342)*UV(91)+JVS(343)*UV(95)+JVS(344)*UV(96)+JVS(345)*UV(101)+JVS(346)& &*UV(102)+JVS(347)*UV(103) JUV(73) = JVS(348)*UV(73)+JVS(349)*UV(89)+JVS(350)*UV(95)+JVS(351)*UV(96)+JVS(352)*UV(100)+JVS(353)*UV(101)+JVS(354)& &*UV(102)+JVS(355)*UV(103) JUV(74) = JVS(356)*UV(62)+JVS(358)*UV(74)+JVS(361)*UV(95)+JVS(362)*UV(96)+JVS(363)*UV(101)+JVS(364)*UV(102)+JVS(365)& &*UV(103) JUV(75) = JVS(366)*UV(55)+JVS(367)*UV(75)+JVS(368)*UV(93)+JVS(369)*UV(95)+JVS(370)*UV(96)+JVS(371)*UV(101)+JVS(372)& &*UV(102)+JVS(373)*UV(103) JUV(76) = JVS(374)*UV(51)+JVS(375)*UV(61)+JVS(376)*UV(76)+JVS(378)*UV(95)+JVS(379)*UV(96)+JVS(381)*UV(100)+JVS(382)& &*UV(101)+JVS(383)*UV(102)+JVS(384)*UV(103) JUV(77) = JVS(385)*UV(53)+JVS(386)*UV(54)+JVS(387)*UV(58)+JVS(388)*UV(71)+JVS(389)*UV(72)+JVS(390)*UV(74)+JVS(391)& &*UV(77)+JVS(392)*UV(79)+JVS(393)*UV(80)+JVS(394)*UV(83)+JVS(399)*UV(95)+JVS(402)*UV(98)+JVS(404)*UV(101)& &+JVS(405)*UV(102)+JVS(406)*UV(103) JUV(78) = JVS(407)*UV(52)+JVS(408)*UV(59)+JVS(409)*UV(78)+JVS(410)*UV(95)+JVS(411)*UV(96)+JVS(413)*UV(100)+JVS(414)& &*UV(101)+JVS(415)*UV(102)+JVS(416)*UV(103) JUV(79) = JVS(417)*UV(54)+JVS(418)*UV(71)+JVS(419)*UV(79)+JVS(421)*UV(95)+JVS(422)*UV(96)+JVS(424)*UV(101)+JVS(425)& &*UV(102)+JVS(426)*UV(103) JUV(80) = JVS(427)*UV(56)+JVS(430)*UV(80)+JVS(432)*UV(95)+JVS(433)*UV(96)+JVS(435)*UV(101)+JVS(436)*UV(102)+JVS(437)& &*UV(103) JUV(81) = JVS(438)*UV(40)+JVS(439)*UV(41)+JVS(440)*UV(51)+JVS(441)*UV(52)+JVS(442)*UV(55)+JVS(443)*UV(56)+JVS(444)& &*UV(64)+JVS(445)*UV(65)+JVS(446)*UV(68)+JVS(447)*UV(69)+JVS(449)*UV(76)+JVS(450)*UV(78)+JVS(451)*UV(80)& &+JVS(452)*UV(81)+JVS(453)*UV(83)+JVS(454)*UV(85)+JVS(455)*UV(86)+JVS(458)*UV(95)+JVS(461)*UV(100)+JVS(462)& &*UV(101)+JVS(463)*UV(102)+JVS(464)*UV(103) JUV(82) = JVS(465)*UV(1)+JVS(466)*UV(39)+JVS(467)*UV(44)+JVS(468)*UV(45)+JVS(469)*UV(47)+JVS(470)*UV(48)+JVS(471)& &*UV(50)+JVS(472)*UV(53)+JVS(473)*UV(55)+JVS(474)*UV(57)+JVS(475)*UV(58)+JVS(476)*UV(59)+JVS(477)*UV(61)& &+JVS(478)*UV(63)+JVS(479)*UV(64)+JVS(480)*UV(65)+JVS(481)*UV(66)+JVS(482)*UV(68)+JVS(483)*UV(69)+JVS(484)& &*UV(70)+JVS(485)*UV(71)+JVS(486)*UV(72)+JVS(487)*UV(73)+JVS(488)*UV(74)+JVS(489)*UV(75)+JVS(490)*UV(76)& &+JVS(491)*UV(77)+JVS(492)*UV(78)+JVS(493)*UV(79)+JVS(494)*UV(80)+JVS(495)*UV(82)+JVS(496)*UV(83)+JVS(497)& &*UV(84)+JVS(498)*UV(85)+JVS(500)*UV(87)+JVS(502)*UV(90)+JVS(503)*UV(91)+JVS(504)*UV(92)+JVS(505)*UV(93)& &+JVS(506)*UV(94)+JVS(507)*UV(95)+JVS(508)*UV(96)+JVS(509)*UV(97)+JVS(510)*UV(98)+JVS(512)*UV(100)+JVS(513)& &*UV(101)+JVS(514)*UV(102)+JVS(515)*UV(103) JUV(83) = JVS(516)*UV(50)+JVS(517)*UV(61)+JVS(518)*UV(83)+JVS(520)*UV(95)+JVS(521)*UV(96)+JVS(524)*UV(101)+JVS(525)& &*UV(102)+JVS(526)*UV(103) JUV(84) = JVS(527)*UV(43)+JVS(528)*UV(46)+JVS(529)*UV(50)+JVS(530)*UV(61)+JVS(531)*UV(68)+JVS(532)*UV(73)+JVS(533)& &*UV(76)+JVS(534)*UV(83)+JVS(535)*UV(84)+JVS(536)*UV(85)+JVS(537)*UV(86)+JVS(538)*UV(87)+JVS(539)*UV(88)& &+JVS(542)*UV(93)+JVS(543)*UV(95)+JVS(545)*UV(98)+JVS(546)*UV(100)+JVS(547)*UV(101)+JVS(548)*UV(102)+JVS(549)& &*UV(103) JUV(85) = JVS(550)*UV(35)+JVS(551)*UV(40)+JVS(552)*UV(68)+JVS(553)*UV(85)+JVS(555)*UV(95)+JVS(556)*UV(96)+JVS(557)& &*UV(100)+JVS(558)*UV(101)+JVS(559)*UV(102)+JVS(560)*UV(103) JUV(86) = JVS(561)*UV(68)+JVS(562)*UV(73)+JVS(563)*UV(74)+JVS(564)*UV(85)+JVS(565)*UV(86)+JVS(568)*UV(96)+JVS(571)& &*UV(102)+JVS(572)*UV(103) JUV(87) = JVS(573)*UV(36)+JVS(574)*UV(46)+JVS(575)*UV(81)+JVS(577)*UV(85)+JVS(578)*UV(86)+JVS(579)*UV(87)+JVS(580)& &*UV(89)+JVS(581)*UV(90)+JVS(584)*UV(95)+JVS(585)*UV(96)+JVS(587)*UV(100)+JVS(588)*UV(101)+JVS(589)*UV(102)& &+JVS(590)*UV(103) JUV(88) = JVS(591)*UV(54)+JVS(592)*UV(55)+JVS(593)*UV(64)+JVS(594)*UV(69)+JVS(595)*UV(71)+JVS(596)*UV(72)+JVS(597)& &*UV(73)+JVS(598)*UV(74)+JVS(599)*UV(75)+JVS(600)*UV(77)+JVS(601)*UV(79)+JVS(606)*UV(88)+JVS(608)*UV(91)& &+JVS(610)*UV(93)+JVS(611)*UV(95)+JVS(612)*UV(96)+JVS(614)*UV(98)+JVS(616)*UV(100)+JVS(617)*UV(101)+JVS(618)& &*UV(102)+JVS(619)*UV(103) JUV(89) = JVS(620)*UV(70)+JVS(621)*UV(72)+JVS(622)*UV(73)+JVS(623)*UV(74)+JVS(624)*UV(75)+JVS(625)*UV(79)+JVS(626)& &*UV(85)+JVS(627)*UV(86)+JVS(628)*UV(89)+JVS(630)*UV(92)+JVS(632)*UV(95)+JVS(635)*UV(100)+JVS(636)*UV(101)& &+JVS(637)*UV(102)+JVS(638)*UV(103) JUV(90) = JVS(639)*UV(27)+JVS(640)*UV(43)+JVS(641)*UV(81)+JVS(645)*UV(89)+JVS(646)*UV(90)+JVS(650)*UV(95)+JVS(651)& &*UV(96)+JVS(653)*UV(99)+JVS(654)*UV(100)+JVS(655)*UV(101)+JVS(656)*UV(102)+JVS(657)*UV(103) JUV(91) = JVS(658)*UV(57)+JVS(659)*UV(59)+JVS(661)*UV(78)+JVS(662)*UV(91)+JVS(663)*UV(92)+JVS(664)*UV(95)+JVS(666)& &*UV(98)+JVS(667)*UV(100)+JVS(668)*UV(101)+JVS(669)*UV(102)+JVS(670)*UV(103) JUV(92) = JVS(671)*UV(57)+JVS(672)*UV(59)+JVS(674)*UV(92)+JVS(675)*UV(95)+JVS(676)*UV(96)+JVS(679)*UV(101)+JVS(680)& &*UV(102)+JVS(681)*UV(103) JUV(93) = JVS(682)*UV(57)+JVS(683)*UV(59)+JVS(685)*UV(78)+JVS(686)*UV(92)+JVS(687)*UV(93)+JVS(688)*UV(95)+JVS(690)& &*UV(98)+JVS(692)*UV(101)+JVS(693)*UV(102)+JVS(694)*UV(103) JUV(94) = JVS(695)*UV(28)+JVS(696)*UV(47)+JVS(697)*UV(63)+JVS(706)*UV(94)+JVS(707)*UV(95)+JVS(708)*UV(96)+JVS(710)& &*UV(99)+JVS(712)*UV(101)+JVS(713)*UV(102)+JVS(714)*UV(103) JUV(95) = JVS(715)*UV(34)+JVS(716)*UV(39)+JVS(717)*UV(44)+JVS(718)*UV(48)+JVS(719)*UV(61)+JVS(720)*UV(62)+JVS(721)& &*UV(64)+JVS(722)*UV(65)+JVS(723)*UV(66)+JVS(724)*UV(68)+JVS(725)*UV(69)+JVS(726)*UV(70)+JVS(727)*UV(72)& &+JVS(728)*UV(73)+JVS(729)*UV(74)+JVS(730)*UV(75)+JVS(731)*UV(76)+JVS(732)*UV(78)+JVS(733)*UV(79)+JVS(734)& &*UV(80)+JVS(735)*UV(83)+JVS(736)*UV(84)+JVS(737)*UV(85)+JVS(738)*UV(86)+JVS(739)*UV(87)+JVS(741)*UV(89)& &+JVS(742)*UV(90)+JVS(744)*UV(92)+JVS(745)*UV(93)+JVS(746)*UV(94)+JVS(747)*UV(95)+JVS(748)*UV(96)+JVS(749)& &*UV(97)+JVS(750)*UV(98)+JVS(752)*UV(100)+JVS(753)*UV(101)+JVS(754)*UV(102)+JVS(755)*UV(103) JUV(96) = JVS(756)*UV(29)+JVS(757)*UV(30)+JVS(758)*UV(31)+JVS(759)*UV(32)+JVS(760)*UV(33)+JVS(761)*UV(39)+JVS(762)& &*UV(40)+JVS(763)*UV(41)+JVS(764)*UV(42)+JVS(765)*UV(43)+JVS(767)*UV(46)+JVS(768)*UV(47)+JVS(769)*UV(49)& &+JVS(770)*UV(50)+JVS(771)*UV(51)+JVS(772)*UV(52)+JVS(773)*UV(53)+JVS(774)*UV(54)+JVS(775)*UV(55)+JVS(776)& &*UV(56)+JVS(777)*UV(57)+JVS(778)*UV(58)+JVS(779)*UV(59)+JVS(780)*UV(60)+JVS(781)*UV(61)+JVS(783)*UV(63)& &+JVS(784)*UV(64)+JVS(785)*UV(65)+JVS(786)*UV(66)+JVS(787)*UV(68)+JVS(788)*UV(69)+JVS(789)*UV(70)+JVS(790)& &*UV(71)+JVS(791)*UV(72)+JVS(792)*UV(73)+JVS(793)*UV(74)+JVS(794)*UV(75)+JVS(795)*UV(76)+JVS(796)*UV(77)& &+JVS(797)*UV(78)+JVS(798)*UV(79)+JVS(799)*UV(80)+JVS(800)*UV(81)+JVS(801)*UV(82)+JVS(802)*UV(83)+JVS(803)& &*UV(84)+JVS(804)*UV(85)+JVS(805)*UV(86)+JVS(806)*UV(87)+JVS(807)*UV(88)+JVS(809)*UV(90)+JVS(810)*UV(91)& &+JVS(811)*UV(92)+JVS(812)*UV(93)+JVS(813)*UV(94)+JVS(814)*UV(95)+JVS(815)*UV(96)+JVS(816)*UV(97)+JVS(817)& &*UV(98)+JVS(818)*UV(99)+JVS(819)*UV(100)+JVS(820)*UV(101)+JVS(821)*UV(102)+JVS(822)*UV(103) JUV(97) = JVS(823)*UV(45)+JVS(824)*UV(58)+JVS(825)*UV(71)+JVS(826)*UV(91)+JVS(828)*UV(93)+JVS(829)*UV(95)+JVS(830)& &*UV(96)+JVS(831)*UV(97)+JVS(833)*UV(99)+JVS(834)*UV(100)+JVS(835)*UV(101)+JVS(836)*UV(102)+JVS(837)*UV(103) JUV(98) = JVS(838)*UV(38)+JVS(839)*UV(58)+JVS(840)*UV(59)+JVS(841)*UV(61)+JVS(842)*UV(71)+JVS(843)*UV(90)+JVS(844)& &*UV(91)+JVS(846)*UV(93)+JVS(847)*UV(94)+JVS(848)*UV(95)+JVS(849)*UV(96)+JVS(850)*UV(97)+JVS(851)*UV(98)& &+JVS(852)*UV(99)+JVS(853)*UV(100)+JVS(854)*UV(101)+JVS(855)*UV(102)+JVS(856)*UV(103) JUV(99) = JVS(857)*UV(27)+JVS(858)*UV(28)+JVS(859)*UV(34)+JVS(860)*UV(37)+JVS(861)*UV(38)+JVS(862)*UV(49)+JVS(863)& &*UV(51)+JVS(864)*UV(52)+JVS(865)*UV(56)+JVS(866)*UV(58)+JVS(867)*UV(64)+JVS(868)*UV(65)+JVS(869)*UV(66)& &+JVS(870)*UV(67)+JVS(871)*UV(68)+JVS(872)*UV(69)+JVS(873)*UV(70)+JVS(874)*UV(72)+JVS(875)*UV(73)+JVS(876)& &*UV(74)+JVS(877)*UV(75)+JVS(878)*UV(76)+JVS(879)*UV(78)+JVS(880)*UV(79)+JVS(881)*UV(80)+JVS(884)*UV(83)& &+JVS(886)*UV(85)+JVS(887)*UV(86)+JVS(888)*UV(87)+JVS(891)*UV(90)+JVS(893)*UV(92)+JVS(895)*UV(94)+JVS(896)& &*UV(95)+JVS(897)*UV(96)+JVS(898)*UV(97)+JVS(899)*UV(98)+JVS(900)*UV(99)+JVS(901)*UV(100)+JVS(902)*UV(101)& &+JVS(903)*UV(102)+JVS(904)*UV(103) JUV(100) = JVS(905)*UV(34)+JVS(906)*UV(35)+JVS(907)*UV(36)+JVS(908)*UV(38)+JVS(909)*UV(44)+JVS(910)*UV(49)+JVS(911)& &*UV(59)+JVS(912)*UV(61)+JVS(913)*UV(67)+JVS(920)*UV(81)+JVS(921)*UV(82)+JVS(923)*UV(84)+JVS(927)*UV(88)& &+JVS(928)*UV(89)+JVS(930)*UV(91)+JVS(935)*UV(96)+JVS(937)*UV(98)+JVS(938)*UV(99)+JVS(939)*UV(100)+JVS(941)& &*UV(102)+JVS(942)*UV(103) JUV(101) = JVS(943)*UV(34)+JVS(944)*UV(45)+JVS(945)*UV(48)+JVS(946)*UV(55)+JVS(947)*UV(62)+JVS(948)*UV(64)+JVS(949)& &*UV(65)+JVS(950)*UV(66)+JVS(951)*UV(68)+JVS(952)*UV(69)+JVS(953)*UV(70)+JVS(954)*UV(72)+JVS(955)*UV(73)& &+JVS(956)*UV(74)+JVS(957)*UV(75)+JVS(958)*UV(76)+JVS(959)*UV(77)+JVS(960)*UV(78)+JVS(961)*UV(79)+JVS(962)& &*UV(80)+JVS(963)*UV(83)+JVS(964)*UV(84)+JVS(965)*UV(85)+JVS(967)*UV(87)+JVS(968)*UV(88)+JVS(969)*UV(89)& &+JVS(970)*UV(90)+JVS(971)*UV(91)+JVS(972)*UV(92)+JVS(973)*UV(93)+JVS(974)*UV(94)+JVS(975)*UV(95)+JVS(976)& &*UV(96)+JVS(977)*UV(97)+JVS(979)*UV(99)+JVS(980)*UV(100)+JVS(981)*UV(101)+JVS(982)*UV(102)+JVS(983)*UV(103) JUV(102) = JVS(984)*UV(23)+JVS(985)*UV(24)+JVS(986)*UV(25)+JVS(987)*UV(26)+JVS(988)*UV(29)+JVS(989)*UV(30)+JVS(990)& &*UV(35)+JVS(991)*UV(36)+JVS(992)*UV(37)+JVS(993)*UV(39)+JVS(994)*UV(40)+JVS(995)*UV(41)+JVS(996)*UV(42)& &+JVS(997)*UV(43)+JVS(998)*UV(44)+JVS(999)*UV(45)+JVS(1000)*UV(46)+JVS(1001)*UV(47)+JVS(1002)*UV(48)& &+JVS(1003)*UV(49)+JVS(1004)*UV(50)+JVS(1005)*UV(51)+JVS(1006)*UV(52)+JVS(1007)*UV(53)+JVS(1008)*UV(54)& &+JVS(1009)*UV(55)+JVS(1010)*UV(56)+JVS(1011)*UV(57)+JVS(1012)*UV(58)+JVS(1013)*UV(59)+JVS(1014)*UV(60)& &+JVS(1015)*UV(61)+JVS(1016)*UV(62)+JVS(1017)*UV(63)+JVS(1021)*UV(67)+JVS(1024)*UV(71)+JVS(1028)*UV(77)& &+JVS(1032)*UV(81)+JVS(1033)*UV(82)+JVS(1035)*UV(84)+JVS(1037)*UV(86)+JVS(1039)*UV(88)+JVS(1040)*UV(89)& &+JVS(1042)*UV(91)+JVS(1044)*UV(93)+JVS(1047)*UV(96)+JVS(1049)*UV(98)+JVS(1050)*UV(99)+JVS(1051)*UV(100)& &+JVS(1052)*UV(101)+JVS(1053)*UV(102)+JVS(1054)*UV(103) JUV(103) = JVS(1055)*UV(31)+JVS(1056)*UV(32)+JVS(1057)*UV(33)+JVS(1058)*UV(37)+JVS(1059)*UV(38)+JVS(1060)*UV(64)& &+JVS(1061)*UV(65)+JVS(1062)*UV(66)+JVS(1063)*UV(68)+JVS(1064)*UV(69)+JVS(1065)*UV(70)+JVS(1066)*UV(72)& &+JVS(1067)*UV(73)+JVS(1068)*UV(74)+JVS(1069)*UV(75)+JVS(1070)*UV(76)+JVS(1071)*UV(78)+JVS(1072)*UV(79)& &+JVS(1073)*UV(80)+JVS(1074)*UV(83)+JVS(1075)*UV(85)+JVS(1077)*UV(87)+JVS(1079)*UV(90)+JVS(1081)*UV(92)& &+JVS(1083)*UV(94)+JVS(1084)*UV(95)+JVS(1085)*UV(96)+JVS(1086)*UV(97)+JVS(1087)*UV(98)+JVS(1088)*UV(99)& &+JVS(1089)*UV(100)+JVS(1090)*UV(101)+JVS(1091)*UV(102)+JVS(1092)*UV(103) END SUBROUTINE Jac_SP_Vec ! End of Jac_SP_Vec function ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! JacTR_SP_Vec - sparse multiplication: sparse Jacobian transposed times vector ! Arguments : ! JVS - sparse Jacobian of variables ! UV - User vector for variables ! JTUV - Jacobian transposed times user vector ! ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ SUBROUTINE JacTR_SP_Vec ( JVS, UV, JTUV ) ! JVS - sparse Jacobian of variables REAL(kind=dp) :: JVS(LU_NONZERO) ! UV - User vector for variables REAL(kind=dp) :: UV(NVAR) ! JTUV - Jacobian transposed times user vector REAL(kind=dp) :: JTUV(NVAR) JTUV(1) = JVS(1)*UV(1)+JVS(465)*UV(82) JTUV(2) = JVS(2)*UV(2) JTUV(3) = JVS(4)*UV(3) JTUV(4) = JVS(6)*UV(4) JTUV(5) = JVS(8)*UV(5) JTUV(6) = JVS(10)*UV(6) JTUV(7) = JVS(12)*UV(7) JTUV(8) = JVS(14)*UV(8) JTUV(9) = JVS(16)*UV(9) JTUV(10) = JVS(18)*UV(10) JTUV(11) = JVS(20)*UV(11) JTUV(12) = JVS(22)*UV(12) JTUV(13) = JVS(25)*UV(13) JTUV(14) = JVS(28)*UV(14) JTUV(15) = JVS(31)*UV(15) JTUV(16) = JVS(34)*UV(16) JTUV(17) = JVS(37)*UV(17) JTUV(18) = JVS(40)*UV(18) JTUV(19) = JVS(43)*UV(19) JTUV(20) = JVS(46)*UV(20) JTUV(21) = JVS(49)*UV(21) JTUV(22) = JVS(57)*UV(22) JTUV(23) = JVS(68)*UV(23)+JVS(94)*UV(32)+JVS(984)*UV(102) JTUV(24) = JVS(70)*UV(24)+JVS(99)*UV(33)+JVS(985)*UV(102) JTUV(25) = JVS(72)*UV(25)+JVS(89)*UV(31)+JVS(986)*UV(102) JTUV(26) = JVS(74)*UV(26)+JVS(265)*UV(65)+JVS(275)*UV(66)+JVS(987)*UV(102) JTUV(27) = JVS(19)*UV(10)+JVS(58)*UV(22)+JVS(76)*UV(27)+JVS(639)*UV(90)+JVS(857)*UV(99) JTUV(28) = JVS(79)*UV(28)+JVS(695)*UV(94)+JVS(858)*UV(99) JTUV(29) = JVS(5)*UV(3)+JVS(59)*UV(22)+JVS(82)*UV(29)+JVS(756)*UV(96)+JVS(988)*UV(102) JTUV(30) = JVS(23)*UV(12)+JVS(85)*UV(30)+JVS(757)*UV(96)+JVS(989)*UV(102) JTUV(31) = JVS(32)*UV(15)+JVS(35)*UV(16)+JVS(90)*UV(31)+JVS(758)*UV(96)+JVS(1055)*UV(103) JTUV(32) = JVS(38)*UV(17)+JVS(41)*UV(18)+JVS(95)*UV(32)+JVS(759)*UV(96)+JVS(1056)*UV(103) JTUV(33) = JVS(44)*UV(19)+JVS(47)*UV(20)+JVS(100)*UV(33)+JVS(760)*UV(96)+JVS(1057)*UV(103) JTUV(34) = JVS(15)*UV(8)+JVS(60)*UV(22)+JVS(104)*UV(34)+JVS(715)*UV(95)+JVS(859)*UV(99)+JVS(905)*UV(100)+JVS(943)& &*UV(101) JTUV(35) = JVS(107)*UV(35)+JVS(285)*UV(67)+JVS(550)*UV(85)+JVS(906)*UV(100)+JVS(990)*UV(102) JTUV(36) = JVS(110)*UV(36)+JVS(286)*UV(67)+JVS(573)*UV(87)+JVS(907)*UV(100)+JVS(991)*UV(102) JTUV(37) = JVS(113)*UV(37)+JVS(860)*UV(99)+JVS(992)*UV(102)+JVS(1058)*UV(103) JTUV(38) = JVS(9)*UV(5)+JVS(61)*UV(22)+JVS(117)*UV(38)+JVS(287)*UV(67)+JVS(838)*UV(98)+JVS(861)*UV(99)+JVS(908)& &*UV(100)+JVS(1059)*UV(103) JTUV(39) = JVS(120)*UV(39)+JVS(466)*UV(82)+JVS(716)*UV(95)+JVS(761)*UV(96)+JVS(993)*UV(102) JTUV(40) = JVS(124)*UV(40)+JVS(438)*UV(81)+JVS(551)*UV(85)+JVS(762)*UV(96)+JVS(994)*UV(102) JTUV(41) = JVS(128)*UV(41)+JVS(266)*UV(65)+JVS(439)*UV(81)+JVS(763)*UV(96)+JVS(995)*UV(102) JTUV(42) = JVS(132)*UV(42)+JVS(236)*UV(62)+JVS(276)*UV(66)+JVS(764)*UV(96)+JVS(996)*UV(102) JTUV(43) = JVS(136)*UV(43)+JVS(527)*UV(84)+JVS(640)*UV(90)+JVS(765)*UV(96)+JVS(997)*UV(102) JTUV(44) = JVS(26)*UV(13)+JVS(86)*UV(30)+JVS(140)*UV(44)+JVS(288)*UV(67)+JVS(467)*UV(82)+JVS(717)*UV(95)+JVS(909)& &*UV(100)+JVS(998)*UV(102) JTUV(45) = JVS(143)*UV(45)+JVS(468)*UV(82)+JVS(823)*UV(97)+JVS(944)*UV(101)+JVS(999)*UV(102) JTUV(46) = JVS(147)*UV(46)+JVS(528)*UV(84)+JVS(574)*UV(87)+JVS(767)*UV(96)+JVS(1000)*UV(102) JTUV(47) = JVS(151)*UV(47)+JVS(469)*UV(82)+JVS(696)*UV(94)+JVS(768)*UV(96)+JVS(1001)*UV(102) JTUV(48) = JVS(155)*UV(48)+JVS(470)*UV(82)+JVS(718)*UV(95)+JVS(945)*UV(101)+JVS(1002)*UV(102) JTUV(49) = JVS(159)*UV(49)+JVS(769)*UV(96)+JVS(862)*UV(99)+JVS(910)*UV(100)+JVS(1003)*UV(102) JTUV(50) = JVS(163)*UV(50)+JVS(471)*UV(82)+JVS(516)*UV(83)+JVS(529)*UV(84)+JVS(770)*UV(96)+JVS(1004)*UV(102) JTUV(51) = JVS(167)*UV(51)+JVS(374)*UV(76)+JVS(440)*UV(81)+JVS(771)*UV(96)+JVS(863)*UV(99)+JVS(1005)*UV(102) JTUV(52) = JVS(171)*UV(52)+JVS(407)*UV(78)+JVS(441)*UV(81)+JVS(772)*UV(96)+JVS(864)*UV(99)+JVS(1006)*UV(102) JTUV(53) = JVS(175)*UV(53)+JVS(207)*UV(60)+JVS(340)*UV(72)+JVS(385)*UV(77)+JVS(472)*UV(82)+JVS(773)*UV(96)+JVS(1007)& &*UV(102) JTUV(54) = JVS(179)*UV(54)+JVS(208)*UV(60)+JVS(246)*UV(63)+JVS(386)*UV(77)+JVS(417)*UV(79)+JVS(591)*UV(88)+JVS(774)& &*UV(96)+JVS(1008)*UV(102) JTUV(55) = JVS(183)*UV(55)+JVS(247)*UV(63)+JVS(366)*UV(75)+JVS(442)*UV(81)+JVS(473)*UV(82)+JVS(592)*UV(88)+JVS(775)& &*UV(96)+JVS(946)*UV(101)+JVS(1009)*UV(102) JTUV(56) = JVS(187)*UV(56)+JVS(427)*UV(80)+JVS(443)*UV(81)+JVS(776)*UV(96)+JVS(865)*UV(99)+JVS(1010)*UV(102) JTUV(57) = JVS(193)*UV(57)+JVS(330)*UV(71)+JVS(474)*UV(82)+JVS(658)*UV(91)+JVS(671)*UV(92)+JVS(682)*UV(93)+JVS(777)& &*UV(96)+JVS(1011)*UV(102) JTUV(58) = JVS(17)*UV(9)+JVS(62)*UV(22)+JVS(198)*UV(58)+JVS(387)*UV(77)+JVS(475)*UV(82)+JVS(778)*UV(96)+JVS(824)& &*UV(97)+JVS(839)*UV(98)+JVS(866)*UV(99)+JVS(1012)*UV(102) JTUV(59) = JVS(29)*UV(14)+JVS(50)*UV(21)+JVS(203)*UV(59)+JVS(209)*UV(60)+JVS(230)*UV(61)+JVS(408)*UV(78)+JVS(476)& &*UV(82)+JVS(659)*UV(91)+JVS(672)*UV(92)+JVS(683)*UV(93)+JVS(779)*UV(96)+JVS(840)*UV(98)+JVS(911)*UV(100)& &+JVS(1013)*UV(102) JTUV(60) = JVS(51)*UV(21)+JVS(210)*UV(60)+JVS(780)*UV(96)+JVS(1014)*UV(102) JTUV(61) = JVS(211)*UV(60)+JVS(231)*UV(61)+JVS(375)*UV(76)+JVS(477)*UV(82)+JVS(517)*UV(83)+JVS(530)*UV(84)+JVS(719)& &*UV(95)+JVS(781)*UV(96)+JVS(841)*UV(98)+JVS(912)*UV(100)+JVS(1015)*UV(102) JTUV(62) = JVS(212)*UV(60)+JVS(237)*UV(62)+JVS(356)*UV(74)+JVS(720)*UV(95)+JVS(947)*UV(101)+JVS(1016)*UV(102) JTUV(63) = JVS(213)*UV(60)+JVS(248)*UV(63)+JVS(478)*UV(82)+JVS(697)*UV(94)+JVS(783)*UV(96)+JVS(1017)*UV(102) JTUV(64) = JVS(188)*UV(56)+JVS(249)*UV(63)+JVS(260)*UV(64)+JVS(289)*UV(67)+JVS(444)*UV(81)+JVS(479)*UV(82)+JVS(593)& &*UV(88)+JVS(721)*UV(95)+JVS(784)*UV(96)+JVS(867)*UV(99)+JVS(948)*UV(101)+JVS(1060)*UV(103) JTUV(65) = JVS(129)*UV(41)+JVS(267)*UV(65)+JVS(445)*UV(81)+JVS(480)*UV(82)+JVS(722)*UV(95)+JVS(785)*UV(96)+JVS(868)& &*UV(99)+JVS(949)*UV(101)+JVS(1061)*UV(103) JTUV(66) = JVS(133)*UV(42)+JVS(238)*UV(62)+JVS(277)*UV(66)+JVS(481)*UV(82)+JVS(723)*UV(95)+JVS(786)*UV(96)+JVS(869)& &*UV(99)+JVS(950)*UV(101)+JVS(1062)*UV(103) JTUV(67) = JVS(7)*UV(4)+JVS(63)*UV(22)+JVS(290)*UV(67)+JVS(870)*UV(99)+JVS(913)*UV(100)+JVS(1021)*UV(102) JTUV(68) = JVS(310)*UV(68)+JVS(446)*UV(81)+JVS(482)*UV(82)+JVS(531)*UV(84)+JVS(552)*UV(85)+JVS(561)*UV(86)+JVS(724)& &*UV(95)+JVS(787)*UV(96)+JVS(871)*UV(99)+JVS(951)*UV(101)+JVS(1063)*UV(103) JTUV(69) = JVS(189)*UV(56)+JVS(317)*UV(69)+JVS(447)*UV(81)+JVS(483)*UV(82)+JVS(594)*UV(88)+JVS(725)*UV(95)+JVS(788)& &*UV(96)+JVS(872)*UV(99)+JVS(952)*UV(101)+JVS(1064)*UV(103) JTUV(70) = JVS(194)*UV(57)+JVS(291)*UV(67)+JVS(324)*UV(70)+JVS(331)*UV(71)+JVS(484)*UV(82)+JVS(620)*UV(89)+JVS(726)& &*UV(95)+JVS(789)*UV(96)+JVS(873)*UV(99)+JVS(953)*UV(101)+JVS(1065)*UV(103) JTUV(71) = JVS(214)*UV(60)+JVS(250)*UV(63)+JVS(332)*UV(71)+JVS(388)*UV(77)+JVS(418)*UV(79)+JVS(485)*UV(82)+JVS(595)& &*UV(88)+JVS(790)*UV(96)+JVS(825)*UV(97)+JVS(842)*UV(98)+JVS(1024)*UV(102) JTUV(72) = JVS(176)*UV(53)+JVS(215)*UV(60)+JVS(292)*UV(67)+JVS(341)*UV(72)+JVS(389)*UV(77)+JVS(486)*UV(82)+JVS(596)& &*UV(88)+JVS(621)*UV(89)+JVS(727)*UV(95)+JVS(791)*UV(96)+JVS(874)*UV(99)+JVS(954)*UV(101)+JVS(1066)*UV(103) JTUV(73) = JVS(348)*UV(73)+JVS(487)*UV(82)+JVS(532)*UV(84)+JVS(562)*UV(86)+JVS(597)*UV(88)+JVS(622)*UV(89)+JVS(728)& &*UV(95)+JVS(792)*UV(96)+JVS(875)*UV(99)+JVS(955)*UV(101)+JVS(1067)*UV(103) JTUV(74) = JVS(358)*UV(74)+JVS(390)*UV(77)+JVS(488)*UV(82)+JVS(563)*UV(86)+JVS(598)*UV(88)+JVS(623)*UV(89)+JVS(729)& &*UV(95)+JVS(793)*UV(96)+JVS(876)*UV(99)+JVS(956)*UV(101)+JVS(1068)*UV(103) JTUV(75) = JVS(184)*UV(55)+JVS(251)*UV(63)+JVS(293)*UV(67)+JVS(367)*UV(75)+JVS(489)*UV(82)+JVS(599)*UV(88)+JVS(624)& &*UV(89)+JVS(730)*UV(95)+JVS(794)*UV(96)+JVS(877)*UV(99)+JVS(957)*UV(101)+JVS(1069)*UV(103) JTUV(76) = JVS(168)*UV(51)+JVS(376)*UV(76)+JVS(449)*UV(81)+JVS(490)*UV(82)+JVS(533)*UV(84)+JVS(731)*UV(95)+JVS(795)& &*UV(96)+JVS(878)*UV(99)+JVS(958)*UV(101)+JVS(1070)*UV(103) JTUV(77) = JVS(391)*UV(77)+JVS(491)*UV(82)+JVS(600)*UV(88)+JVS(796)*UV(96)+JVS(959)*UV(101)+JVS(1028)*UV(102) JTUV(78) = JVS(172)*UV(52)+JVS(294)*UV(67)+JVS(409)*UV(78)+JVS(450)*UV(81)+JVS(492)*UV(82)+JVS(661)*UV(91)+JVS(685)& &*UV(93)+JVS(732)*UV(95)+JVS(797)*UV(96)+JVS(879)*UV(99)+JVS(960)*UV(101)+JVS(1071)*UV(103) JTUV(79) = JVS(180)*UV(54)+JVS(216)*UV(60)+JVS(252)*UV(63)+JVS(295)*UV(67)+JVS(392)*UV(77)+JVS(419)*UV(79)+JVS(493)& &*UV(82)+JVS(601)*UV(88)+JVS(625)*UV(89)+JVS(733)*UV(95)+JVS(798)*UV(96)+JVS(880)*UV(99)+JVS(961)*UV(101)& &+JVS(1072)*UV(103) JTUV(80) = JVS(190)*UV(56)+JVS(253)*UV(63)+JVS(296)*UV(67)+JVS(393)*UV(77)+JVS(430)*UV(80)+JVS(451)*UV(81)+JVS(494)& &*UV(82)+JVS(734)*UV(95)+JVS(799)*UV(96)+JVS(881)*UV(99)+JVS(962)*UV(101)+JVS(1073)*UV(103) JTUV(81) = JVS(217)*UV(60)+JVS(297)*UV(67)+JVS(452)*UV(81)+JVS(575)*UV(87)+JVS(641)*UV(90)+JVS(800)*UV(96)+JVS(920)& &*UV(100)+JVS(1032)*UV(102) JTUV(82) = JVS(3)*UV(2)+JVS(64)*UV(22)+JVS(218)*UV(60)+JVS(298)*UV(67)+JVS(495)*UV(82)+JVS(801)*UV(96)+JVS(921)& &*UV(100)+JVS(1033)*UV(102) JTUV(83) = JVS(164)*UV(50)+JVS(394)*UV(77)+JVS(453)*UV(81)+JVS(496)*UV(82)+JVS(518)*UV(83)+JVS(534)*UV(84)+JVS(735)& &*UV(95)+JVS(802)*UV(96)+JVS(884)*UV(99)+JVS(963)*UV(101)+JVS(1074)*UV(103) JTUV(84) = JVS(219)*UV(60)+JVS(299)*UV(67)+JVS(497)*UV(82)+JVS(535)*UV(84)+JVS(736)*UV(95)+JVS(803)*UV(96)+JVS(923)& &*UV(100)+JVS(964)*UV(101)+JVS(1035)*UV(102) JTUV(85) = JVS(125)*UV(40)+JVS(239)*UV(62)+JVS(268)*UV(65)+JVS(278)*UV(66)+JVS(454)*UV(81)+JVS(498)*UV(82)+JVS(536)& &*UV(84)+JVS(553)*UV(85)+JVS(564)*UV(86)+JVS(577)*UV(87)+JVS(626)*UV(89)+JVS(737)*UV(95)+JVS(804)*UV(96)& &+JVS(886)*UV(99)+JVS(965)*UV(101)+JVS(1075)*UV(103) JTUV(86) = JVS(21)*UV(11)+JVS(65)*UV(22)+JVS(240)*UV(62)+JVS(269)*UV(65)+JVS(279)*UV(66)+JVS(311)*UV(68)+JVS(455)& &*UV(81)+JVS(537)*UV(84)+JVS(565)*UV(86)+JVS(578)*UV(87)+JVS(627)*UV(89)+JVS(738)*UV(95)+JVS(805)*UV(96)& &+JVS(887)*UV(99)+JVS(1037)*UV(102) JTUV(87) = JVS(148)*UV(46)+JVS(500)*UV(82)+JVS(538)*UV(84)+JVS(579)*UV(87)+JVS(739)*UV(95)+JVS(806)*UV(96)+JVS(888)& &*UV(99)+JVS(967)*UV(101)+JVS(1077)*UV(103) JTUV(88) = JVS(220)*UV(60)+JVS(300)*UV(67)+JVS(539)*UV(84)+JVS(606)*UV(88)+JVS(807)*UV(96)+JVS(927)*UV(100)+JVS(968)& &*UV(101)+JVS(1039)*UV(102) JTUV(89) = JVS(301)*UV(67)+JVS(349)*UV(73)+JVS(580)*UV(87)+JVS(628)*UV(89)+JVS(645)*UV(90)+JVS(741)*UV(95)+JVS(928)& &*UV(100)+JVS(969)*UV(101)+JVS(1040)*UV(102) JTUV(90) = JVS(77)*UV(27)+JVS(137)*UV(43)+JVS(502)*UV(82)+JVS(581)*UV(87)+JVS(646)*UV(90)+JVS(742)*UV(95)+JVS(809)& &*UV(96)+JVS(843)*UV(98)+JVS(891)*UV(99)+JVS(970)*UV(101)+JVS(1079)*UV(103) JTUV(91) = JVS(52)*UV(21)+JVS(221)*UV(60)+JVS(302)*UV(67)+JVS(318)*UV(69)+JVS(342)*UV(72)+JVS(503)*UV(82)+JVS(608)& &*UV(88)+JVS(662)*UV(91)+JVS(810)*UV(96)+JVS(826)*UV(97)+JVS(844)*UV(98)+JVS(930)*UV(100)+JVS(971)*UV(101)& &+JVS(1042)*UV(102) JTUV(92) = JVS(195)*UV(57)+JVS(325)*UV(70)+JVS(333)*UV(71)+JVS(504)*UV(82)+JVS(630)*UV(89)+JVS(663)*UV(91)+JVS(674)& &*UV(92)+JVS(686)*UV(93)+JVS(744)*UV(95)+JVS(811)*UV(96)+JVS(893)*UV(99)+JVS(972)*UV(101)+JVS(1081)*UV(103) JTUV(93) = JVS(222)*UV(60)+JVS(232)*UV(61)+JVS(368)*UV(75)+JVS(505)*UV(82)+JVS(542)*UV(84)+JVS(610)*UV(88)+JVS(687)& &*UV(93)+JVS(745)*UV(95)+JVS(812)*UV(96)+JVS(828)*UV(97)+JVS(846)*UV(98)+JVS(973)*UV(101)+JVS(1044)*UV(102) JTUV(94) = JVS(80)*UV(28)+JVS(152)*UV(47)+JVS(506)*UV(82)+JVS(706)*UV(94)+JVS(746)*UV(95)+JVS(813)*UV(96)+JVS(847)& &*UV(98)+JVS(895)*UV(99)+JVS(974)*UV(101)+JVS(1083)*UV(103) JTUV(95) = JVS(121)*UV(39)+JVS(223)*UV(60)+JVS(241)*UV(62)+JVS(254)*UV(63)+JVS(261)*UV(64)+JVS(270)*UV(65)+JVS(280)& &*UV(66)+JVS(303)*UV(67)+JVS(312)*UV(68)+JVS(319)*UV(69)+JVS(326)*UV(70)+JVS(334)*UV(71)+JVS(343)*UV(72)& &+JVS(350)*UV(73)+JVS(361)*UV(74)+JVS(369)*UV(75)+JVS(378)*UV(76)+JVS(399)*UV(77)+JVS(410)*UV(78)+JVS(421)& &*UV(79)+JVS(432)*UV(80)+JVS(458)*UV(81)+JVS(507)*UV(82)+JVS(520)*UV(83)+JVS(543)*UV(84)+JVS(555)*UV(85)& &+JVS(584)*UV(87)+JVS(611)*UV(88)+JVS(632)*UV(89)+JVS(650)*UV(90)+JVS(664)*UV(91)+JVS(675)*UV(92)+JVS(688)& &*UV(93)+JVS(707)*UV(94)+JVS(747)*UV(95)+JVS(814)*UV(96)+JVS(829)*UV(97)+JVS(848)*UV(98)+JVS(896)*UV(99)& &+JVS(975)*UV(101)+JVS(1084)*UV(103) JTUV(96) = JVS(33)*UV(15)+JVS(39)*UV(17)+JVS(45)*UV(19)+JVS(83)*UV(29)+JVS(91)*UV(31)+JVS(96)*UV(32)+JVS(101)*UV(33)& &+JVS(122)*UV(39)+JVS(126)*UV(40)+JVS(130)*UV(41)+JVS(134)*UV(42)+JVS(138)*UV(43)+JVS(144)*UV(45)+JVS(149)& &*UV(46)+JVS(153)*UV(47)+JVS(156)*UV(48)+JVS(160)*UV(49)+JVS(165)*UV(50)+JVS(169)*UV(51)+JVS(173)*UV(52)& &+JVS(177)*UV(53)+JVS(181)*UV(54)+JVS(185)*UV(55)+JVS(191)*UV(56)+JVS(196)*UV(57)+JVS(262)*UV(64)+JVS(271)& &*UV(65)+JVS(281)*UV(66)+JVS(313)*UV(68)+JVS(320)*UV(69)+JVS(327)*UV(70)+JVS(344)*UV(72)+JVS(351)*UV(73)& &+JVS(362)*UV(74)+JVS(370)*UV(75)+JVS(379)*UV(76)+JVS(411)*UV(78)+JVS(422)*UV(79)+JVS(433)*UV(80)+JVS(508)& &*UV(82)+JVS(521)*UV(83)+JVS(556)*UV(85)+JVS(568)*UV(86)+JVS(585)*UV(87)+JVS(612)*UV(88)+JVS(651)*UV(90)& &+JVS(676)*UV(92)+JVS(708)*UV(94)+JVS(748)*UV(95)+JVS(815)*UV(96)+JVS(830)*UV(97)+JVS(849)*UV(98)+JVS(897)& &*UV(99)+JVS(935)*UV(100)+JVS(976)*UV(101)+JVS(1047)*UV(102)+JVS(1085)*UV(103) JTUV(97) = JVS(145)*UV(45)+JVS(199)*UV(58)+JVS(509)*UV(82)+JVS(749)*UV(95)+JVS(816)*UV(96)+JVS(831)*UV(97)+JVS(850)& &*UV(98)+JVS(898)*UV(99)+JVS(977)*UV(101)+JVS(1086)*UV(103) JTUV(98) = JVS(13)*UV(7)+JVS(53)*UV(21)+JVS(66)*UV(22)+JVS(200)*UV(58)+JVS(204)*UV(59)+JVS(225)*UV(60)+JVS(233)*UV(61)& &+JVS(256)*UV(63)+JVS(336)*UV(71)+JVS(402)*UV(77)+JVS(510)*UV(82)+JVS(545)*UV(84)+JVS(614)*UV(88)+JVS(666)& &*UV(91)+JVS(690)*UV(93)+JVS(750)*UV(95)+JVS(817)*UV(96)+JVS(851)*UV(98)+JVS(899)*UV(99)+JVS(937)*UV(100)& &+JVS(1049)*UV(102)+JVS(1087)*UV(103) JTUV(99) = JVS(11)*UV(6)+JVS(67)*UV(22)+JVS(78)*UV(27)+JVS(81)*UV(28)+JVS(105)*UV(34)+JVS(114)*UV(37)+JVS(118)*UV(38)& &+JVS(161)*UV(49)+JVS(201)*UV(58)+JVS(305)*UV(67)+JVS(653)*UV(90)+JVS(710)*UV(94)+JVS(818)*UV(96)+JVS(833)& &*UV(97)+JVS(852)*UV(98)+JVS(900)*UV(99)+JVS(938)*UV(100)+JVS(979)*UV(101)+JVS(1050)*UV(102)+JVS(1088)& &*UV(103) JTUV(100) = JVS(87)*UV(30)+JVS(108)*UV(35)+JVS(111)*UV(36)+JVS(119)*UV(38)+JVS(141)*UV(44)+JVS(205)*UV(59)+JVS(226)& &*UV(60)+JVS(234)*UV(61)+JVS(306)*UV(67)+JVS(321)*UV(69)+JVS(352)*UV(73)+JVS(381)*UV(76)+JVS(413)*UV(78)& &+JVS(461)*UV(81)+JVS(512)*UV(82)+JVS(546)*UV(84)+JVS(557)*UV(85)+JVS(587)*UV(87)+JVS(616)*UV(88)+JVS(635)& &*UV(89)+JVS(654)*UV(90)+JVS(667)*UV(91)+JVS(752)*UV(95)+JVS(819)*UV(96)+JVS(834)*UV(97)+JVS(853)*UV(98)& &+JVS(901)*UV(99)+JVS(939)*UV(100)+JVS(980)*UV(101)+JVS(1051)*UV(102)+JVS(1089)*UV(103) JTUV(101) = JVS(54)*UV(21)+JVS(106)*UV(34)+JVS(157)*UV(48)+JVS(227)*UV(60)+JVS(243)*UV(62)+JVS(257)*UV(63)+JVS(263)& &*UV(64)+JVS(272)*UV(65)+JVS(282)*UV(66)+JVS(307)*UV(67)+JVS(314)*UV(68)+JVS(322)*UV(69)+JVS(328)*UV(70)& &+JVS(337)*UV(71)+JVS(345)*UV(72)+JVS(353)*UV(73)+JVS(363)*UV(74)+JVS(371)*UV(75)+JVS(382)*UV(76)+JVS(404)& &*UV(77)+JVS(414)*UV(78)+JVS(424)*UV(79)+JVS(435)*UV(80)+JVS(462)*UV(81)+JVS(513)*UV(82)+JVS(524)*UV(83)& &+JVS(547)*UV(84)+JVS(558)*UV(85)+JVS(588)*UV(87)+JVS(617)*UV(88)+JVS(636)*UV(89)+JVS(655)*UV(90)+JVS(668)& &*UV(91)+JVS(679)*UV(92)+JVS(692)*UV(93)+JVS(712)*UV(94)+JVS(753)*UV(95)+JVS(820)*UV(96)+JVS(835)*UV(97)& &+JVS(854)*UV(98)+JVS(902)*UV(99)+JVS(981)*UV(101)+JVS(1052)*UV(102)+JVS(1090)*UV(103) JTUV(102) = JVS(24)*UV(12)+JVS(27)*UV(13)+JVS(30)*UV(14)+JVS(55)*UV(21)+JVS(69)*UV(23)+JVS(71)*UV(24)+JVS(73)*UV(25)& &+JVS(75)*UV(26)+JVS(84)*UV(29)+JVS(88)*UV(30)+JVS(92)*UV(31)+JVS(97)*UV(32)+JVS(102)*UV(33)+JVS(109)*UV(35)& &+JVS(112)*UV(36)+JVS(115)*UV(37)+JVS(123)*UV(39)+JVS(127)*UV(40)+JVS(131)*UV(41)+JVS(135)*UV(42)+JVS(139)& &*UV(43)+JVS(142)*UV(44)+JVS(146)*UV(45)+JVS(150)*UV(46)+JVS(154)*UV(47)+JVS(158)*UV(48)+JVS(162)*UV(49)& &+JVS(166)*UV(50)+JVS(170)*UV(51)+JVS(174)*UV(52)+JVS(178)*UV(53)+JVS(182)*UV(54)+JVS(186)*UV(55)+JVS(192)& &*UV(56)+JVS(197)*UV(57)+JVS(202)*UV(58)+JVS(206)*UV(59)+JVS(228)*UV(60)+JVS(235)*UV(61)+JVS(244)*UV(62)& &+JVS(258)*UV(63)+JVS(273)*UV(65)+JVS(283)*UV(66)+JVS(308)*UV(67)+JVS(315)*UV(68)+JVS(338)*UV(71)+JVS(346)& &*UV(72)+JVS(354)*UV(73)+JVS(364)*UV(74)+JVS(372)*UV(75)+JVS(383)*UV(76)+JVS(405)*UV(77)+JVS(415)*UV(78)& &+JVS(425)*UV(79)+JVS(436)*UV(80)+JVS(463)*UV(81)+JVS(514)*UV(82)+JVS(525)*UV(83)+JVS(548)*UV(84)+JVS(559)& &*UV(85)+JVS(571)*UV(86)+JVS(589)*UV(87)+JVS(618)*UV(88)+JVS(637)*UV(89)+JVS(656)*UV(90)+JVS(669)*UV(91)& &+JVS(680)*UV(92)+JVS(693)*UV(93)+JVS(713)*UV(94)+JVS(754)*UV(95)+JVS(821)*UV(96)+JVS(836)*UV(97)+JVS(855)& &*UV(98)+JVS(903)*UV(99)+JVS(941)*UV(100)+JVS(982)*UV(101)+JVS(1053)*UV(102)+JVS(1091)*UV(103) JTUV(103) = JVS(36)*UV(16)+JVS(42)*UV(18)+JVS(48)*UV(20)+JVS(56)*UV(21)+JVS(93)*UV(31)+JVS(98)*UV(32)+JVS(103)*UV(33)& &+JVS(116)*UV(37)+JVS(229)*UV(60)+JVS(245)*UV(62)+JVS(259)*UV(63)+JVS(264)*UV(64)+JVS(274)*UV(65)+JVS(284)& &*UV(66)+JVS(309)*UV(67)+JVS(316)*UV(68)+JVS(323)*UV(69)+JVS(329)*UV(70)+JVS(339)*UV(71)+JVS(347)*UV(72)& &+JVS(355)*UV(73)+JVS(365)*UV(74)+JVS(373)*UV(75)+JVS(384)*UV(76)+JVS(406)*UV(77)+JVS(416)*UV(78)+JVS(426)& &*UV(79)+JVS(437)*UV(80)+JVS(464)*UV(81)+JVS(515)*UV(82)+JVS(526)*UV(83)+JVS(549)*UV(84)+JVS(560)*UV(85)& &+JVS(572)*UV(86)+JVS(590)*UV(87)+JVS(619)*UV(88)+JVS(638)*UV(89)+JVS(657)*UV(90)+JVS(670)*UV(91)+JVS(681)& &*UV(92)+JVS(694)*UV(93)+JVS(714)*UV(94)+JVS(755)*UV(95)+JVS(822)*UV(96)+JVS(837)*UV(97)+JVS(856)*UV(98)& &+JVS(904)*UV(99)+JVS(942)*UV(100)+JVS(983)*UV(101)+JVS(1054)*UV(102)+JVS(1092)*UV(103) END SUBROUTINE JacTR_SP_Vec ! End of JacTR_SP_Vec function ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ END MODULE gckpp_Jacobian