cacm3_Function.f90.org

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! The ODE Function of Chemical Model File ! ! Generated by KPP-2.2.3 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 : cacm3_Function.f90 ! Time : Thu Jan 23 11:51:16 2020 ! Working directory : /home/dandanwe/boxmox/tmp_cacm3_4449_BOXMOX ! Equation file : cacm3.kpp ! Output root filename : cacm3 ! ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

MODULE cacm3_Function

USE cacm3_Parameters IMPLICIT NONE

! A - Rate for each equation REAL(kind=dp) :: A(NREACT)

CONTAINS

! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ! ! Fun - time derivatives of variables - Agregate form ! Arguments : ! V - Concentrations of variable species (local) ! F - Concentrations of fixed species (local) ! RCT - Rate constants (local) ! Vdot - Time derivative of variable species concentrations ! ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

SUBROUTINE Fun ( V, F, RCT, Vdot )

! 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) ! Vdot - Time derivative of variable species concentrations REAL(kind=dp) :: Vdot(NVAR)

! Computation of equation rates A(1) = RCT(1)*V(397)*F(1) A(2) = RCT(2)*V(391) A(3) = RCT(3)*V(391)*V(397) A(4) = RCT(4)*V(391)*V(397) A(5) = RCT(5)*V(401)*V(402) A(6) = RCT(6)*V(397)*V(401) A(7) = RCT(7)*V(400)*V(402) A(8) = RCT(8)*V(402)*V(402) A(9) = RCT(9)*V(397)*V(400) A(10) = RCT(10)*V(154) A(11) = RCT(11)*V(154)*F(3) A(12) = RCT(12)*V(397)*V(400) A(13) = RCT(13)*V(400)*F(1) A(14) = RCT(14)*V(400)*F(1) A(15) = RCT(15)*V(401)*F(1) A(16) = RCT(16)*V(401)*F(1) A(17) = RCT(17)*V(134)*F(3) A(18) = RCT(18)*V(134) A(19) = RCT(19)*V(402)*V(404) A(20) = RCT(20)*V(191)*F(1) A(21) = RCT(21)*V(397)*F(3) A(22) = RCT(22)*V(397)*V(404) A(23) = RCT(23)*V(270)*V(404) A(24) = RCT(24)*V(373)*V(404) A(25) = RCT(25)*V(401)*V(404) A(26) = RCT(26)*V(402)*V(407) A(27) = RCT(27)*V(397)*V(407) A(28) = RCT(28)*V(177) A(29) = RCT(29)*V(177)*V(404) A(30) = RCT(30)*V(401)*V(407) A(31) = RCT(31)*V(407)*V(407) A(32) = RCT(32)*V(407)*V(407)*F(3) A(33) = RCT(33)*V(400)*V(407) A(34) = RCT(34)*V(391)*V(401) A(35) = RCT(35)*V(133)*V(404) A(36) = RCT(36)*V(243)*F(1) A(37) = RCT(37)*V(243)*V(404) A(38) = RCT(38)*V(391)*V(402) A(39) = RCT(39)*V(191)*V(404) A(40) = RCT(40)*V(400)*V(404) A(41) = RCT(41)*V(400)*V(400) A(42) = RCT(42)*V(404)*V(407) A(43) = 5.42e-06*V(369)*F(1) A(44) = RCT(44)*V(369)*V(404) A(45) = RCT(45)*V(369)*V(404) A(46) = 5.42e-06*V(344)*F(1) A(47) = RCT(47)*V(344)*V(404) A(48) = RCT(48)*V(344)*V(404) A(49) = RCT(49)*V(393)*V(402) A(50) = RCT(50)*V(393)*V(393) A(51) = RCT(51)*V(393)*V(407) A(52) = RCT(52)*V(135)*V(404) A(53) = RCT(53)*V(320)*V(402) A(54) = RCT(54)*V(320)*V(393) A(55) = RCT(55)*V(320)*V(407) A(56) = RCT(56)*V(136)*V(404) A(57) = RCT(57)*V(237)*V(402) A(58) = RCT(58)*V(237)*V(393) A(59) = RCT(59)*V(237)*V(407) A(60) = RCT(60)*V(246)*V(404) A(61) = RCT(61)*V(371)*V(402) A(62) = RCT(62)*V(371)*V(402) A(63) = RCT(63)*V(371)*V(393) A(64) = RCT(64)*V(371)*V(407) A(65) = RCT(65)*V(203)*V(404) A(66) = RCT(66)*V(370)*V(402) A(67) = RCT(67)*V(370)*V(402) A(68) = RCT(68)*V(370)*V(393) A(69) = RCT(69)*V(370)*V(407) A(70) = RCT(70)*V(281)*V(402) A(71) = RCT(71)*V(281)*V(402) A(72) = RCT(72)*V(281)*V(393) A(73) = RCT(73)*V(281)*V(407) A(74) = RCT(74)*V(173)*V(404) A(75) = RCT(75)*V(155)*V(404) A(76) = RCT(76)*V(140)*V(404) A(77) = RCT(77)*V(273)*V(402) A(78) = RCT(78)*V(273)*V(402) A(79) = RCT(79)*V(273)*V(393) A(80) = RCT(80)*V(273)*V(407) A(81) = RCT(81)*V(242)*V(402) A(82) = RCT(82)*V(242)*V(402) A(83) = RCT(83)*V(242)*V(393) A(84) = RCT(84)*V(242)*V(407) A(85) = RCT(85)*V(169)*V(404) A(86) = RCT(86)*V(170)*V(404) A(87) = RCT(87)*V(345)*V(404) A(88) = RCT(88)*V(345)*V(400) A(89) = RCT(89)*V(345)*V(401) A(90) = RCT(90)*V(345)*V(391) A(91) = RCT(91)*V(283)*V(402) A(92) = RCT(92)*V(283)*V(393) A(93) = RCT(93)*V(283)*V(407) A(94) = RCT(94)*V(261)*V(402) A(95) = RCT(95)*V(261)*V(393) A(96) = RCT(96)*V(261)*V(407) A(97) = RCT(97)*V(337)*V(402) A(98) = RCT(98)*V(337)*V(393) A(99) = RCT(99)*V(337)*V(407) A(100) = RCT(100)*V(367)*V(404) A(101) = RCT(101)*V(367)*V(400) A(102) = RCT(102)*V(367)*V(401) A(103) = RCT(103)*V(367)*V(391) A(104) = RCT(104)*V(288)*V(402) A(105) = RCT(105)*V(288)*V(393) A(106) = RCT(106)*V(288)*V(407) A(107) = RCT(107)*V(293)*V(402) A(108) = RCT(108)*V(293)*V(393) A(109) = RCT(109)*V(293)*V(407) A(110) = RCT(110)*V(358)*V(404) A(111) = RCT(111)*V(358)*V(400) A(112) = RCT(112)*V(358)*V(401) A(113) = RCT(113)*V(358)*V(391) A(114) = RCT(114)*V(295)*V(402) A(115) = RCT(115)*V(295)*V(393) A(116) = RCT(116)*V(295)*V(407) A(117) = RCT(117)*V(375)*V(404) A(118) = RCT(118)*V(375)*V(400) A(119) = 7.35e-06*V(375)*F(1) A(120) = RCT(120)*V(336)*V(402) A(121) = RCT(121)*V(336)*V(397) A(122) = RCT(122)*V(144) A(123) = RCT(123)*V(336)*V(407) A(124) = RCT(124)*V(336)*V(393) A(125) = 3.02e-05*V(383)*F(1) A(126) = 4.63e-05*V(383)*F(1) A(127) = RCT(127)*V(383)*V(404) A(128) = RCT(128)*V(383)*V(400) A(129) = 4.87e-06*V(384)*F(1) A(130) = RCT(130)*V(384)*V(404) A(131) = RCT(131)*V(384)*V(400) A(132) = RCT(132)*V(392)*V(402) A(133) = RCT(133)*V(392)*V(397) A(134) = RCT(134)*V(146) A(135) = RCT(135)*V(392)*V(407) A(136) = RCT(136)*V(392)*V(393) A(137) = 7.35e-06*V(380)*F(1) A(138) = RCT(138)*V(380)*V(404) A(139) = RCT(139)*V(380)*V(400) A(140) = RCT(140)*V(339)*V(402) A(141) = RCT(141)*V(339)*V(397) A(142) = RCT(142)*V(147) A(143) = RCT(143)*V(339)*V(407) A(144) = RCT(144)*V(339)*V(393) A(145) = RCT(145)*V(192)*V(404) A(146) = RCT(146)*V(149)*V(404) A(147) = RCT(147)*V(149)*V(404) A(148) = RCT(148)*V(137)*V(404) A(149) = RCT(149)*V(363)*V(404) A(150) = 1.58e-06*V(363)*F(1) A(151) = RCT(151)*V(302)*V(402) A(152) = RCT(152)*V(302)*V(393) A(153) = RCT(153)*V(302)*V(407) A(154) = RCT(154)*V(390)*V(404) A(155) = 1.58e-06*V(390)*F(1) A(156) = RCT(156)*V(303)*V(402) A(157) = RCT(157)*V(303)*V(393) A(158) = RCT(158)*V(303)*V(407) A(159) = RCT(159)*V(156)*V(404) A(160) = RCT(160)*V(157) A(161) = RCT(161)*V(157)*V(397) A(162) = RCT(162)*V(247)*V(402) A(163) = RCT(163)*V(247)*V(402) A(164) = RCT(164)*V(247)*V(393) A(165) = RCT(165)*V(247)*V(407) A(166) = RCT(166)*V(321) A(167) = RCT(167)*V(321)*V(402) A(168) = RCT(168)*V(321)*V(393) A(169) = RCT(169)*V(321)*V(407) A(170) = RCT(170)*V(248)*V(402) A(171) = RCT(171)*V(248)*V(393) A(172) = RCT(172)*V(248)*V(407) A(173) = RCT(173)*V(211)*V(404) A(174) = RCT(174)*V(212)*V(404) A(175) = 7.35e-06*V(212)*F(1) A(176) = RCT(176)*V(396)*V(404) A(177) = RCT(177)*V(396)*V(400) A(178) = 0.000142*V(396)*F(1) A(179) = RCT(179)*V(335)*V(402) A(180) = RCT(180)*V(335)*V(397) A(181) = RCT(181)*V(158) A(182) = RCT(182)*V(335)*V(407) A(183) = RCT(183)*V(335)*V(393) A(184) = RCT(184)*V(178)*V(404) A(185) = RCT(185)*V(148)*V(404) A(186) = RCT(186)*V(159) A(187) = RCT(187)*V(159)*V(397) A(188) = RCT(188)*V(323) A(189) = RCT(189)*V(323)*V(402) A(190) = RCT(190)*V(323)*V(393) A(191) = RCT(191)*V(323)*V(407) A(192) = RCT(192)*V(249)*V(402) A(193) = RCT(193)*V(249)*V(393) A(194) = RCT(194)*V(249)*V(407) A(195) = RCT(195)*V(209)*V(400) A(196) = RCT(196)*V(209)*V(404) A(197) = RCT(197)*V(160)*V(397) A(198) = RCT(198)*V(195) A(199) = RCT(199)*V(195)*V(397) A(200) = RCT(200)*V(287)*V(402) A(201) = RCT(201)*V(287)*V(402) A(202) = RCT(202)*V(287)*V(393) A(203) = RCT(203)*V(287)*V(407) A(204) = RCT(204)*V(347) A(205) = RCT(205)*V(347)*V(402) A(206) = RCT(206)*V(347)*V(393) A(207) = RCT(207)*V(347)*V(407) A(208) = RCT(208)*V(250)*V(402) A(209) = RCT(209)*V(250)*V(393) A(210) = RCT(210)*V(250)*V(407) A(211) = RCT(211)*V(196)*V(400) A(212) = RCT(212)*V(138)*V(397) A(213) = RCT(213)*V(213)*V(404) A(214) = RCT(214)*V(204)*V(404) A(215) = RCT(215)*V(161)*V(404) A(216) = RCT(216)*V(179)*V(404) A(217) = RCT(217)*V(354)*V(400) A(218) = RCT(218)*V(354)*V(404) A(219) = RCT(219)*V(354)*V(404) A(220) = RCT(220)*V(354)*V(404) A(221) = RCT(221)*V(354)*V(404) A(222) = RCT(222)*V(214) A(223) = RCT(223)*V(214)*V(397) A(224) = RCT(224)*V(289)*V(402) A(225) = RCT(225)*V(289)*V(402) A(226) = RCT(226)*V(289)*V(393) A(227) = RCT(227)*V(289)*V(407) A(228) = RCT(228)*V(349) A(229) = RCT(229)*V(349)*V(402) A(230) = RCT(230)*V(349)*V(393) A(231) = RCT(231)*V(349)*V(407) A(232) = RCT(232)*V(251)*V(402) A(233) = RCT(233)*V(251)*V(393) A(234) = RCT(234)*V(251)*V(407) A(235) = RCT(235)*V(180)*V(404) A(236) = RCT(236)*V(216)*V(404) A(237) = RCT(237)*V(218)*V(404) A(238) = RCT(238)*V(215)*V(404) A(239) = RCT(239)*V(217)*V(404) A(240) = RCT(240)*V(162)*V(404) A(241) = RCT(241)*V(181)*V(404) A(242) = RCT(242)*V(219)*V(404) A(243) = RCT(243)*V(182) A(244) = RCT(244)*V(182)*V(397) A(245) = RCT(245)*V(316)*V(402) A(246) = RCT(246)*V(316)*V(402) A(247) = RCT(247)*V(316)*V(393) A(248) = RCT(248)*V(316)*V(407) A(249) = RCT(249)*V(355) A(250) = RCT(250)*V(355)*V(402) A(251) = RCT(251)*V(355)*V(393) A(252) = RCT(252)*V(355)*V(407) A(253) = RCT(253)*V(252)*V(402) A(254) = RCT(254)*V(252)*V(393) A(255) = RCT(255)*V(252)*V(407) A(256) = RCT(256)*V(183)*V(404) A(257) = RCT(257)*V(150)*V(404) A(258) = RCT(258)*V(163) A(259) = RCT(259)*V(163)*V(397) A(260) = RCT(260)*V(241)*V(402) A(261) = RCT(261)*V(241)*V(402) A(262) = RCT(262)*V(241)*V(393) A(263) = RCT(263)*V(241)*V(407) A(264) = RCT(264)*V(330) A(265) = RCT(265)*V(330)*V(402) A(266) = RCT(266)*V(330)*V(393) A(267) = RCT(267)*V(330)*V(407) A(268) = RCT(268)*V(272)*V(402) A(269) = RCT(269)*V(272)*V(393) A(270) = RCT(270)*V(272)*V(407) A(271) = RCT(271)*V(205)*V(404) A(272) = RCT(272)*V(206)*V(404) A(273) = RCT(273)*V(164)*V(404) A(274) = RCT(274)*V(165)*V(404) A(275) = RCT(275)*V(139)*V(404) A(276) = RCT(276)*V(331)*V(402) A(277) = RCT(277)*V(331)*V(393) A(278) = RCT(278)*V(331)*V(407) A(279) = RCT(279)*V(372)*V(404) A(280) = RCT(280)*V(372)*V(400) A(281) = RCT(281)*V(372)*V(401) A(282) = RCT(282)*V(372)*V(391) A(283) = RCT(283)*V(313)*V(402) A(284) = RCT(284)*V(313)*V(393) A(285) = RCT(285)*V(313)*V(407) A(286) = RCT(286)*V(334)*V(402) A(287) = RCT(287)*V(334)*V(393) A(288) = RCT(288)*V(334)*V(407) A(289) = RCT(289)*V(310)*V(402) A(290) = RCT(290)*V(310)*V(402) A(291) = RCT(291)*V(310)*V(393) A(292) = RCT(292)*V(310)*V(407) A(293) = RCT(293)*V(389)*V(402) A(294) = RCT(294)*V(389)*V(393) A(295) = RCT(295)*V(389)*V(407) A(296) = RCT(296)*V(343)*V(402) A(297) = RCT(297)*V(343)*V(393) A(298) = RCT(298)*V(343)*V(407) A(299) = 7.35e-06*V(379)*F(1) A(300) = RCT(300)*V(379)*V(404) A(301) = RCT(301)*V(379)*V(400) A(302) = RCT(302)*V(353)*V(402) A(303) = RCT(303)*V(353)*V(397) A(304) = RCT(304)*V(166) A(305) = RCT(305)*V(353)*V(393) A(306) = RCT(306)*V(353)*V(407) A(307) = 7.35e-06*V(388)*F(1) A(308) = RCT(308)*V(388)*V(404) A(309) = RCT(309)*V(388)*V(400) A(310) = RCT(310)*V(351)*V(402) A(311) = RCT(311)*V(351)*V(397) A(312) = RCT(312)*V(167) A(313) = RCT(313)*V(351)*V(393) A(314) = RCT(314)*V(351)*V(407) A(315) = RCT(315)*V(301)*V(402) A(316) = RCT(316)*V(301)*V(402) A(317) = RCT(317)*V(301)*V(393) A(318) = RCT(318)*V(301)*V(407) A(319) = RCT(319)*V(376)*V(402) A(320) = RCT(320)*V(376)*V(397) A(321) = RCT(321)*V(168) A(322) = RCT(322)*V(376)*V(393) A(323) = RCT(323)*V(376)*V(407) A(324) = RCT(324)*V(184)*V(404) A(325) = RCT(325)*V(185)*V(404) A(326) = RCT(326)*V(186)*V(404) A(327) = RCT(327)*V(184) A(328) = RCT(328)*V(31) A(329) = RCT(329)*V(32) A(330) = RCT(330)*V(33) A(331) = RCT(331)*V(34) A(332) = RCT(332)*V(35) A(333) = RCT(333)*V(359)*V(404) A(334) = RCT(334)*V(359)*V(400) A(335) = RCT(335)*V(359)*V(401) A(336) = RCT(336)*V(359)*V(391) A(337) = RCT(337)*V(278)*V(402) A(338) = RCT(338)*V(278)*V(393) A(339) = RCT(339)*V(278)*V(407) A(340) = RCT(340)*V(342)*V(402) A(341) = RCT(341)*V(342)*V(393) A(342) = RCT(342)*V(342)*V(407) A(343) = RCT(343)*V(279)*V(404) A(344) = RCT(344)*V(365)*V(402) A(345) = RCT(345)*V(365)*V(393) A(346) = RCT(346)*V(365)*V(407) A(347) = RCT(347)*V(368)*V(404) A(348) = RCT(348)*V(368)*V(400) A(349) = RCT(349)*V(368)*V(401) A(350) = RCT(350)*V(368)*V(391) A(351) = RCT(351)*V(265)*V(402) A(352) = RCT(352)*V(265)*V(393) A(353) = RCT(353)*V(265)*V(407) A(354) = RCT(354)*V(296)*V(402) A(355) = RCT(355)*V(296)*V(393) A(356) = RCT(356)*V(296)*V(407) A(357) = RCT(357)*V(332)*V(402) A(358) = RCT(358)*V(332)*V(393) A(359) = RCT(359)*V(332)*V(407) A(360) = RCT(360)*V(309)*V(402) A(361) = RCT(361)*V(309)*V(393) A(362) = RCT(362)*V(309)*V(407) A(363) = RCT(363)*V(341)*V(402) A(364) = RCT(364)*V(341)*V(393) A(365) = RCT(365)*V(341)*V(407) A(366) = RCT(366)*V(312)*V(402) A(367) = RCT(367)*V(312)*V(393) A(368) = RCT(368)*V(312)*V(407) A(369) = RCT(369)*V(382)*V(402) A(370) = RCT(370)*V(382)*V(393) A(371) = RCT(371)*V(382)*V(407) A(372) = RCT(372)*V(311)*V(402) A(373) = RCT(373)*V(311)*V(393) A(374) = RCT(374)*V(311)*V(407) A(375) = RCT(375)*V(260)*V(402) A(376) = RCT(376)*V(260)*V(397) A(377) = RCT(377)*V(141) A(378) = RCT(378)*V(260)*V(393) A(379) = RCT(379)*V(260)*V(407) A(380) = RCT(380)*V(385)*V(404) A(381) = RCT(381)*V(385)*V(400) A(382) = RCT(382)*V(385)*V(401) A(383) = 7.35e-06*V(385)*F(1) A(384) = RCT(384)*V(282)*V(402) A(385) = RCT(385)*V(282)*V(393) A(386) = RCT(386)*V(282)*V(407) A(387) = RCT(387)*V(306)*V(402) A(388) = RCT(388)*V(306)*V(393) A(389) = RCT(389)*V(306)*V(407) A(390) = RCT(390)*V(322)*V(402) A(391) = RCT(391)*V(322)*V(393) A(392) = RCT(392)*V(322)*V(407) A(393) = 1.58e-06*V(378)*F(1) A(394) = RCT(394)*V(378)*V(404) A(395) = RCT(395)*V(378)*V(400) A(396) = RCT(396)*V(378)*V(401) A(397) = RCT(397)*V(378)*V(391) A(398) = RCT(398)*V(258)*V(402) A(399) = RCT(399)*V(258)*V(393) A(400) = RCT(400)*V(258)*V(407) A(401) = RCT(401)*V(284)*V(402) A(402) = RCT(402)*V(284)*V(393) A(403) = RCT(403)*V(284)*V(407) A(404) = RCT(404)*V(315)*V(402) A(405) = RCT(405)*V(315)*V(393) A(406) = RCT(406)*V(315)*V(407) A(407) = RCT(407)*V(350)*V(402) A(408) = RCT(408)*V(350)*V(393) A(409) = RCT(409)*V(350)*V(407) A(410) = RCT(410)*V(305)*V(402) A(411) = RCT(411)*V(305)*V(393) A(412) = RCT(412)*V(305)*V(407) A(413) = RCT(413)*V(264)*V(402) A(414) = RCT(414)*V(264)*V(397) A(415) = RCT(415)*V(142) A(416) = RCT(416)*V(264)*V(393) A(417) = RCT(417)*V(264)*V(407) A(418) = RCT(418)*V(328)*V(404) A(419) = RCT(419)*V(254)*V(404) A(420) = RCT(420)*V(256)*V(404) A(421) = RCT(421)*V(174) A(422) = RCT(422)*V(194) A(423) = RCT(423)*V(193) A(424) = RCT(424)*V(175) A(425) = RCT(425)*V(36) A(426) = RCT(426)*V(37) A(427) = RCT(427)*V(38) A(428) = RCT(428)*V(39) A(429) = RCT(429)*V(40) A(430) = RCT(430)*V(41) A(431) = RCT(431)*V(42) A(432) = RCT(432)*V(43) A(433) = RCT(433)*V(44) A(434) = RCT(434)*V(357)*V(404) A(435) = RCT(435)*V(357)*V(400) A(436) = RCT(436)*V(357)*V(401) A(437) = RCT(437)*V(357)*V(391) A(438) = RCT(438)*V(266)*V(402) A(439) = RCT(439)*V(266)*V(402) A(440) = RCT(440)*V(266)*V(393) A(441) = RCT(441)*V(266)*V(407) A(442) = RCT(442)*V(267)*V(402) A(443) = RCT(443)*V(267)*V(393) A(444) = RCT(444)*V(267)*V(407) A(445) = RCT(445)*V(297)*V(402) A(446) = RCT(446)*V(297)*V(393) A(447) = RCT(447)*V(297)*V(407) A(448) = RCT(448)*V(366)*V(404) A(449) = RCT(449)*V(366)*V(400) A(450) = 7.35e-06*V(366)*F(1) A(451) = RCT(451)*V(269)*V(402) A(452) = RCT(452)*V(269)*V(397) A(453) = RCT(453)*V(143) A(454) = RCT(454)*V(269)*V(407) A(455) = RCT(455)*V(269)*V(393) A(456) = RCT(456)*V(45) A(457) = RCT(457)*V(360)*V(404) A(458) = RCT(458)*V(360)*V(400) A(459) = RCT(459)*V(360)*V(401) A(460) = RCT(460)*V(360)*V(391) A(461) = RCT(461)*V(238)*V(402) A(462) = RCT(462)*V(238)*V(402) A(463) = RCT(463)*V(238)*V(393) A(464) = RCT(464)*V(238)*V(407) A(465) = RCT(465)*V(239)*V(402) A(466) = RCT(466)*V(239)*V(393) A(467) = RCT(467)*V(239)*V(407) A(468) = RCT(468)*V(304)*V(402) A(469) = RCT(469)*V(304)*V(402) A(470) = RCT(470)*V(304)*V(393) A(471) = RCT(471)*V(304)*V(407) A(472) = RCT(472)*V(286)*V(402) A(473) = RCT(473)*V(286)*V(393) A(474) = RCT(474)*V(286)*V(407) A(475) = RCT(475)*V(271)*V(402) A(476) = RCT(476)*V(271)*V(393) A(477) = RCT(477)*V(271)*V(407) A(478) = RCT(478)*V(374)*V(404) A(479) = RCT(479)*V(374)*V(400) A(480) = 7.35e-06*V(374)*F(1) A(481) = 7.35e-06*V(374)*F(1) A(482) = RCT(482)*V(338)*V(402) A(483) = RCT(483)*V(338)*V(397) A(484) = RCT(484)*V(151) A(485) = RCT(485)*V(338)*V(407) A(486) = RCT(486)*V(338)*V(393) A(487) = RCT(487)*V(132) A(488) = RCT(488)*V(46) A(489) = RCT(489)*V(47) A(490) = RCT(490)*V(340)*V(404) A(491) = RCT(491)*V(340)*V(401) A(492) = RCT(492)*V(340)*V(400) A(493) = RCT(493)*V(348)*V(402) A(494) = RCT(494)*V(348)*V(393) A(495) = RCT(495)*V(348)*V(407) A(496) = RCT(496)*V(307)*V(402) A(497) = RCT(497)*V(307)*V(402) A(498) = RCT(498)*V(307)*V(393) A(499) = RCT(499)*V(307)*V(407) A(500) = RCT(500)*V(308)*V(402) A(501) = RCT(501)*V(308)*V(393) A(502) = RCT(502)*V(308)*V(407) A(503) = 6.4e-06*V(236)*F(1) A(504) = 6.4e-06*V(197)*F(1) A(505) = RCT(505)*V(236)*V(404) A(506) = RCT(506)*V(197)*V(404) A(507) = RCT(507)*V(236)*V(404) A(508) = RCT(508)*V(197)*V(404) A(509) = RCT(509)*V(176)*V(404) A(510) = 2e-10*V(300)*V(404) A(511) = 1.2e-14*V(300)*V(401) A(512) = 2.2e-11*V(300)*V(400) A(513) = RCT(513)*V(262)*V(402) A(514) = RCT(514)*V(262)*V(402) A(515) = RCT(515)*V(262)*V(407) A(516) = 9e-14*V(262)*V(393) A(517) = RCT(517)*V(290)*V(402) A(518) = RCT(518)*V(290)*V(407) A(519) = 9e-14*V(290)*V(393) A(520) = RCT(520)*V(240)*V(402) A(521) = RCT(521)*V(240)*V(407) A(522) = 9e-14*V(240)*V(393) A(523) = 3.2e-10*V(198)*V(404) A(524) = 1e-15*V(198)*V(401) A(525) = 2.2e-11*V(198)*V(400) A(526) = 2e-10*V(199)*V(404) A(527) = 1.2e-14*V(199)*V(401) A(528) = 2.2e-11*V(199)*V(400) A(529) = RCT(529)*V(291)*V(402) A(530) = RCT(530)*V(291)*V(407) A(531) = 9e-14*V(291)*V(393) A(532) = RCT(532)*V(292)*V(402) A(533) = RCT(533)*V(292)*V(407) A(534) = 9e-14*V(292)*V(393) A(535) = RCT(535)*V(188)*V(404) A(536) = RCT(536)*V(356)*V(404) A(537) = RCT(537)*V(356)*V(404) A(538) = RCT(538)*V(356)*V(404) A(539) = RCT(539)*V(356)*V(404) A(540) = RCT(540)*V(399)*V(399) A(541) = RCT(541)*V(398)*V(398) A(542) = RCT(542)*V(398)*V(399) A(543) = RCT(543)*V(399)*V(399) A(544) = RCT(544)*V(398)*V(398) A(545) = RCT(545)*V(398)*V(399) A(546) = RCT(546)*V(399)*V(402) A(547) = RCT(547)*V(399)*V(402) A(548) = RCT(548)*V(399)*V(402) A(549) = RCT(549)*V(399)*V(402) A(550) = RCT(550)*V(398)*V(402) A(551) = RCT(551)*V(398)*V(402) A(552) = RCT(552)*V(398)*V(402) A(553) = RCT(553)*V(398)*V(402) A(554) = RCT(554)*V(399)*V(407) A(555) = RCT(555)*V(399)*V(407) A(556) = RCT(556)*V(398)*V(407) A(557) = RCT(557)*V(398)*V(407) A(558) = RCT(558)*V(377)*V(399) A(559) = RCT(559)*V(377)*V(399) A(560) = RCT(560)*V(377)*V(398) A(561) = RCT(561)*V(377)*V(398) A(562) = RCT(562)*V(399) A(563) = RCT(563)*V(399) A(564) = RCT(564)*V(398) A(565) = RCT(565)*V(398) A(566) = RCT(566)*V(220)*V(404) A(567) = RCT(567)*V(220) A(568) = RCT(568)*V(221)*V(404) A(569) = RCT(569)*V(221) A(570) = RCT(570)*V(294)*V(404) A(571) = RCT(571)*V(298)*V(404) A(572) = RCT(572)*V(314)*V(404) A(573) = RCT(573)*V(314)*V(404) A(574) = RCT(574)*V(314)*V(404) A(575) = RCT(575)*V(314) A(576) = RCT(576)*V(299)*V(404) A(577) = RCT(577)*V(299)*V(404) A(578) = RCT(578)*V(299)*V(404) A(579) = RCT(579)*V(299) A(580) = RCT(580)*V(145)*V(404) A(581) = RCT(581)*V(352)*V(407) A(582) = RCT(582)*V(352)*V(402) A(583) = RCT(583)*V(352)*V(402) A(584) = RCT(584)*V(352) A(585) = RCT(585)*V(326)*V(404) A(586) = RCT(586)*V(319)*V(404) A(587) = RCT(587)*V(268) A(588) = RCT(588)*V(268)*V(402) A(589) = RCT(589)*V(268)*V(402) A(590) = RCT(590)*V(268)*V(407) A(591) = RCT(591)*V(259) A(592) = RCT(592)*V(259)*V(402) A(593) = RCT(593)*V(259)*V(402) A(594) = RCT(594)*V(259)*V(407) A(595) = RCT(595)*V(222) A(596) = RCT(596)*V(222)*V(402) A(597) = RCT(597)*V(222)*V(402) A(598) = RCT(598)*V(222)*V(407) A(599) = RCT(599)*V(223)*V(404) A(600) = RCT(600)*V(276)*V(404) A(601) = RCT(601)*V(276)*V(404) A(602) = RCT(602)*V(280)*V(404) A(603) = RCT(603)*V(280)*V(404) A(604) = RCT(604)*V(327) A(605) = RCT(605)*V(327) A(606) = RCT(606)*V(327)*V(402) A(607) = RCT(607)*V(327)*V(402) A(608) = RCT(608)*V(327)*V(407) A(609) = RCT(609)*V(346) A(610) = RCT(610)*V(346) A(611) = RCT(611)*V(346)*V(407) A(612) = RCT(612)*V(346)*V(402) A(613) = RCT(613)*V(346)*V(402) A(614) = RCT(614)*V(245)*V(404) A(615) = RCT(615)*V(190)*V(404) A(616) = RCT(616)*V(190) A(617) = RCT(617)*V(277)*V(404) A(618) = RCT(618)*V(277) A(619) = RCT(619)*V(207)*V(404) A(620) = RCT(620)*V(153)*V(404) A(621) = RCT(621)*V(275)*V(404) A(622) = RCT(622)*V(395)*V(407) A(623) = RCT(623)*V(395)*V(402) A(624) = RCT(624)*V(395)*V(402) A(625) = RCT(625)*V(364)*V(404) A(626) = RCT(626)*V(364)*V(404) A(627) = RCT(627)*V(387)*V(407) A(628) = RCT(628)*V(394)*V(407) A(629) = RCT(629)*V(387) A(630) = RCT(630)*V(387)*V(402) A(631) = RCT(631)*V(387)*V(402) A(632) = RCT(632)*V(394)*V(402) A(633) = RCT(633)*V(394)*V(397) A(634) = RCT(634)*V(202) A(635) = 2.9e-11*V(202)*V(404) A(636) = RCT(636)*V(324) A(637) = RCT(637)*V(329) A(638) = RCT(638)*V(324)*V(404) A(639) = RCT(639)*V(329)*V(404) A(640) = RCT(640)*V(263)*V(402) A(641) = RCT(641)*V(263)*V(407) A(642) = 9e-12*V(263)*V(397) A(643) = RCT(643)*V(253)*V(402) A(644) = RCT(644)*V(253)*V(407) A(645) = 9e-12*V(253)*V(397) A(646) = RCT(646)*V(208) A(647) = RCT(647)*V(235) A(648) = RCT(648)*V(208)*V(404) A(649) = RCT(649)*V(235)*V(404) A(650) = RCT(650)*V(257)*V(404) A(651) = RCT(651)*V(257) A(652) = RCT(652)*V(224)*V(402) A(653) = RCT(653)*V(224)*V(407) A(654) = RCT(654)*V(225)*V(402) A(655) = RCT(655)*V(225)*V(407) A(656) = RCT(656)*V(200)*V(404) A(657) = RCT(657)*V(200)*V(404) A(658) = RCT(658)*V(201)*V(404) A(659) = RCT(659)*V(201)*V(404) A(660) = RCT(660)*V(226) A(661) = RCT(661)*V(226)*V(407) A(662) = RCT(662)*V(226)*V(402) A(663) = RCT(663)*V(226)*V(402) A(664) = RCT(664)*V(227) A(665) = RCT(665)*V(227)*V(407) A(666) = RCT(666)*V(227)*V(402) A(667) = RCT(667)*V(227)*V(402) A(668) = 1.3e-17*V(356)*V(401) A(669) = 1.7e-15*V(187)*F(3) A(670) = RCT(670)*V(187)*F(3)*F(3) A(671) = RCT(671)*V(187)*V(401) A(672) = RCT(672)*V(152)*V(404) A(673) = RCT(673)*V(356)*V(400) A(674) = RCT(674)*V(406)*V(407) A(675) = RCT(675)*V(405)*V(407) A(676) = 2.56e-12*V(405)*V(406) A(677) = 3.71e-12*V(405)*V(405) A(678) = 1.61e-12*V(406)*V(406) A(679) = 2.8e-13*V(377)*V(406) A(680) = 1.18e-12*V(377)*V(405) A(681) = 1.92e-12*V(403)*V(406) A(682) = 7.71e-12*V(403)*V(405) A(683) = 2.3e-12*V(400)*V(406) A(684) = 2.3e-12*V(400)*V(405) A(685) = RCT(685)*V(402)*V(406) A(686) = RCT(686)*V(402)*V(406) A(687) = RCT(687)*V(402)*V(405) A(688) = RCT(688)*V(402)*V(405) A(689) = RCT(689)*V(210)*F(2) A(690) = RCT(690)*V(210) A(691) = RCT(691)*V(171)*V(404) A(692) = RCT(692)*V(362)*V(404) A(693) = RCT(693)*V(362)*V(404) A(694) = RCT(694)*V(362)*V(404) A(695) = RCT(695)*V(361)*V(404) A(696) = RCT(696)*V(361)*V(404) A(697) = RCT(697)*V(361)*V(404) A(698) = RCT(698)*V(229)*V(407) A(699) = RCT(699)*V(229)*V(402) A(700) = RCT(700)*V(229)*V(402) A(701) = RCT(701)*V(229) A(702) = RCT(702)*V(230)*V(407) A(703) = RCT(703)*V(230)*V(402) A(704) = RCT(704)*V(230)*V(402) A(705) = RCT(705)*V(230) A(706) = RCT(706)*V(333)*V(407) A(707) = RCT(707)*V(333)*V(402) A(708) = RCT(708)*V(333)*V(402) A(709) = RCT(709)*V(189)*V(404) A(710) = RCT(710)*V(285)*V(404) A(711) = RCT(711)*V(189)*V(404) A(712) = RCT(712)*V(285)*V(404) A(713) = RCT(713)*V(285)*V(404) A(714) = RCT(714)*V(189)*V(404) A(715) = RCT(715)*V(285)*V(404) A(716) = RCT(716)*V(285)*V(404) A(717) = RCT(717)*V(232) A(718) = RCT(718)*V(231) A(719) = RCT(719)*V(231)*V(407) A(720) = RCT(720)*V(232)*V(407) A(721) = RCT(721)*V(231)*V(402) A(722) = RCT(722)*V(231)*V(402) A(723) = RCT(723)*V(232)*V(402) A(724) = RCT(724)*V(232)*V(402) A(725) = RCT(725)*V(317)*V(404) A(726) = RCT(726)*V(317)*V(404) A(727) = RCT(727)*V(244)*V(407) A(728) = RCT(728)*V(244)*V(402) A(729) = RCT(729)*V(244)*V(402) A(730) = RCT(730)*V(386)*V(404) A(731) = RCT(731)*V(386)*V(404) A(732) = RCT(732)*V(381)*V(404) A(733) = RCT(733)*V(381)*V(404) A(734) = RCT(734)*V(228)*V(404) A(735) = RCT(735)*V(255) A(736) = RCT(736)*V(255)*V(407) A(737) = RCT(737)*V(255)*V(402) A(738) = RCT(738)*V(255)*V(402) A(739) = RCT(739)*V(274) A(740) = RCT(740)*V(274)*V(407) A(741) = RCT(741)*V(274)*V(402) A(742) = RCT(742)*V(274)*V(402) A(743) = RCT(743)*V(318)*V(407) A(744) = RCT(744)*V(318)*V(402) A(745) = RCT(745)*V(318)*V(402) A(746) = RCT(746)*V(386) A(747) = RCT(747)*V(381) A(748) = RCT(748)*V(228) A(749) = RCT(749)*V(386)*V(404) A(750) = RCT(750)*V(381)*V(404) A(751) = RCT(751)*V(233)*V(402) A(752) = RCT(752)*V(233)*V(407) A(753) = RCT(753)*V(234)*V(402) A(754) = RCT(754)*V(234)*V(407) A(755) = RCT(755)*V(325)*V(404) A(756) = RCT(756)*V(172)*V(404)

