3583 lines
122 KiB
Fortran
3583 lines
122 KiB
Fortran
! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
|
|
!
|
|
! 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
|
|
|