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