! ddw - altered production and loss output files; added to pass chemical reaction rates to output ! do ii=1,NREACT ! chemR8(ii)=A(ii) ! enddo

! Aggregate function Vdot(1) = A(24)+0.2*A(102)+0.2*A(112)+A(120)+A(124)+A(132)+A(136)+A(140)+A(144)+A(179)+A(183)+A(302)+0.8*A(305)+A(310)& &+0.8*A(313)+A(319)+0.6*A(322)+0.16*A(335)+A(375)+0.5*A(378)+A(413)+0.5*A(416)+0.9*A(418)+0.9*A(419)+0.9& &*A(420)+A(451)+A(455)+A(482)+A(486)+0.091*A(491)+0.355*A(570)+0.355*A(571)+A(615)+0.63*A(619)+A(620)+0.5& &*A(628)+A(632)+0.13*A(648)+A(652)+A(653)+A(654)+A(655)+0.407*A(668)+A(681)+A(682)+0.33*A(751)+0.33*A(752)& &+0.33*A(753)+0.33*A(754) Vdot(2) = A(35) Vdot(3) = A(71)+A(72)+A(73)+A(74)+A(75)+A(123) Vdot(4) = A(82)+A(83)+A(84)+A(85)+A(86) Vdot(5) = 0.12*A(112) Vdot(6) = 0.37*A(89)+0.24*A(102)+0.12*A(112)+A(135)+A(143)+0.11*A(335)+0.036*A(349)+0.09*A(382)+0.285*A(491)+0.2& &*A(510) Vdot(7) = A(161) Vdot(8) = 0.67*A(173) Vdot(9) = 0.7*A(167)+0.7*A(168)+0.7*A(169)+0.33*A(173)+A(174)+0.5*A(189)+0.5*A(190)+0.5*A(191)+0.7*A(205)+0.7*A(206)& &+0.7*A(207)+0.7*A(229)+0.7*A(230)+0.7*A(231)+0.7*A(250)+0.7*A(251)+0.7*A(252) Vdot(10) = A(175) Vdot(11) = A(187) Vdot(12) = A(212) Vdot(13) = A(213)+0.16*A(242) Vdot(14) = A(215) Vdot(15) = A(235)+A(244) Vdot(16) = A(240) Vdot(17) = A(238) Vdot(18) = 0.16*A(257) Vdot(19) = A(259) Vdot(20) = A(261)+A(262)+A(263)+A(274) Vdot(21) = A(272) Vdot(22) = A(273) Vdot(23) = 0.2*A(283)+0.135*A(287) Vdot(24) = 0.175*A(286)+0.07*A(287)+0.175*A(288) Vdot(25) = A(289) Vdot(26) = A(315) Vdot(27) = 0.75*A(282) Vdot(28) = 0.3*A(284) Vdot(29) = 0.05*A(297) Vdot(30) = A(327)+A(328)+A(329)+A(330)+A(331)+A(332)+A(425)+A(426)+A(427)+A(428)+A(429)+A(430)+A(431)+A(432)+A(433)& &+A(456)+A(488)+A(489) Vdot(31) = 0.05*A(281)+0.2*A(305)+A(306)-A(328) Vdot(32) = 0.2*A(313)+A(314)-A(329) Vdot(33) = 0.1*A(322)+A(323)+A(324)-A(330) Vdot(34) = 0.1*A(297)+A(326)-A(331) Vdot(35) = A(325)-A(332) Vdot(36) = 0.033*A(349)-A(425) Vdot(37) = 0.055*A(396)+0.1*A(418)-A(426) Vdot(38) = 0.1*A(367)-A(427) Vdot(39) = 0.3*A(378)-A(428) Vdot(40) = 0.2*A(378)+A(379)-A(429) Vdot(41) = 0.1*A(411)-A(430) Vdot(42) = 0.3*A(416)-A(431) Vdot(43) = 0.2*A(416)+A(417)-A(432) Vdot(44) = 0.1*A(419)+0.1*A(420)-A(433) Vdot(45) = 0.055*A(436)+A(454)-A(456) Vdot(46) = 0.055*A(459)-A(488) Vdot(47) = A(485)-A(489) Vdot(48) = 0.35*A(337)+0.28*A(341) Vdot(49) = 0.1*A(340)+0.142*A(341)+0.17*A(342) Vdot(50) = 0.75*A(336) Vdot(51) = 0.3*A(338) Vdot(52) = 0.4*A(340)+0.26*A(341) Vdot(53) = 0.35*A(345) Vdot(54) = 0.3*A(345) Vdot(55) = 0.12*A(369) Vdot(56) = 0.36*A(384)+0.28*A(388) Vdot(57) = 0.17*A(387)+0.142*A(388)+0.17*A(389) Vdot(58) = 0.35*A(398)+0.265*A(402) Vdot(59) = 0.175*A(401)+0.175*A(402)+0.175*A(403) Vdot(60) = 0.1*A(404) Vdot(61) = 0.75*A(350) Vdot(62) = 0.25*A(350) Vdot(63) = 0.75*A(397) Vdot(64) = 0.25*A(397) Vdot(65) = 0.3*A(352) Vdot(66) = 0.3*A(355) Vdot(67) = 0.2*A(364)+0.1*A(367) Vdot(68) = 0.2*A(364)+0.1*A(367) Vdot(69) = 0.3*A(385) Vdot(70) = 0.3*A(391) Vdot(71) = 0.2*A(408)+0.1*A(411) Vdot(72) = 0.3*A(399) Vdot(73) = 0.2*A(405) Vdot(74) = 0.2*A(405) Vdot(75) = 0.2*A(408)+0.1*A(411) Vdot(76) = A(438) Vdot(77) = 0.445*A(436)+A(450)+A(451)+A(455) Vdot(78) = 0.75*A(437) Vdot(79) = 0.25*A(437) Vdot(80) = A(445)+A(446)+A(447) Vdot(81) = 0.055*A(459)+0.3*A(462)+0.3*A(463)+0.3*A(464)+0.3*A(465)+0.3*A(466)+0.3*A(467) Vdot(82) = 0.75*A(460)+A(487) Vdot(83) = 0.25*A(460)+0.7*A(462)+0.7*A(463)+0.7*A(464)+0.7*A(465)+0.7*A(466)+0.7*A(467) Vdot(84) = A(472)+A(473)+A(474) Vdot(85) = A(468) Vdot(86) = 0.5*A(493)+0.49*A(494)+0.5*A(497)+0.49*A(498)+A(500)+A(501)+A(502) Vdot(87) = 0.49*A(494)+0.49*A(498) Vdot(88) = A(505)+A(506) Vdot(89) = 0.03*A(494)+0.03*A(498)+0.3*A(507)+0.3*A(508)+A(509) Vdot(90) = A(514)+A(515)+A(516)+0.5*A(517)+A(520)+A(522)+A(529)+A(531)+A(532)+A(534)+A(535) Vdot(91) = 0.4*A(517)+A(518)+A(519) Vdot(92) = 0.7*A(516)+A(521)+A(530)+A(533) Vdot(93) = 0.259*A(491)+0.04*A(510)+0.75*A(649)+0.21*A(669)+0.54*A(670)+0.5*A(672) Vdot(94) = 0.52*A(581)+0.8*A(583)+A(584)+0.83*A(588)+0.378*A(590)+A(596)+0.65*A(598)+0.8*A(606)+0.52*A(608)+0.21& &*A(619)+0.507*A(627)+A(629)+0.86*A(630)+0.25*A(635)+0.459*A(661)+0.728*A(662)+0.441*A(702)+0.858*A(703)& &+0.033*A(706)+0.071*A(707)+0.042*A(719)+0.051*A(721)+0.638*A(743)+0.25*A(745)+0.5*A(755) Vdot(95) = 0.3*A(538)+0.3*A(563)+0.018*A(585)+0.83*A(588)+0.378*A(590)+0.249*A(719)+0.471*A(720)+0.303*A(721)+0.59& &*A(723) Vdot(96) = 0.311*A(490)+0.5*A(493)+0.12*A(495)+0.5*A(497)+0.12*A(499)+0.5*A(500)+0.5*A(501)+0.153*A(592)+0.07*A(594)& &+A(596)+0.65*A(598)+0.52*A(608)+0.52*A(611)+0.8*A(612)+0.36*A(622)+0.758*A(623)+0.459*A(661)+0.728*A(662)& &+0.551*A(698)+0.935*A(699)+0.242*A(706)+0.546*A(707)+0.34*A(719)+0.384*A(721)+0.187*A(736) Vdot(97) = 0.16*A(97)+0.8*A(606)+0.5*A(616)+0.5*A(647)+0.429*A(651)+0.25*A(737) Vdot(98) = 0.3*A(539)+0.3*A(565)+0.018*A(586)+0.153*A(592)+0.07*A(594)+0.036*A(625)+0.009*A(719)+0.054*A(720)+0.013& &*A(721)+0.07*A(723) Vdot(99) = A(646)+0.87*A(648) Vdot(100) = 0.551*A(698)+0.935*A(699)+0.242*A(706)+0.546*A(707)+0.589*A(719)+0.471*A(720)+0.687*A(721)+0.59*A(723)+0.2& &*A(727)+A(729)+0.187*A(736)+0.25*A(737)+0.33*A(751)+0.33*A(752)+0.35*A(756) Vdot(101) = 0.441*A(702)+0.858*A(703)+0.033*A(706)+0.071*A(707)+0.051*A(719)+0.054*A(720)+0.064*A(721)+0.07*A(723)& &+0.212*A(740)+0.25*A(741)+0.33*A(753)+0.33*A(754) Vdot(102) = 0.325*A(628)+0.65*A(632)+0.407*A(668)+A(681)+A(682) Vdot(103) = 0.75*A(635) Vdot(104) = 0.112*A(636)+0.182*A(637)+0.15*A(638)+0.15*A(639)+0.182*A(746)+0.112*A(747) Vdot(105) = A(686)+A(688) Vdot(106) = 0.5*A(718)& -3.2E-12*2.46E13*V(404)*V(106)-RCT(567)/4.6*V(106) ! modified by DDW to account for the loss of CHNP according to RCIM_plus ! RCT(567) is 3.0E-5 ! 2.46E13 converts the unit to ppm s-1 Vdot(107) = A(660)+A(664)+A(701)+A(705)+A(717)+0.5*A(718)+0.8*A(727)+0.25*A(736)+0.15*A(740)+0.15*A(743)& -1.0E-11*2.46E13*V(404)*V(107)-RCT(567)*V(107)-RCT(567)/2.0*V(107)& +2.0E-12*2.46E13*V(404)*V(108)+3.0E-12*2.46E13*V(404)*V(110)*0.67 ! modified by DDW to account for the loss of HNDP according to RCIM_plus ! RCT(567) is 3.0E-5 ! 2.46E13 converts the unit to ppm s-1 Vdot(108) = 0.234*A(719)+0.387*A(720)& -2.0E-12*2.46E13*V(404)*V(108)-RCT(567)/6.0*V(108)-RCT(567)/2.0*V(108) ! modified by DDW to account for the loss of DHDN according to RCIM_plus ! RCT(567) is 3.0E-5 ! 2.46E13 converts the unit to ppm s-1 Vdot(109) = A(663)+A(667)+A(700)+A(704)+A(708)& -2.0E-12*2.46E13*V(404)*V(109)-RCT(567)/3.0*V(109) ! modified by DDW to account for the loss of DHPN according to RCIM_plus ! RCT(567) is 3.0E-5 ! 2.46E13 converts the unit to ppm s-1 Vdot(110) = A(589)+A(593)+A(597)+0.482*A(661)+0.401*A(665)+0.418*A(698)+0.494*A(702)+0.621*A(706)& -3.0E-12*2.46E13*V(404)*V(110)-RCT(567)/6.0*V(110)-RCT(567)/4.6*V(110) ! modified by DDW to account for the loss of DHCN according to RCIM_plus ! RCT(567) is 3.0E-5 ! 2.46E13 converts the unit to ppm s-1 Vdot(111) = A(582)+A(607)+A(613)& -1.0E-11*2.46E13*V(404)*V(111)-RCT(567)/6.0*V(111)-RCT(567)*V(111) ! modified by DDW to account for the loss of CPDH according to RCIM_plus Vdot(112) = 0.35*A(581)+0.32*A(591)+0.35*A(608)+0.35*A(611)& -1.0E-11*2.46E13*V(404)*V(112)-RCT(567)/6.0*V(112)-RCT(567)*V(112) ! modified by DDW to account for the loss of DHDP according to RCIM_plus Vdot(113) = 0.59*A(590)+0.76*A(594)+0.35*A(598)& -3.0E-12*2.46E13*V(404)*V(113)-RCT(567)/2.0*V(113) ! modified by DDW to account for the loss of DHPE according to RCIM_plus Vdot(114) = A(587)+0.68*A(591)+A(595)& -3.0E-12*2.46E13*V(404)*V(114)-RCT(567)/4.6*V(114) ! modified by DDW to account for the loss of HPDN according to RCIM_plus Vdot(115) = A(722)+A(724)& -1.0E-12*2.46E13*V(404)*V(115)-RCT(567)/4.6*V(115)-RCT(567)/6.0*V(115) ! modified by DDW to account for the loss of CHDN according to RCIM_plus Vdot(116) = A(728)+A(738)+A(742)+A(744)& -4.0E-12*2.46E13*V(404)*V(116)-RCT(567)/3.0*V(116)-RCT(567)*V(116) ! modified by DDW to account for the loss of DCHP according to RCIM_plus Vdot(117) = A(604)+A(609)& -2.25E-11*2.46E13*V(404)*V(117)-RCT(567)*2.2*V(117) Vdot(118) = 0.212*A(743)+0.75*A(745) Vdot(119) = A(543)+0.5*A(545)+0.55*A(549)+0.5*A(559)+0.193*A(570)+0.125*A(576)+0.847*A(592)+0.17*A(594)+0.135*A(622)& &+0.054*A(719)+0.059*A(721) Vdot(120) = A(544)+0.5*A(545)+0.55*A(553)+0.5*A(561)+0.452*A(571)+0.165*A(574)+0.17*A(588)+0.032*A(590)+0.41*A(627)& &+0.004*A(719)+0.006*A(721) Vdot(121) = 0.19*A(566)+0.2*A(622)+0.35*A(638) Vdot(122) = 0.438*A(624)+0.065*A(702)+0.142*A(703)+0.01*A(706)+0.028*A(707)+0.01*A(720)+0.026*A(723)+A(739)+0.638& &*A(740)+0.75*A(741) Vdot(123) = 0.562*A(624)+A(642)+0.599*A(665)+A(666) Vdot(124) = 0.37*A(628) Vdot(125) = 0.13*A(628) Vdot(126) = A(631)+A(645)+0.059*A(661)+0.272*A(662) Vdot(127) = 0.6*A(697)+0.008*A(719)+0.005*A(720)+0.014*A(721)+0.023*A(723)+0.281*A(732) Vdot(128) = 0.6*A(694)+0.06*A(719)+0.05*A(720)+0.17*A(721)+0.221*A(723)+0.66*A(731) Vdot(129) = 0.031*A(698)+0.065*A(699)+0.094*A(706)+0.355*A(707)+0.023*A(720)+0.07*A(723)+A(735)+0.563*A(736)+0.75& &*A(737) Vdot(130) = 0.09*A(568)+0.15*A(639) Vdot(131) = A(182) Vdot(132) = A(461)-A(487) Vdot(133) = -A(35) Vdot(134) = A(16)-A(17)-A(18) Vdot(135) = -A(52) Vdot(136) = -A(56) Vdot(137) = -A(148) Vdot(138) = A(211)-A(212) Vdot(139) = -A(275) Vdot(140) = -A(76) Vdot(141) = A(376)-A(377) Vdot(142) = A(414)-A(415) Vdot(143) = A(452)-A(453) Vdot(144) = A(121)-A(122) Vdot(145) = 0.16*A(574)-A(580) Vdot(146) = A(133)-A(134) Vdot(147) = A(141)-A(142) Vdot(148) = -A(185) Vdot(149) = -A(146)-A(147) Vdot(150) = -A(257) Vdot(151) = A(483)-A(484) Vdot(152) = 0.73*A(669)+0.4*A(670)-A(672) Vdot(153) = 0.16*A(619)-A(620) Vdot(154) = A(9)-A(10)-A(11) Vdot(155) = A(70)-A(75) Vdot(156) = -A(159) Vdot(157) = 0.74*A(159)-A(160)-A(161) Vdot(158) = A(180)-A(181) Vdot(159) = 0.84*A(185)-A(186)-A(187) Vdot(160) = A(195)-A(197) Vdot(161) = A(214)-A(215) Vdot(162) = A(239)-A(240) Vdot(163) = 0.74*A(257)-A(258)-A(259) Vdot(164) = A(271)-A(273) Vdot(165) = A(260)-A(274) Vdot(166) = A(303)-A(304) Vdot(167) = A(311)-A(312) Vdot(168) = A(320)-A(321) Vdot(169) = A(77)-A(85) Vdot(170) = A(81)-A(86) Vdot(171) = 0.399*A(676)+0.07*A(678)+0.355*A(679)-A(691) Vdot(172) = 0.599*A(726)-A(756) Vdot(173) = A(61)+A(66)-A(74) Vdot(174) = 0.36*A(351)+0.265*A(358)-A(421) Vdot(175) = 0.17*A(360)+0.142*A(361)+0.17*A(362)-A(424) Vdot(176) = A(496)-A(509) Vdot(177) = A(27)-A(28)-A(29) Vdot(178) = A(162)-A(184) Vdot(179) = A(200)-A(216) Vdot(180) = A(223)-A(235) Vdot(181) = A(224)-A(241) Vdot(182) = 0.74*A(242)-A(243)-A(244) Vdot(183) = A(245)-A(256) Vdot(184) = 0.3*A(322)-A(324)-A(327) Vdot(185) = 0.3*A(294)-A(325) Vdot(186) = 0.05*A(297)-A(326) Vdot(187) = 0.58*A(668)-A(669)-A(670)-A(671) Vdot(188) = A(513)+0.1*A(517)-A(535) Vdot(189) = 0.473*A(674)-A(709)-A(711)-A(714) Vdot(190) = 0.4*A(614)-A(615)-A(616)+0.25*A(649) Vdot(191) = A(19)-A(20)+A(21)-A(39) Vdot(192) = -A(145)+0.034*A(679)+0.355*A(680) Vdot(193) = 0.175*A(357)+0.175*A(358)+0.175*A(359)-A(423) Vdot(194) = 0.36*A(354)+0.28*A(361)-A(422) Vdot(195) = 0.74*A(196)-A(198)-A(199) Vdot(196) = A(197)+A(199)-A(211) Vdot(197) = 0.88*A(499)-A(504)-A(506)-A(508) Vdot(198) = -A(523)-A(524)-A(525) Vdot(199) = -A(526)-A(527)-A(528) Vdot(200) = A(546)-A(656)-A(657) Vdot(201) = A(550)-A(658)-A(659) Vdot(202) = A(633)-A(634)-A(635) Vdot(203) = -A(65)+0.1*A(112)+0.5*A(113)+0.16*A(335) Vdot(204) = 0.3*A(205)+0.3*A(206)+0.3*A(207)-A(214) Vdot(205) = A(265)+A(266)+A(267)-A(271) Vdot(206) = A(268)+A(269)+A(270)-A(272) Vdot(207) = A(605)+0.2*A(606)+0.13*A(608)-A(619) Vdot(208) = 0.455*A(637)-A(646)-A(648)+0.026*A(725)+0.455*A(746) Vdot(209) = 0.16*A(159)+0.16*A(185)-A(195)-0.84*A(196) Vdot(210) = 0.474*A(676)+0.34*A(677)+0.244*A(680)+0.841*A(682)+0.841*A(684)+0.841*A(687)-A(689)-A(690) Vdot(211) = 0.3*A(167)+0.3*A(168)+0.3*A(169)-A(173)+0.5*A(189)+0.5*A(190)+0.5*A(191) Vdot(212) = A(170)+A(171)+A(172)-A(174)-A(175)+A(192)+A(193)+A(194)+A(208)+A(209)+A(210)+A(232)+A(233)+A(234)+A(253)& &+A(254)+A(255) Vdot(213) = A(201)+A(202)+A(203)-A(213)+A(216)+A(218) Vdot(214) = A(220)-A(222)-A(223) Vdot(215) = A(236)-A(238)+A(246)+A(247)+A(248)+A(256) Vdot(216) = A(225)+A(226)+A(227)-A(236)+A(241) Vdot(217) = A(237)-A(239)+0.3*A(250)+0.3*A(251)+0.3*A(252) Vdot(218) = 0.3*A(229)+0.3*A(230)+0.3*A(231)-A(237) Vdot(219) = A(217)+A(219)-A(242) Vdot(220) = 0.15*A(538)+0.15*A(563)-A(566)-A(567)+0.009*A(585) Vdot(221) = 0.15*A(539)+0.15*A(565)-A(568)-A(569)+0.009*A(586) Vdot(222) = 0.655*A(572)+0.655*A(577)-A(595)-A(596)-A(597)-A(598) Vdot(223) = 0.255*A(585)+0.595*A(586)-A(599)+A(692)+A(695) Vdot(224) = 0.315*A(638)-A(652)-A(653) Vdot(225) = 0.315*A(639)-A(654)-A(655) Vdot(226) = A(656)-A(660)-A(661)-A(662)-A(663) Vdot(227) = A(658)-A(664)-A(665)-A(666)-A(667) Vdot(228) = 0.038*A(676)+0.07*A(678)+0.034*A(679)-A(734)-A(748) Vdot(229) = A(693)-A(698)-A(699)-A(700)-A(701) Vdot(230) = A(696)-A(702)-A(703)-A(704)-A(705) Vdot(231) = A(709)-A(718)-A(719)-A(721)-A(722) Vdot(232) = A(710)-A(717)-A(720)-A(723)-A(724) Vdot(233) = A(749)-0.33*A(751)-0.33*A(752) Vdot(234) = A(750)-0.33*A(753)-0.33*A(754) Vdot(235) = 0.552*A(636)-A(647)-A(649)+0.187*A(725)+0.552*A(747) Vdot(236) = 0.88*A(495)+A(502)-A(503)-A(505)-A(507) Vdot(237) = A(56)-A(57)-A(58)-A(59) Vdot(238) = A(457)-A(461)-A(462)-A(463)-A(464) Vdot(239) = A(458)-A(465)-A(466)-A(467) Vdot(240) = A(512)-A(520)-A(521)-A(522) Vdot(241) = 0.1*A(257)-A(260)-A(261)-A(262)-A(263) Vdot(242) = A(78)+A(79)+A(80)-A(81)-A(82)-A(83)-A(84) Vdot(243) = A(31)+A(32)-A(36)-A(37)+0.15*A(281)+0.35*A(335)+0.069*A(349)+0.09*A(382)+0.055*A(396)+0.055*A(436)+0.055& &*A(459)+0.013*A(668)+0.06*A(669)+0.06*A(670) Vdot(244) = 0.201*A(725)+0.201*A(726)-A(727)-A(728)-A(729) Vdot(245) = A(610)+0.13*A(611)+0.2*A(612)-A(614) Vdot(246) = -A(60)+0.1*A(102)+0.5*A(103)+0.2*A(276)+0.2*A(277)+0.2*A(278) Vdot(247) = 0.1*A(159)-A(162)-A(163)-A(164)-A(165) Vdot(248) = A(166)-A(170)-A(171)-A(172) Vdot(249) = A(188)-A(192)-A(193)-A(194) Vdot(250) = A(204)-A(208)-A(209)-A(210) Vdot(251) = A(228)-A(232)-A(233)-A(234) Vdot(252) = A(249)-A(253)-A(254)-A(255) Vdot(253) = 0.182*A(637)-A(643)-A(644)-A(645)+0.182*A(746) Vdot(254) = A(390)+0.7*A(391)+A(392)-0.23*A(419) Vdot(255) = 0.388*A(725)+0.2*A(726)+0.34*A(731)-A(735)-A(736)-A(737)-A(738) Vdot(256) = A(407)+0.6*A(408)+A(409)-0.23*A(420) Vdot(257) = 0.112*A(636)+0.182*A(637)-A(650)-A(651)+0.182*A(746)+0.112*A(747) Vdot(258) = A(394)-A(398)-A(399)-A(400) Vdot(259) = 0.345*A(572)+0.255*A(586)-A(591)-A(592)-A(593)-A(594) Vdot(260) = A(366)+0.7*A(367)+A(368)-A(375)-A(376)+A(377)-A(378)-A(379) Vdot(261) = A(88)-A(94)-A(95)-A(96) Vdot(262) = A(510)-A(513)-A(514)-A(515)-A(516) Vdot(263) = 0.224*A(636)-A(640)-A(641)-A(642)+0.224*A(747) Vdot(264) = A(410)+0.7*A(411)+A(412)-A(413)-A(414)+A(415)-A(416)-A(417) Vdot(265) = 0.6*A(347)-A(351)-A(352)-A(353) Vdot(266) = A(434)-A(438)-A(439)-A(440)-A(441) Vdot(267) = A(435)-A(442)-A(443)-A(444) Vdot(268) = 0.345*A(577)+0.595*A(585)-A(587)-A(588)-A(589)-A(590) Vdot(269) = A(448)+A(449)-A(451)-A(452)+A(453)-A(454)-A(455) Vdot(270) = 2*A(11)+A(21)+A(22)-A(23)+0.2*A(33)+A(118)+A(128)+A(131)+A(139)+A(177)+A(195)+A(211)+A(217)+A(301)+A(309)& &+A(449)+A(479) Vdot(271) = A(469)+A(470)+A(471)-A(475)-A(476)-A(477) Vdot(272) = A(264)-A(268)-A(269)-A(270) Vdot(273) = A(76)-A(77)-A(78)-A(79)-A(80) Vdot(274) = 0.026*A(725)+0.719*A(732)-A(739)-A(740)-A(741)-A(742) Vdot(275) = 2*A(540)+A(542)+A(547)+0.063*A(554)+A(558)+A(562)+0.23*A(566)+A(567)+0.125*A(576)+A(579)-A(621)+0.035& &*A(638)+A(652)+A(653)+0.177*A(668)+0.479*A(674)+0.544*A(676)+1.737*A(678)+0.583*A(679)+0.903*A(681)+0.903& &*A(683)+0.904*A(685) Vdot(276) = 0.68*A(573)+0.67*A(578)-A(600)-A(601)+0.67*A(657)+0.67*A(659) Vdot(277) = 0.13*A(581)+0.2*A(583)+0.6*A(614)-A(617)-A(618)+0.255*A(622) Vdot(278) = A(333)-A(337)-A(338)-A(339) Vdot(279) = 0.25*A(282)+0.51*A(335)+0.25*A(336)+0.65*A(337)+0.7*A(338)+A(339)+0.5*A(340)+0.318*A(341)+0.83*A(342)& &-A(343) Vdot(280) = 0.32*A(573)+0.33*A(578)-A(602)-A(603)+0.33*A(657)+0.33*A(659) Vdot(281) = A(67)+A(68)+A(69)-A(70)-A(71)-A(72)-A(73)+A(110) Vdot(282) = A(380)-A(384)-A(385)-A(386) Vdot(283) = A(87)-A(91)-A(92)-A(93)+A(147) Vdot(284) = A(395)-A(401)-A(402)-A(403) Vdot(285) = A(675)-A(710)-A(712)-A(713)-A(715)-A(716) Vdot(286) = 0.445*A(459)-A(472)-A(473)-A(474) Vdot(287) = 0.1*A(196)-A(200)-A(201)-A(202)-A(203) Vdot(288) = A(100)-A(104)-A(105)-A(106)+A(148) Vdot(289) = A(221)-A(224)-A(225)-A(226)-A(227) Vdot(290) = A(511)-A(517)-A(518)-A(519) Vdot(291) = A(523)+A(524)+A(525)-A(529)-A(530)-A(531) Vdot(292) = A(526)+A(527)+A(528)-A(532)-A(533)-A(534) Vdot(293) = A(101)-A(107)-A(108)-A(109) Vdot(294) = A(543)+0.5*A(545)+0.45*A(549)+0.5*A(559)-A(570)+0.06*A(586) Vdot(295) = A(111)-A(114)-A(115)-A(116) Vdot(296) = 0.4*A(347)-A(354)-A(355)-A(356) Vdot(297) = 0.445*A(436)-A(445)-A(446)-A(447) Vdot(298) = A(544)+0.5*A(545)+0.45*A(553)+0.5*A(561)-A(571)+0.06*A(585) Vdot(299) = 0.937*A(554)-A(576)-A(577)-A(578)-A(579) Vdot(300) = -A(510)-A(511)-A(512) Vdot(301) = A(293)+0.7*A(294)+A(295)+A(307)+A(310)+0.8*A(313)-A(315)-A(316)-A(317)-A(318)+A(319)+0.6*A(322) Vdot(302) = A(149)-A(151)-A(152)-A(153) Vdot(303) = A(154)-A(156)-A(157)-A(158) Vdot(304) = 0.445*A(459)-A(468)-A(469)-A(470)-A(471) Vdot(305) = 0.222*A(396)-A(410)-A(411)-A(412) Vdot(306) = A(381)-A(387)-A(388)-A(389) Vdot(307) = 0.065*A(490)-A(496)-A(497)-A(498)-A(499)+0.7*A(508) Vdot(308) = A(492)-A(500)-A(501)-A(502) Vdot(309) = 0.4*A(348)-A(360)-A(361)-A(362) Vdot(310) = 0.2*A(281)-A(289)-A(290)-A(291)-A(292)+A(299)+A(302)+0.8*A(305) Vdot(311) = 0.164*A(349)-A(372)-A(373)-A(374) Vdot(312) = 0.133*A(349)-A(366)-A(367)-A(368) Vdot(313) = A(279)-A(283)-A(284)-A(285) Vdot(314) = 0.937*A(556)-A(572)-A(573)-A(574)-A(575) Vdot(315) = 0.445*A(396)-A(404)-A(405)-A(406)+0.875*A(413)+0.4375*A(416) Vdot(316) = 0.1*A(242)-A(245)-A(246)-A(247)-A(248) Vdot(317) = A(711)+A(712)-A(725)-A(726) Vdot(318) = 0.172*A(725)+A(734)-A(743)-A(744)-A(745) Vdot(319) = A(557)-A(586) Vdot(320) = A(52)-A(53)-A(54)-A(55)+0.6*A(90)+A(129)+A(132)+A(136) Vdot(321) = A(160)-A(166)-A(167)-A(168)-A(169) Vdot(322) = 0.41*A(382)-A(390)-A(391)-A(392) Vdot(323) = A(186)-A(188)-A(189)-A(190)-A(191) Vdot(324) = 0.25*A(538)+0.25*A(563)+0.015*A(585)-A(636)-A(638) Vdot(325) = 0.58*A(566)+0.77*A(568)+0.75*A(599)+A(600)+A(602)+0.075*A(638)+0.175*A(639)+A(713)+A(730)+A(733)-A(755) Vdot(326) = A(555)-A(585) Vdot(327) = 0.193*A(571)+0.33*A(601)+0.19*A(603)-A(604)-A(605)-A(606)-A(607)-A(608) Vdot(328) = 0.59*A(382)+0.64*A(384)+0.7*A(385)+A(386)+0.83*A(387)+0.578*A(388)+0.83*A(389)+0.055*A(396)+0.65*A(398)& &+0.7*A(399)+A(400)+0.825*A(401)+0.56*A(402)+0.825*A(403)-0.2*A(418) Vdot(329) = 0.25*A(539)+0.25*A(565)+0.015*A(586)-A(637)-A(639) Vdot(330) = A(258)-A(264)-A(265)-A(266)-A(267) Vdot(331) = A(275)-A(276)-A(277)-A(278) Vdot(332) = 0.6*A(348)-A(357)-A(358)-A(359) Vdot(333) = A(690)+A(691)-A(706)-A(707)-A(708) Vdot(334) = A(280)-A(286)-A(287)-A(288) Vdot(335) = A(176)+A(177)-A(179)-A(180)+A(181)-A(182)-A(183) Vdot(336) = A(117)+A(118)-A(120)-A(121)+A(122)-A(123)-A(124) Vdot(337) = 0.4*A(90)-A(97)-A(98)-A(99)+0.1*A(103)+0.1*A(113) Vdot(338) = A(478)+A(479)+A(481)-A(482)-A(483)+A(484)-A(485)-A(486) Vdot(339) = A(138)+A(139)-A(140)-A(141)+A(142)-A(143)-A(144) Vdot(340) = -A(490)-A(491)-A(492) Vdot(341) = 0.134*A(349)-A(363)-A(364)-A(365)+0.79*A(369)+0.9*A(370)+0.9*A(371) Vdot(342) = A(334)-A(340)-A(341)-A(342) Vdot(343) = 0.27*A(281)-A(296)-A(297)-A(298) Vdot(344) = -A(46)-A(47)-A(48)+A(69)+A(80)+A(96)+A(116)+A(153)+A(158)+0.825*A(288)+A(295)+A(298)+0.83*A(342)+A(346)& &+0.825*A(359)+0.83*A(362)+A(365)+A(368)+A(371)+A(374)+0.83*A(389)+A(392)+0.825*A(403)+A(406)+A(409)+A(412)& &+A(444)+A(447)+A(467)+A(471)+A(474)+A(477) Vdot(345) = -A(87)-A(88)-A(89)-A(90) Vdot(346) = 0.452*A(570)+0.67*A(601)+0.81*A(603)-A(609)-A(610)-A(611)-A(612)-A(613) Vdot(347) = A(198)-A(204)-A(205)-A(206)-A(207) Vdot(348) = 0.935*A(490)-A(493)-A(494)-A(495)+0.7*A(507) Vdot(349) = A(222)-A(228)-A(229)-A(230)-A(231) Vdot(350) = 0.223*A(396)+0.9*A(404)+0.53*A(405)+0.88*A(406)-A(407)-A(408)-A(409) Vdot(351) = 0.8*A(308)+A(309)-A(310)-A(311)+A(312)-A(313)-A(314) Vdot(352) = A(580)-A(581)-A(582)-A(583)-A(584)+0.25*A(599) Vdot(353) = 0.8*A(300)+A(301)-A(302)-A(303)+A(304)-A(305)-A(306) Vdot(354) = A(163)+A(164)+A(165)+A(184)-A(217)-A(218)-A(219)-A(220)-A(221) Vdot(355) = A(243)-A(249)-A(250)-A(251)-A(252) Vdot(356) = -A(536)-A(537)-A(538)-A(539)-A(668)-A(673) Vdot(357) = -A(434)-A(435)-A(436)-A(437) Vdot(358) = -A(110)-A(111)-A(112)-A(113) Vdot(359) = -A(333)-A(334)-A(335)-A(336) Vdot(360) = -A(457)-A(458)-A(459)-A(460) Vdot(361) = A(548)+0.006*A(676)+0.127*A(677)+0.057*A(680)-A(695)-A(696)-A(697) Vdot(362) = A(552)+0.032*A(676)+0.671*A(677)+0.298*A(680)-A(692)-A(693)-A(694) Vdot(363) = -A(149)-A(150)+0.2*A(276)+0.2*A(277)+0.2*A(278)+0.624*A(490)+0.015*A(491)+0.01*A(510) Vdot(364) = 2*A(541)+A(542)+A(551)+0.063*A(556)+A(560)+A(564)+0.14*A(568)+A(569)+0.165*A(574)+A(575)-A(625)-A(626)& &+0.035*A(639)+A(654)+A(655)+0.416*A(668)+0.048*A(674)+0.019*A(676)+0.123*A(678)+0.028*A(679)+0.097*A(681)& &+0.097*A(683)+0.096*A(685) Vdot(365) = 0.33*A(335)+A(343)-A(344)-A(345)-A(346) Vdot(366) = 0.055*A(436)+A(439)+A(440)+A(441)+A(442)+A(443)+A(444)-A(448)-A(449)-A(450) Vdot(367) = -A(100)-A(101)-A(102)-A(103) Vdot(368) = -A(347)-A(348)-A(349)-A(350) Vdot(369) = -A(43)-A(44)-A(45)+A(55)+A(59)+A(64)+A(73)+A(84)+A(93)+A(99)+A(106)+A(109)+A(165)+A(169)+A(172)+A(191)& &+A(194)+A(203)+A(207)+A(210)+A(227)+A(231)+A(234)+A(248)+A(252)+A(255)+A(263)+A(267)+A(270)+A(278)+A(285)& &+0.175*A(288)+A(292)+A(318)+A(339)+0.17*A(342)+A(353)+A(356)+0.175*A(359)+0.17*A(362)+A(386)+0.17*A(389)& &+A(400)+0.175*A(403)+A(441)+A(464)+A(515)+A(518)+A(521)+A(530)+A(533) Vdot(370) = A(65)-A(66)-A(67)-A(68)-A(69)+0.28*A(112)+0.1*A(113)+A(119)+A(120)+A(124) Vdot(371) = A(60)-A(61)-A(62)-A(63)-A(64)+0.28*A(102)+0.1*A(103)+A(137)+A(140)+A(144)+A(150)+A(155) Vdot(372) = -A(279)-A(280)-A(281)-A(282) Vdot(373) = -A(24)+0.63*A(89)+0.6*A(90)+A(97)+A(98)+A(99)+0.56*A(102)+0.56*A(112)+A(119)+A(125)+A(126)+A(127)+A(128)& &+A(129)+A(137)+A(178)+0.2*A(281)+A(299)+A(307)+0.03*A(335)+0.431*A(349)+0.41*A(382)+A(383)+0.445*A(396)& &+0.445*A(436)+A(450)+0.445*A(459)+A(480)+0.365*A(491)+0.6*A(538)+0.9*A(539)+A(543)+A(544)+A(545)+0.55& &*A(549)+0.55*A(553)+0.5*A(559)+0.5*A(561)+0.6*A(563)+0.9*A(565)+0.42*A(566)+A(567)+0.23*A(568)+A(569)+0.548& &*A(570)+0.807*A(571)+0.33*A(574)+0.25*A(576)+0.52*A(581)+0.8*A(583)+2*A(584)+0.036*A(585)+0.054*A(586)& &+A(605)+0.2*A(606)+0.13*A(608)+A(610)+0.13*A(611)+0.2*A(612)+0.5*A(616)+2*A(617)+A(618)+0.84*A(619)+A(620)& &+0.036*A(625)+0.507*A(627)+0.325*A(628)+A(629)+0.86*A(630)+0.65*A(632)+0.25*A(635)+0.888*A(636)+0.818& &*A(637)+0.46*A(638)+0.36*A(639)+A(646)+0.5*A(647)+A(648)+A(650)+0.571*A(651)+0.407*A(668)+0.6*A(694)+0.6& &*A(697)+0.2*A(727)+A(729)+0.66*A(731)+0.281*A(732)+A(735)+0.563*A(736)+0.75*A(737)+A(739)+0.638*A(740)+0.75& &*A(741)+0.818*A(746)+0.888*A(747)+0.33*A(751)+0.33*A(752)+1.5*A(755)+2*A(756) Vdot(374) = A(475)+A(476)+A(477)-A(478)-A(479)-A(480)-A(481) Vdot(375) = 0.5*A(112)+0.4*A(113)+A(114)+A(115)+A(116)-A(117)-A(118)-A(119) Vdot(376) = A(296)+0.8*A(297)+A(298)-A(319)-A(320)+A(321)-A(322)-A(323)+A(344)+0.35*A(345)+A(346) Vdot(377) = -A(558)-A(559)-A(560)-A(561)-A(679)-A(680) Vdot(378) = 0.236*A(349)+0.64*A(354)+0.7*A(355)+A(356)+0.83*A(360)+0.578*A(361)+0.83*A(362)+0.77*A(372)+0.77*A(373)& &+0.77*A(374)-0.23*A(393)-A(394)-A(395)-A(396)-A(397)+A(422)+A(424) Vdot(379) = 0.15*A(281)+0.6*A(283)+0.7*A(284)+0.8*A(285)+0.6*A(286)+0.795*A(287)+0.6*A(288)-A(299)-A(300)-A(301) Vdot(380) = 0.4*A(103)+A(104)+A(105)+A(106)+A(107)+A(108)+A(109)-A(137)-A(138)-A(139)+A(151)+A(152)+A(153)+A(156)& &+A(157)+A(158)+0.2*A(276)+0.2*A(277)+0.2*A(278)+0.311*A(490)+0.5*A(493) Vdot(381) = 0.152*A(676)+0.191*A(677)+0.102*A(680)+0.159*A(682)+0.159*A(684)+0.159*A(687)+0.4*A(697)+0.157*A(716)& &-A(732)-A(733)-A(747)-A(750) Vdot(382) = 0.267*A(349)-A(369)-A(370)-A(371)+0.88*A(375)+0.44*A(378)+A(383)+A(423) Vdot(383) = A(43)+A(46)+A(53)+A(54)+A(89)+1.61*A(91)+1.61*A(92)+2*A(94)+2*A(95)+A(97)+A(98)+0.5*A(102)+A(104)+A(105)& &+A(107)+A(108)+0.5*A(112)+A(114)+A(115)-A(125)-A(126)-A(127)-A(128)+A(145)+0.4*A(276)+0.4*A(277)+0.2*A(283)& &+0.225*A(286)+A(296)+0.8*A(297)+0.84*A(335)+0.65*A(337)+0.7*A(338)+0.5*A(340)+0.318*A(341)+0.2*A(349)+0.64& &*A(354)+0.7*A(355)+0.83*A(360)+0.578*A(361)+A(366)+0.7*A(367)+0.5*A(382)+0.64*A(384)+0.7*A(385)+0.83*A(387)& &+0.578*A(388)+A(410)+0.7*A(411)+A(422)+0.311*A(490)+0.5*A(493)+0.97*A(494)+0.12*A(495)+0.97*A(498)+0.12& &*A(499)+0.5*A(500)+A(503)+A(504)+0.76*A(510)+0.14*A(511)+0.3*A(538)+0.3*A(539)+2*A(540)+2*A(541)+2*A(542)& &+A(547)+A(551)+0.063*A(554)+0.063*A(556)+2*A(558)+A(559)+2*A(560)+A(561)+A(562)+0.3*A(563)+A(564)+0.3& &*A(565)+A(575)+A(579)+0.65*A(581)+A(583)+0.018*A(585)+0.018*A(586)+0.17*A(588)+0.032*A(590)+0.847*A(592)& &+0.17*A(594)+A(618)+0.05*A(622)+0.242*A(623)+0.083*A(627)+0.5*A(628)+0.14*A(630)+A(632)+0.075*A(638)+0.175& &*A(639)+0.059*A(661)+0.272*A(662)+0.599*A(665)+A(666)+0.827*A(668)+0.06*A(669)+0.06*A(670)+A(671)+0.5& &*A(672)+0.527*A(674)+0.563*A(676)+1.86*A(678)+1.577*A(679)+0.645*A(680)+A(681)+A(683)+A(685)+0.031*A(698)& &+0.065*A(699)+0.065*A(702)+0.142*A(703)+0.104*A(706)+0.383*A(707)+0.126*A(719)+0.088*A(720)+0.249*A(721)& &+0.34*A(723)+0.213*A(725)+0.2*A(727)+A(729)+0.212*A(743)+0.75*A(745)+A(748)+0.5*A(755) Vdot(384) = A(57)+A(58)+A(59)+A(62)+A(63)+A(64)+0.195*A(91)+0.195*A(92)+0.195*A(93)+0.5*A(102)-A(129)-A(130)-A(131)& &+A(146)+0.311*A(490)+0.01*A(510) Vdot(385) = 0.033*A(349)+0.64*A(351)+0.7*A(352)+A(353)+0.825*A(357)+0.56*A(358)+0.825*A(359)+0.8*A(363)+0.48*A(364)& &+0.8*A(365)-A(380)-A(381)-A(382)-A(383)+A(421) Vdot(386) = 0.336*A(676)+0.67*A(677)+0.299*A(680)+A(689)+0.4*A(694)+0.843*A(716)-A(730)-A(731)-A(746)-A(749) Vdot(387) = 0.964*A(625)-A(627)-A(629)-A(630)-A(631)+A(643)+A(644) Vdot(388) = A(290)+A(291)+A(292)-A(307)-A(308)-A(309) Vdot(389) = 0.33*A(281)-A(293)-A(294)-A(295)+0.2*A(300)+0.2*A(308) Vdot(390) = -A(154)-A(155)+0.2*A(283)+0.2*A(285)+0.225*A(286)+0.225*A(288)+A(316)+A(317)+A(318) Vdot(391) = A(1)-A(2)-A(3)-A(4)+A(14)+A(15)+A(18)-A(34)-A(38)-A(90)-A(103)-A(113)-A(282)-A(336)-A(350)-A(397)-A(437)& &-A(460) Vdot(392) = A(130)+A(131)-A(132)-A(133)+A(134)-A(135)-A(136)+A(150)+A(151)+A(152)+A(153)+A(155)+A(156)+A(157)+A(158)& &+A(178)+A(179)+A(183)+A(293)+0.7*A(294)+A(295)+A(363)+0.6*A(364)+A(365)+A(372)+A(373)+A(374)+A(390)+0.7& &*A(391)+A(392)+A(393)+A(407)+0.6*A(408)+A(409)+A(445)+A(446)+A(447)+A(472)+A(473)+A(474)+A(475)+A(476)& &+A(477) Vdot(393) = -A(49)-A(50)-A(51)+A(52)+A(56)+A(60)+A(65)+A(67)+A(68)+A(69)+A(76)+A(78)+A(79)+A(80)+A(87)+A(88)+A(90)& &+A(100)+A(101)+0.28*A(102)+0.2*A(103)+A(110)+A(111)+0.28*A(112)+0.2*A(113)+A(117)+A(118)+A(119)+A(120)& &+A(129)+A(130)+A(131)+A(132)+A(136)+A(137)+A(138)+A(139)+A(144)+A(147)+A(148)+A(149)+2*A(150)+A(151)+A(153)& &+A(154)+2*A(155)+A(156)+A(158)+0.06*A(159)+A(160)+A(166)+A(176)+A(177)+A(178)+A(179)+A(181)+A(183)+A(186)& &+A(188)+0.1*A(196)+A(198)+A(204)+A(221)+A(222)+A(228)+0.1*A(242)+A(243)+A(249)+0.1*A(257)+A(258)+A(264)& &+A(275)+A(279)+A(280)+0.8*A(281)+2*A(293)+1.4*A(294)+2*A(295)+A(296)+0.8*A(297)+A(298)+A(300)+A(301)+A(302)& &+A(304)+0.8*A(305)+A(307)+A(308)+A(309)+A(310)+A(312)+0.8*A(313)+A(319)+A(321)+0.6*A(322)+A(333)+A(334)& &+0.33*A(335)+A(343)+A(344)+0.35*A(345)+A(346)+A(347)+A(348)+0.698*A(349)+A(363)+0.6*A(364)+A(365)+A(366)& &+0.7*A(367)+A(368)+0.79*A(369)+0.9*A(370)+0.9*A(371)+A(372)+A(373)+A(374)+0.88*A(375)+0.44*A(378)+A(380)& &+A(381)+0.41*A(382)+A(383)+A(390)+0.7*A(391)+A(392)+A(393)+A(394)+A(395)+0.89*A(396)+0.9*A(404)+0.53*A(405)& &+0.88*A(406)+A(407)+0.6*A(408)+A(409)+A(410)+0.7*A(411)+A(412)+0.875*A(413)+0.4375*A(416)+A(423)+A(434)& &+A(435)+0.445*A(436)+A(446)+A(447)+A(448)+A(449)+A(457)+A(458)+0.89*A(459)+A(469)+A(470)+A(471)+A(472)& &+A(473)+A(474)+A(475)+A(476)+A(477)+A(478)+A(479)+A(480)+A(481)+A(482)+A(484)+A(486)+A(490)+A(492)-0.03& &*A(494)-0.03*A(498)+A(510)+A(511)+A(512)-A(519)+A(524)+A(525)+A(526)+A(527)+A(528)-A(531)-A(534) Vdot(394) = A(626)-A(628)-A(632)-A(633)+A(634)+A(748) Vdot(395) = A(621)-A(622)-A(623)-A(624)+A(640)+A(641) Vdot(396) = A(170)+A(171)+A(172)-A(176)-A(177)-A(178)+A(192)+A(193)+A(194)+A(208)+A(209)+A(210)+A(232)+A(233)+A(234)& &+A(253)+A(254)+A(255)+A(268)+A(269)+A(270)+0.5*A(501)+A(503)+A(504)+0.3*A(538)+0.3*A(539)+0.3*A(563)+0.3& &*A(565)+0.355*A(570)+0.355*A(571)+0.018*A(585)+0.018*A(586)+0.52*A(611)+0.8*A(612)+A(615)+0.5*A(616)+0.05& &*A(622)+0.242*A(623)+0.083*A(627)+0.14*A(630)+0.075*A(638)+0.175*A(639)+0.5*A(647)+0.75*A(649)+A(650)+0.571& &*A(651)+0.212*A(740)+0.25*A(741)+0.5*A(755)+0.65*A(756) Vdot(397) = -A(1)-A(3)-A(4)+A(5)-A(6)+2*A(7)+2*A(8)-A(9)+A(10)+A(14)+0.1*A(20)-2*A(21)-A(22)+A(26)-A(27)+A(28)+A(29)& &+0.8*A(33)+A(38)+A(39)+A(40)+2*A(41)+A(53)+A(57)+A(62)+A(67)+A(71)+A(74)+A(75)+A(78)+A(82)+A(85)+A(86)& &+A(91)+2*A(94)+A(95)+A(96)+A(97)+A(104)+2*A(107)+A(108)+A(109)+2*A(114)+A(115)+A(116)+A(120)-A(121)+A(122)& &+A(132)-A(133)+A(134)+A(140)-A(141)+A(142)+A(151)+A(156)-A(161)+A(163)+A(167)+A(170)+A(179)-A(180)+A(181)& &+A(184)-A(187)+A(189)+A(192)-A(197)-A(199)+A(201)+A(205)+A(208)-A(212)+A(216)-A(223)+A(225)+A(229)+A(232)& &+A(241)-A(244)+A(246)+A(250)+A(253)+A(256)-A(259)+A(261)+A(265)+A(268)+A(274)+A(276)+0.8*A(283)+1.825& &*A(286)+0.795*A(287)+0.825*A(288)+A(290)+A(293)+A(296)+A(302)-A(303)+A(304)+A(310)-A(311)+A(312)+A(316)& &+A(319)-A(320)+A(321)+0.65*A(337)+1.1*A(340)+0.058*A(341)+0.83*A(342)+A(344)+0.64*A(351)+0.64*A(354)+1.825& &*A(357)+0.56*A(358)+0.825*A(359)+1.83*A(360)+0.578*A(361)+0.83*A(362)+A(363)+A(366)+0.88*A(369)+A(372)& &+A(375)-A(376)+A(377)+0.64*A(384)+1.83*A(387)+0.578*A(388)+0.83*A(389)+A(390)+0.65*A(398)+1.825*A(401)+0.56& &*A(402)+0.825*A(403)+0.9*A(404)+A(407)+A(410)+A(413)-A(414)+A(415)+A(418)+A(419)+A(420)+A(421)+A(422)& &+A(423)+A(424)+A(439)+2*A(442)+A(443)+A(444)+A(445)+A(451)-A(452)+A(453)+A(462)+2*A(465)+A(466)+A(467)& &+A(469)+A(472)+A(475)+A(482)-A(483)+A(484)+A(487)+A(493)+A(497)+A(500)+A(501)+A(502)+A(509)+A(514)+A(517)& &+A(520)+A(529)+A(532)+A(535)+A(547)+A(549)+A(551)+A(553)+A(583)+A(588)+A(592)+A(596)+A(606)+A(612)+A(623)& &+A(630)+A(632)-A(633)+A(634)+A(640)-A(642)+A(643)-A(645)+A(652)+A(654)+A(657)+A(659)+0.459*A(661)+1.728& &*A(662)+A(666)+0.527*A(674)+0.563*A(676)+1.86*A(678)+0.611*A(679)+A(681)+2*A(683)+A(684)+2*A(685)+A(687)& &+A(692)+A(695)+A(699)+A(703)+A(707)+A(713)+0.058*A(719)+1.065*A(721)+A(723)+0.213*A(725)+A(729)+A(730)& &+A(733)+A(737)+A(741)+0.638*A(743)+1.25*A(745)+A(746)+A(747)+A(748)+A(751)+A(753)+0.65*A(756) Vdot(398) = A(537)-2*A(541)-A(542)-2*A(544)-A(545)-A(550)-A(551)-A(552)-A(553)-A(556)-A(557)-A(560)-A(561)-A(564)& &-A(565)+0.51*A(574)+0.03*A(586) Vdot(399) = A(480)+A(482)+A(486)+A(536)-2*A(540)-A(542)-2*A(543)-A(545)-A(546)-A(547)-A(548)-A(549)-A(554)-A(555)& &-A(558)-A(559)-A(562)-A(563)+0.75*A(576)+0.03*A(585) Vdot(400) = A(4)+A(6)-A(7)-A(9)+A(10)-A(12)-A(13)-A(14)+A(23)-A(33)-A(40)-2*A(41)-A(88)-A(101)-A(111)-A(118)-A(128)& &-A(131)-A(139)-A(177)-A(195)-A(211)-A(217)-A(280)-A(301)-A(309)-A(334)-A(348)-A(381)-A(395)-A(435)-A(449)& &-A(458)-A(479)-A(492)-A(512)-A(525)-A(528)+A(635)-A(673)-A(683)-A(684) Vdot(401) = A(2)-A(5)-A(6)-A(15)-A(16)-A(25)-A(30)-A(34)-A(89)-A(102)-A(112)+A(123)+A(135)+A(143)+A(173)+A(182)+A(213)& &+A(214)+A(215)+A(217)+A(219)+A(235)+A(236)+A(237)+A(238)+A(239)+A(240)+A(271)+A(272)+A(273)-A(281)+A(306)& &+A(314)+A(324)+A(325)+A(326)-A(335)-A(349)-A(382)-A(396)+0.1*A(418)+0.1*A(419)+0.1*A(420)-A(436)+A(454)& &-A(459)+A(485)-A(491)-A(511)-A(524)-A(527)+0.13*A(628)-A(668)-A(671) Vdot(402) = A(1)+A(3)-A(5)-A(7)-2*A(8)+A(12)+A(13)-A(19)+0.9*A(20)-A(26)-A(38)-A(53)-A(57)-A(61)-A(62)-A(66)-A(67)& &-A(70)-A(71)-A(77)-A(78)-A(81)-A(82)-A(91)-A(94)-A(97)-A(104)-A(107)-A(114)-A(120)-A(132)-A(140)-A(151)& &-A(156)-A(162)-A(163)-A(167)-A(170)-A(179)-A(189)-A(192)-A(200)-A(201)-A(205)-A(208)-A(224)-A(225)-A(229)& &-A(232)-A(245)-A(246)-A(250)-A(253)-A(260)-A(261)-A(265)-A(268)-A(276)-A(283)-A(286)-A(289)-A(290)-A(293)& &-A(296)-A(302)-A(310)-A(315)-A(316)-A(319)-A(337)-A(340)-A(344)-A(351)-A(354)-A(357)-A(360)-A(363)-A(366)& &-A(369)-A(372)-A(375)-A(384)-A(387)-A(390)-A(398)-A(401)-A(404)-A(407)-A(410)-A(413)-A(418)-A(419)-A(420)& &-A(438)-A(439)-A(442)-A(445)-A(451)-A(461)-A(462)-A(465)-A(468)-A(469)-A(472)-A(475)-A(482)-A(493)-A(496)& &-A(497)-A(500)-A(513)-A(514)-A(517)-A(520)-A(529)-A(532)-A(546)-A(547)-A(548)-A(549)-A(550)-A(551)-A(552)& &-A(553)-A(582)-A(583)-A(588)-A(589)-A(592)-A(593)-A(596)-A(597)-A(606)-A(607)-A(612)-A(613)-A(623)-A(624)& &-A(630)-A(631)-A(632)-A(640)-A(643)-A(652)-A(654)-A(662)-A(663)-A(666)-A(667)-A(685)-A(686)-A(687)-A(688)& &-A(699)-A(700)-A(703)-A(704)-A(707)-A(708)-A(721)-A(722)-A(723)-A(724)-A(728)-A(729)-A(737)-A(738)-A(741)& &-A(742)-A(744)-A(745)-A(751)-A(753) Vdot(403) = 0.3*A(538)+0.3*A(563)+0.018*A(585)+0.5*A(616)+A(617)+A(618)+0.63*A(619)+A(620)+0.36*A(622)+0.758*A(623)& &+0.175*A(628)+0.35*A(632)+0.5*A(647)+0.13*A(648)+0.429*A(651)-A(681)-A(682)+0.638*A(743)+0.25*A(745)+0.33& &*A(753)+0.33*A(754) Vdot(404) = 2*A(17)-A(19)+0.9*A(20)-A(22)-A(23)-A(24)-A(25)+A(26)-A(29)+A(30)+0.8*A(33)-A(35)+2*A(36)-A(37)-A(39)& &-A(40)-A(42)+A(43)-A(44)+A(46)-A(47)-A(52)-A(56)-A(60)-A(65)-A(74)-A(75)-A(76)-A(85)-A(86)-A(87)+0.12*A(89)& &-A(100)+0.36*A(102)-A(110)+0.36*A(112)-A(117)-A(127)-A(130)-A(138)-A(145)-A(146)-A(147)-A(148)-A(149)& &-A(154)-A(159)-A(173)-A(174)-A(176)-A(184)-A(185)-A(196)-A(213)-A(214)-A(215)-A(216)-A(218)-A(219)-A(220)& &-A(221)-A(235)-A(236)-A(237)-A(238)-A(239)-A(240)-A(241)-A(242)-A(256)-A(257)-A(271)-A(272)-A(273)-A(274)& &-A(275)-A(279)+0.8*A(281)-A(300)-A(308)-A(324)-A(325)-A(326)-A(333)+0.33*A(335)-A(343)-A(347)+0.698*A(349)& &-A(380)+0.41*A(382)-A(394)+0.89*A(396)-A(418)-A(419)-A(420)-A(434)+0.89*A(436)-A(448)-A(457)+0.89*A(459)& &-A(478)-A(490)+0.099*A(491)+0.12*A(495)+0.12*A(499)+A(503)+A(504)-0.7*A(507)-0.7*A(508)-A(509)-A(510)+0.06& &*A(511)-A(523)+0.06*A(524)-A(526)+0.06*A(527)-A(535)-A(536)-A(537)+0.5*A(538)+0.5*A(539)+A(543)+A(544)& &+A(545)+0.55*A(549)+0.55*A(553)+0.063*A(554)+0.063*A(556)+0.5*A(559)+0.5*A(561)+A(562)+1.5*A(563)+A(564)& &+1.5*A(565)+A(567)+A(569)+0.065*A(570)+0.065*A(571)-A(572)-0.84*A(574)+A(575)-A(576)-A(577)+A(579)-A(580)& &+0.65*A(581)+A(584)-0.595*A(585)-0.255*A(586)+A(587)+0.41*A(590)+A(591)+0.24*A(594)+A(595)+A(596)+1.3& &*A(598)-A(599)-A(600)-A(601)-A(602)-A(603)+A(605)+0.65*A(608)+A(610)+0.65*A(611)-A(614)-A(617)-0.37*A(619)& &-A(621)+0.665*A(622)-A(625)-A(626)+0.59*A(627)+0.5*A(628)+A(629)-A(635)+1.662*A(636)+1.637*A(637)+0.075& &*A(638)+0.175*A(639)+A(641)+A(644)+2*A(646)+A(647)-0.87*A(648)-0.25*A(649)+A(651)+A(652)+2*A(653)+A(654)+2& &*A(655)-A(656)-A(657)-A(658)-A(659)+0.518*A(661)+0.599*A(665)+0.28*A(668)-0.5*A(672)+0.527*A(674)-A(691)& &-A(692)-A(693)-0.4*A(694)-A(695)-A(696)-0.4*A(697)+0.582*A(698)+0.506*A(702)+0.379*A(706)-A(709)-A(710)& &-A(713)-A(714)-A(715)+A(717)+A(718)+1.147*A(719)+0.646*A(720)+0.435*A(721)+0.096*A(723)-A(725)-A(726)+0.2& &*A(727)-A(730)-A(731)-A(732)-A(733)-A(734)+A(735)+0.75*A(736)+A(739)+0.85*A(740)+0.85*A(743)+0.637*A(746)& &+0.662*A(747)-A(749)-A(750)+A(752)+A(754) Vdot(405) = 0.535*A(673)-A(675)-A(676)-2*A(677)-A(680)-A(682)-A(684)-A(687)-A(688)+A(715) Vdot(406) = 0.465*A(673)-A(674)-A(676)-2*A(678)-A(679)-A(681)-A(683)-A(685)-A(686)+A(714) Vdot(407) = 0.1*A(20)+A(24)+A(25)-A(26)-A(27)+A(28)-A(30)-2*A(31)-2*A(32)-A(33)+A(35)+A(37)+A(40)-A(42)+A(43)+A(44)& &+A(53)+A(54)-A(55)+A(57)+A(58)-A(59)+A(62)+A(63)-A(64)-A(69)+A(71)+A(72)-A(73)-A(80)+A(82)+A(83)-A(84)+0.12& &*A(89)+A(90)+A(91)+A(92)-A(93)-A(96)+A(97)+A(98)-A(99)+0.28*A(102)+A(104)+A(105)-A(106)+A(108)-A(109)+0.28& &*A(112)-A(116)+A(119)-A(123)+2*A(125)+A(127)+A(128)+A(129)-A(135)+A(137)-A(143)+A(145)+A(146)-A(153)-A(158)& &+0.16*A(159)+A(163)+A(164)-A(165)+A(167)+A(168)-A(169)+A(170)+A(171)-A(172)-A(173)+A(174)+A(178)-A(182)& &+0.16*A(185)+A(189)+A(190)-A(191)+A(192)+A(193)-A(194)+0.26*A(196)+A(201)+A(202)-A(203)+A(205)+A(206)& &-A(207)+A(208)+A(209)-A(210)-A(213)-A(214)-A(215)-A(217)+A(218)-A(219)+A(225)+A(226)-A(227)+A(229)+A(230)& &-A(231)+A(232)+A(233)-A(234)-A(235)-A(236)-A(237)-A(238)-A(239)-A(240)+0.16*A(242)+A(246)+A(247)-A(248)& &+A(250)+A(251)-A(252)+A(253)+A(254)-A(255)+0.16*A(257)+A(261)+A(262)-A(263)+A(265)+A(266)-A(267)+A(268)& &+A(269)-A(270)-A(271)-A(272)-A(273)+A(276)+A(277)-A(278)+0.8*A(283)+A(284)-A(285)+0.4*A(286)-0.775*A(288)& &+A(290)+A(291)-A(292)-A(295)-A(298)+A(299)-A(306)+A(307)-A(314)+A(316)+A(317)-A(318)-A(323)-A(324)-A(325)& &-A(326)+0.65*A(337)+0.7*A(338)-A(339)+0.1*A(340)+0.102*A(341)-A(342)-A(346)+0.64*A(351)+0.7*A(352)-A(353)& &+0.64*A(354)+0.7*A(355)-A(356)+0.175*A(357)+0.105*A(358)-A(359)+0.17*A(360)+0.102*A(361)-A(362)-A(365)& &-A(368)-A(371)-A(374)-A(379)+A(383)+0.64*A(384)+0.7*A(385)-A(386)+0.17*A(387)+0.102*A(388)-A(389)-A(392)& &+0.65*A(398)+0.7*A(399)-A(400)+0.175*A(401)+0.105*A(402)-A(403)-A(406)-A(409)-A(412)-A(417)+0.8*A(418)+0.8& &*A(419)+0.8*A(420)+0.445*A(436)+A(439)+A(440)-A(441)-A(444)-A(447)+2*A(450)+A(451)-A(454)+A(455)+A(462)& &+A(463)-A(464)-A(467)-A(471)-A(474)-A(477)+A(480)+A(481)-A(485)+0.084*A(491)+A(493)+0.97*A(494)-0.88*A(495)& &+0.97*A(498)-0.88*A(499)-A(502)+A(503)+A(504)+A(509)+A(514)-A(515)+0.7*A(516)+A(517)-A(518)+0.7*A(519)& &+A(520)-A(521)+0.7*A(522)+A(529)-A(530)+0.7*A(531)+A(532)-A(533)+0.7*A(534)+A(535)+0.4*A(538)+0.7*A(539)+2& &*A(540)+2*A(541)+2*A(542)+A(543)+A(544)+A(545)+A(547)+0.45*A(549)+A(551)+0.45*A(553)-0.937*A(554)-A(555)& &-0.937*A(556)-A(557)+2*A(558)+1.5*A(559)+2*A(560)+1.5*A(561)+0.4*A(563)+0.7*A(565)+A(567)+A(569)+0.193& &*A(570)+0.452*A(571)+0.33*A(574)+A(575)+0.25*A(576)+A(579)-0.35*A(581)+A(583)+A(584)+0.024*A(585)+0.042& &*A(586)+A(588)-0.59*A(590)+A(592)-0.76*A(594)-A(598)+0.75*A(599)+A(600)+A(602)+A(604)+A(606)-0.35*A(608)& &+A(609)-0.35*A(611)+A(612)+A(614)+1.5*A(616)+A(617)+A(618)+A(619)+A(620)-0.695*A(622)+0.242*A(623)+0.036& &*A(625)-0.493*A(627)-A(628)+0.86*A(630)+0.112*A(636)+0.182*A(637)+0.075*A(638)+0.175*A(639)-A(641)-A(644)& &+0.5*A(647)+0.87*A(648)+0.25*A(649)+0.571*A(651)-A(653)-A(655)+A(660)-0.941*A(661)+0.272*A(662)+A(664)& &-0.401*A(665)+A(666)+0.16*A(668)+0.5*A(672)-A(674)-A(675)+0.089*A(676)+0.064*A(677)+0.611*A(679)+0.336& &*A(680)+0.159*A(682)+0.159*A(684)+0.159*A(687)+A(689)+0.4*A(694)+0.4*A(697)-0.418*A(698)+A(699)+A(701)& &-0.494*A(702)+A(703)+A(705)-0.621*A(706)+A(707)-0.674*A(719)-0.42*A(720)+0.5*A(721)+0.904*A(723)-0.8*A(727)& &+A(729)+0.66*A(731)+0.281*A(732)-0.25*A(736)+A(737)-0.15*A(740)+A(741)-0.788*A(743)+0.75*A(745)+0.182& &*A(746)+0.112*A(747)+0.33*A(751)-0.67*A(752)-A(754)+0.65*A(756)

END SUBROUTINE Fun

! End of Fun function ! ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

END MODULE cacm3_Function