      subroutine lsrxn5(H2O,M,O2,CH4,H2,y,dbrk,r)
      implicit none
c
c----CAMx v7Beta6 190902
c
c     LSRXN5 computes double precision fluxes for each reaction
c
c     Copyright 1996 - 2018
c     Ramboll
c     Created by the CMC version 5.2.6
c
c     Routines Called:
c        none
c
c     Called by:
c        LSRATE5
c
      include "camx.prm"
      include "chmdat.inc"
      include "ddmchm.inc"
c
      double precision H2O, M, O2, CH4, H2, N2
      double precision y(*)
      double precision dbrk(*)
      double precision r(*)
c
c --- Entry point
c
      N2  = M - O2
      r(  1) = dbrk(  1)*y(lNO2)
      r(  2) = dbrk(  2)*y(lO3P)*O2*M
      r(  3) = dbrk(  3)*y(lO3P)*y(lO3)
      r(  4) = dbrk(  4)*y(lO3P)*y(lNO)
      r(  5) = dbrk(  5)*y(lO3P)*y(lNO2)
      r(  6) = dbrk(  6)*y(lO3P)*y(lNO2)
      r(  7) = dbrk(  7)*y(lO3)*y(lNO)
      r(  8) = dbrk(  8)*y(lO3)*y(lNO2)
      r(  9) = dbrk(  9)*y(lNO)*y(lNO3)
      r( 10) = dbrk( 10)*y(lNO)*y(lNO)*O2
      r( 11) = dbrk( 11)*y(lNO2)*y(lNO3)
      r( 12) = dbrk( 12)*y(lN2O5)
      r( 13) = dbrk( 13)*y(lN2O5)*H2O
      r( 14) = dbrk( 14)*y(lN2O5)*H2O*H2O
      r( 15) = dbrk( 15)*y(lNO2)*y(lNO3)
      r( 16) = dbrk( 16)*y(lNO3)
      r( 17) = dbrk( 17)*y(lNO3)
      r( 18) = dbrk( 18)*y(lO3)
      r( 19) = dbrk( 19)*y(lO3)
      r( 20) = dbrk( 20)*y(lO1D)*H2O
      r( 21) = dbrk( 21)*y(lO1D)*M
      r( 22) = dbrk( 22)*y(lOH)*y(lNO)
      r( 23) = dbrk( 23)*y(lHONO)
      r( 24) = dbrk( 24)*y(lOH)*y(lHONO)
      r( 25) = dbrk( 25)*y(lOH)*y(lNO2)
      r( 26) = dbrk( 26)*y(lOH)*y(lNO3)
      r( 27) = dbrk( 27)*y(lOH)*y(lHNO3)
      r( 28) = dbrk( 28)*y(lHNO3)
      r( 29) = dbrk( 29)*y(lOH)*y(lCO)
      r( 30) = dbrk( 30)*y(lOH)*y(lO3)
      r( 31) = dbrk( 31)*y(lHO2)*y(lNO)
      r( 32) = dbrk( 32)*y(lHO2)*y(lNO2)
      r( 33) = dbrk( 33)*y(lPNA)
      r( 34) = dbrk( 34)*y(lPNA)
      r( 35) = dbrk( 35)*y(lPNA)*y(lOH)
      r( 36) = dbrk( 36)*y(lHO2)*y(lO3)
      r( 37) = dbrk( 37)*y(lHO2)*y(lHO2)
      r( 38) = dbrk( 38)*y(lHO2)*y(lHO2)*H2O
      r( 39) = dbrk( 39)*y(lNO3)*y(lHO2)
      r( 40) = dbrk( 40)*y(lNO3)*y(lNO3)
      r( 41) = dbrk( 41)*y(lH2O2)
      r( 42) = dbrk( 42)*y(lH2O2)*y(lOH)
      r( 43) = dbrk( 43)*y(lOH)*y(lHO2)
      r( 44) = dbrk( 44)*y(lOH)*y(lSO2)
      r( 45) = dbrk( 45)*y(lOH)*H2
      r( 46) = dbrk( 46)*y(lMEO2)*y(lNO)
      r( 47) = dbrk( 47)*y(lMEO2)*y(lHO2)
      r( 48) = dbrk( 48)*y(lMEO2)*y(lHO2)
      r( 49) = dbrk( 49)*y(lMEO2)*y(lNO3)
      r( 50) = dbrk( 50)*y(lMEO2)*y(lMEO2)
      r( 51) = dbrk( 51)*y(lMEO2)*y(lMEO2)
      r( 52) = dbrk( 52)*y(lRO2C)*y(lNO)
      r( 53) = dbrk( 53)*y(lRO2C)*y(lHO2)
      r( 54) = dbrk( 54)*y(lRO2C)*y(lNO3)
      r( 55) = dbrk( 55)*y(lRO2C)*y(lMEO2)
      r( 56) = dbrk( 56)*y(lRO2C)*y(lRO2C)
      r( 57) = dbrk( 57)*y(lRO2X)*y(lNO)
      r( 58) = dbrk( 58)*y(lRO2X)*y(lHO2)
      r( 59) = dbrk( 59)*y(lRO2X)*y(lNO3)
      r( 60) = dbrk( 60)*y(lRO2X)*y(lMEO2)
      r( 61) = dbrk( 61)*y(lRO2X)*y(lRO2C)
      r( 62) = dbrk( 62)*y(lRO2X)*y(lRO2X)
      r( 63) = dbrk( 63)*y(lMCO3)*y(lNO2)
      r( 64) = dbrk( 64)*y(lPAN)
      r( 65) = dbrk( 65)*y(lPAN)
      r( 66) = dbrk( 66)*y(lMCO3)*y(lNO)
      r( 67) = dbrk( 67)*y(lMCO3)*y(lHO2)
      r( 68) = dbrk( 68)*y(lMCO3)*y(lNO3)
      r( 69) = dbrk( 69)*y(lMCO3)*y(lMEO2)
      r( 70) = dbrk( 70)*y(lMCO3)*y(lRO2C)
      r( 71) = dbrk( 71)*y(lMCO3)*y(lRO2X)
      r( 72) = dbrk( 72)*y(lMCO3)*y(lMCO3)
      r( 73) = dbrk( 73)*y(lRCO3)*y(lNO2)
      r( 74) = dbrk( 74)*y(lPAN2)
      r( 75) = dbrk( 75)*y(lPAN2)
      r( 76) = dbrk( 76)*y(lRCO3)*y(lNO)
      r( 77) = dbrk( 77)*y(lRCO3)*y(lHO2)
      r( 78) = dbrk( 78)*y(lRCO3)*y(lNO3)
      r( 79) = dbrk( 79)*y(lRCO3)*y(lMEO2)
      r( 80) = dbrk( 80)*y(lRCO3)*y(lRO2C)
      r( 81) = dbrk( 81)*y(lRCO3)*y(lRO2X)
      r( 82) = dbrk( 82)*y(lRCO3)*y(lMCO3)
      r( 83) = dbrk( 83)*y(lRCO3)*y(lRCO3)
      r( 84) = dbrk( 84)*y(lBZC3)*y(lNO2)
      r( 85) = dbrk( 85)*y(lPBZN)
      r( 86) = dbrk( 86)*y(lPBZN)
      r( 87) = dbrk( 87)*y(lBZC3)*y(lNO)
      r( 88) = dbrk( 88)*y(lBZC3)*y(lHO2)
      r( 89) = dbrk( 89)*y(lBZC3)*y(lNO3)
      r( 90) = dbrk( 90)*y(lBZC3)*y(lMEO2)
      r( 91) = dbrk( 91)*y(lBZC3)*y(lRO2C)
      r( 92) = dbrk( 92)*y(lBZC3)*y(lRO2X)
      r( 93) = dbrk( 93)*y(lBZC3)*y(lMCO3)
      r( 94) = dbrk( 94)*y(lBZC3)*y(lRCO3)
      r( 95) = dbrk( 95)*y(lBZC3)*y(lBZC3)
      r( 96) = dbrk( 96)*y(lMAC3)*y(lNO2)
      r( 97) = dbrk( 97)*y(lMPAN)
      r( 98) = dbrk( 98)*y(lMPAN)
      r( 99) = dbrk( 99)*y(lMAC3)*y(lNO)
      r(100) = dbrk(100)*y(lMAC3)*y(lHO2)
      r(101) = dbrk(101)*y(lMAC3)*y(lNO3)
      r(102) = dbrk(102)*y(lMAC3)*y(lMEO2)
      r(103) = dbrk(103)*y(lMAC3)*y(lRO2C)
      r(104) = dbrk(104)*y(lMAC3)*y(lRO2X)
      r(105) = dbrk(105)*y(lMAC3)*y(lMCO3)
      r(106) = dbrk(106)*y(lMAC3)*y(lRCO3)
      r(107) = dbrk(107)*y(lMAC3)*y(lBZC3)
      r(108) = dbrk(108)*y(lMAC3)*y(lMAC3)
      r(109) = dbrk(109)*y(lTBUO)*y(lNO2)
      r(110) = dbrk(110)*y(lTBUO)
      r(111) = dbrk(111)*y(lBZO)*y(lNO2)
      r(112) = dbrk(112)*y(lBZO)*y(lHO2)
      r(113) = dbrk(113)*y(lBZO)
      r(114) = dbrk(114)*y(lXHO2)*y(lNO)
      r(115) = dbrk(115)*y(lXHO2)*y(lHO2)
      r(116) = dbrk(116)*y(lXHO2)*y(lNO3)
      r(117) = dbrk(117)*y(lXHO2)*y(lMEO2)
      r(118) = dbrk(118)*y(lXHO2)*y(lRO2C)
      r(119) = dbrk(119)*y(lXHO2)*y(lRO2X)
      r(120) = dbrk(120)*y(lXHO2)*y(lMCO3)
      r(121) = dbrk(121)*y(lXHO2)*y(lRCO3)
      r(122) = dbrk(122)*y(lXHO2)*y(lBZC3)
      r(123) = dbrk(123)*y(lXHO2)*y(lMAC3)
      r(124) = dbrk(124)*y(lXOH)*y(lNO)
      r(125) = dbrk(125)*y(lXOH)*y(lHO2)
      r(126) = dbrk(126)*y(lXOH)*y(lNO3)
      r(127) = dbrk(127)*y(lXOH)*y(lMEO2)
      r(128) = dbrk(128)*y(lXOH)*y(lRO2C)
      r(129) = dbrk(129)*y(lXOH)*y(lRO2X)
      r(130) = dbrk(130)*y(lXOH)*y(lMCO3)
      r(131) = dbrk(131)*y(lXOH)*y(lRCO3)
      r(132) = dbrk(132)*y(lXOH)*y(lBZC3)
      r(133) = dbrk(133)*y(lXOH)*y(lMAC3)
      r(134) = dbrk(134)*y(lXNO2)*y(lNO)
      r(135) = dbrk(135)*y(lXNO2)*y(lHO2)
      r(136) = dbrk(136)*y(lXNO2)*y(lNO3)
      r(137) = dbrk(137)*y(lXNO2)*y(lMEO2)
      r(138) = dbrk(138)*y(lXNO2)*y(lRO2C)
      r(139) = dbrk(139)*y(lXNO2)*y(lRO2X)
      r(140) = dbrk(140)*y(lXNO2)*y(lMCO3)
      r(141) = dbrk(141)*y(lXNO2)*y(lRCO3)
      r(142) = dbrk(142)*y(lXNO2)*y(lBZC3)
      r(143) = dbrk(143)*y(lXNO2)*y(lMAC3)
      r(144) = dbrk(144)*y(lXMEO)*y(lNO)
      r(145) = dbrk(145)*y(lXMEO)*y(lHO2)
      r(146) = dbrk(146)*y(lXMEO)*y(lNO3)
      r(147) = dbrk(147)*y(lXMEO)*y(lMEO2)
      r(148) = dbrk(148)*y(lXMEO)*y(lRO2C)
      r(149) = dbrk(149)*y(lXMEO)*y(lRO2X)
      r(150) = dbrk(150)*y(lXMEO)*y(lMCO3)
      r(151) = dbrk(151)*y(lXMEO)*y(lRCO3)
      r(152) = dbrk(152)*y(lXMEO)*y(lBZC3)
      r(153) = dbrk(153)*y(lXMEO)*y(lMAC3)
      r(154) = dbrk(154)*y(lXMC3)*y(lNO)
      r(155) = dbrk(155)*y(lXMC3)*y(lHO2)
      r(156) = dbrk(156)*y(lXMC3)*y(lNO3)
      r(157) = dbrk(157)*y(lXMC3)*y(lMEO2)
      r(158) = dbrk(158)*y(lXMC3)*y(lRO2C)
      r(159) = dbrk(159)*y(lXMC3)*y(lRO2X)
      r(160) = dbrk(160)*y(lXMC3)*y(lMCO3)
      r(161) = dbrk(161)*y(lXMC3)*y(lRCO3)
      r(162) = dbrk(162)*y(lXMC3)*y(lBZC3)
      r(163) = dbrk(163)*y(lXMC3)*y(lMAC3)
      r(164) = dbrk(164)*y(lXRC3)*y(lNO)
      r(165) = dbrk(165)*y(lXRC3)*y(lHO2)
      r(166) = dbrk(166)*y(lXRC3)*y(lNO3)
      r(167) = dbrk(167)*y(lXRC3)*y(lMEO2)
      r(168) = dbrk(168)*y(lXRC3)*y(lRO2C)
      r(169) = dbrk(169)*y(lXRC3)*y(lRO2X)
      r(170) = dbrk(170)*y(lXRC3)*y(lMCO3)
      r(171) = dbrk(171)*y(lXRC3)*y(lRCO3)
      r(172) = dbrk(172)*y(lXRC3)*y(lBZC3)
      r(173) = dbrk(173)*y(lXRC3)*y(lMAC3)
      r(174) = dbrk(174)*y(lXMA3)*y(lNO)
      r(175) = dbrk(175)*y(lXMA3)*y(lHO2)
      r(176) = dbrk(176)*y(lXMA3)*y(lNO3)
      r(177) = dbrk(177)*y(lXMA3)*y(lMEO2)
      r(178) = dbrk(178)*y(lXMA3)*y(lRO2C)
      r(179) = dbrk(179)*y(lXMA3)*y(lRO2X)
      r(180) = dbrk(180)*y(lXMA3)*y(lMCO3)
      r(181) = dbrk(181)*y(lXMA3)*y(lRCO3)
      r(182) = dbrk(182)*y(lXMA3)*y(lBZC3)
      r(183) = dbrk(183)*y(lXMA3)*y(lMAC3)
      r(184) = dbrk(184)*y(lXTBU)*y(lNO)
      r(185) = dbrk(185)*y(lXTBU)*y(lHO2)
      r(186) = dbrk(186)*y(lXTBU)*y(lNO3)
      r(187) = dbrk(187)*y(lXTBU)*y(lMEO2)
      r(188) = dbrk(188)*y(lXTBU)*y(lRO2C)
      r(189) = dbrk(189)*y(lXTBU)*y(lRO2X)
      r(190) = dbrk(190)*y(lXTBU)*y(lMCO3)
      r(191) = dbrk(191)*y(lXTBU)*y(lRCO3)
      r(192) = dbrk(192)*y(lXTBU)*y(lBZC3)
      r(193) = dbrk(193)*y(lXTBU)*y(lMAC3)
      r(194) = dbrk(194)*y(lXCO)*y(lNO)
      r(195) = dbrk(195)*y(lXCO)*y(lHO2)
      r(196) = dbrk(196)*y(lXCO)*y(lNO3)
      r(197) = dbrk(197)*y(lXCO)*y(lMEO2)
      r(198) = dbrk(198)*y(lXCO)*y(lRO2C)
      r(199) = dbrk(199)*y(lXCO)*y(lRO2X)
      r(200) = dbrk(200)*y(lXCO)*y(lMCO3)
      r(201) = dbrk(201)*y(lXCO)*y(lRCO3)
      r(202) = dbrk(202)*y(lXCO)*y(lBZC3)
      r(203) = dbrk(203)*y(lXCO)*y(lMAC3)
      r(204) = dbrk(204)*y(lHCHO)
      r(205) = dbrk(205)*y(lHCHO)
      r(206) = dbrk(206)*y(lHCHO)*y(lOH)
      r(207) = dbrk(207)*y(lHCHO)*y(lNO3)
      r(208) = dbrk(208)*y(lCCHO)*y(lOH)
      r(209) = dbrk(209)*y(lCCHO)
      r(210) = dbrk(210)*y(lCCHO)*y(lNO3)
      r(211) = dbrk(211)*y(lRCHO)*y(lOH)
      r(212) = dbrk(212)*y(lRCHO)
      r(213) = dbrk(213)*y(lRCHO)*y(lNO3)
      r(214) = dbrk(214)*y(lACET)*y(lOH)
      r(215) = dbrk(215)*y(lACET)
      r(216) = dbrk(216)*y(lMEK)*y(lOH)
      r(217) = dbrk(217)*y(lMEK)
      r(218) = dbrk(218)*y(lMEOH)*y(lOH)
      r(219) = dbrk(219)*y(lFACD)*y(lOH)
      r(220) = dbrk(220)*y(lAACD)*y(lOH)
      r(221) = dbrk(221)*y(lPACD)*y(lOH)
      r(222) = dbrk(222)*y(lCOOH)*y(lOH)
      r(223) = dbrk(223)*y(lCOOH)
      r(224) = dbrk(224)*y(lROOH)*y(lOH)
      r(225) = dbrk(225)*y(lROOH)
      r(226) = dbrk(226)*y(lR6PX)*y(lOH)
      r(227) = dbrk(227)*y(lR6PX)
      r(228) = dbrk(228)*y(lRAPX)*y(lOH)
      r(229) = dbrk(229)*y(lRAPX)
      r(230) = dbrk(230)*y(lGLY)
      r(231) = dbrk(231)*y(lGLY)
      r(232) = dbrk(232)*y(lGLY)*y(lOH)
      r(233) = dbrk(233)*y(lGLY)*y(lNO3)
      r(234) = dbrk(234)*y(lMGLY)
      r(235) = dbrk(235)*y(lMGLY)*y(lOH)
      r(236) = dbrk(236)*y(lMGLY)*y(lNO3)
      r(237) = dbrk(237)*y(lBACL)
      r(238) = dbrk(238)*y(lCRES)*y(lOH)
      r(239) = dbrk(239)*y(lCRES)*y(lNO3)
      r(240) = dbrk(240)*y(lNPHE)*y(lOH)
      r(241) = dbrk(241)*y(lNPHE)
      r(242) = dbrk(242)*y(lNPHE)
      r(243) = dbrk(243)*y(lBALD)*y(lOH)
      r(244) = dbrk(244)*y(lBALD)
      r(245) = dbrk(245)*y(lBALD)*y(lNO3)
      r(246) = dbrk(246)*y(lAFG1)*y(lOH)
      r(247) = dbrk(247)*y(lAFG1)*y(lO3)
      r(248) = dbrk(248)*y(lAFG1)
      r(249) = dbrk(249)*y(lAFG2)*y(lOH)
      r(250) = dbrk(250)*y(lAFG2)*y(lO3)
      r(251) = dbrk(251)*y(lAFG2)
      r(252) = dbrk(252)*y(lAFG3)*y(lOH)
      r(253) = dbrk(253)*y(lAFG3)*y(lO3)
      r(254) = dbrk(254)*y(lMACR)*y(lOH)
      r(255) = dbrk(255)*y(lMACR)*y(lO3)
      r(256) = dbrk(256)*y(lMACR)*y(lNO3)
      r(257) = dbrk(257)*y(lMACR)*y(lO3P)
      r(258) = dbrk(258)*y(lMACR)
      r(259) = dbrk(259)*y(lMVK)*y(lOH)
      r(260) = dbrk(260)*y(lMVK)*y(lO3)
      r(261) = dbrk(261)*y(lMVK)*y(lO3P)
      r(262) = dbrk(262)*y(lMVK)
      r(263) = dbrk(263)*y(lIPRD)*y(lOH)
      r(264) = dbrk(264)*y(lIPRD)*y(lO3)
      r(265) = dbrk(265)*y(lIPRD)*y(lNO3)
      r(266) = dbrk(266)*y(lIPRD)
      r(267) = dbrk(267)*y(lPRD2)*y(lOH)
      r(268) = dbrk(268)*y(lPRD2)
      r(269) = dbrk(269)*y(lRNO3)*y(lOH)
      r(270) = dbrk(270)*y(lRNO3)
      r(271) = dbrk(271)*y(lGLYD)*y(lOH)
      r(272) = dbrk(272)*y(lGLYD)
      r(273) = dbrk(273)*y(lGLYD)*y(lNO3)
      r(274) = dbrk(274)*y(lACRO)*y(lOH)
      r(275) = dbrk(275)*y(lACRO)*y(lO3)
      r(276) = dbrk(276)*y(lACRO)*y(lNO3)
      r(277) = dbrk(277)*y(lACRO)*y(lO3P)
      r(278) = dbrk(278)*y(lACRO)
      r(279) = dbrk(279)*y(lCO3H)*y(lOH)
      r(280) = dbrk(280)*y(lCO3H)
      r(281) = dbrk(281)*y(lRO3H)*y(lOH)
      r(282) = dbrk(282)*y(lRO3H)
      r(283) = dbrk(283)*y(lXHCH)*y(lNO)
      r(284) = dbrk(284)*y(lXHCH)*y(lHO2)
      r(285) = dbrk(285)*y(lXHCH)*y(lNO3)
      r(286) = dbrk(286)*y(lXHCH)*y(lMEO2)
      r(287) = dbrk(287)*y(lXHCH)*y(lRO2C)
      r(288) = dbrk(288)*y(lXHCH)*y(lRO2X)
      r(289) = dbrk(289)*y(lXHCH)*y(lMCO3)
      r(290) = dbrk(290)*y(lXHCH)*y(lRCO3)
      r(291) = dbrk(291)*y(lXHCH)*y(lBZC3)
      r(292) = dbrk(292)*y(lXHCH)*y(lMAC3)
      r(293) = dbrk(293)*y(lXCCH)*y(lNO)
      r(294) = dbrk(294)*y(lXCCH)*y(lHO2)
      r(295) = dbrk(295)*y(lXCCH)*y(lNO3)
      r(296) = dbrk(296)*y(lXCCH)*y(lMEO2)
      r(297) = dbrk(297)*y(lXCCH)*y(lRO2C)
      r(298) = dbrk(298)*y(lXCCH)*y(lRO2X)
      r(299) = dbrk(299)*y(lXCCH)*y(lMCO3)
      r(300) = dbrk(300)*y(lXCCH)*y(lRCO3)
      r(301) = dbrk(301)*y(lXCCH)*y(lBZC3)
      r(302) = dbrk(302)*y(lXCCH)*y(lMAC3)
      r(303) = dbrk(303)*y(lXRCH)*y(lNO)
      r(304) = dbrk(304)*y(lXRCH)*y(lHO2)
      r(305) = dbrk(305)*y(lXRCH)*y(lNO3)
      r(306) = dbrk(306)*y(lXRCH)*y(lMEO2)
      r(307) = dbrk(307)*y(lXRCH)*y(lRO2C)
      r(308) = dbrk(308)*y(lXRCH)*y(lRO2X)
      r(309) = dbrk(309)*y(lXRCH)*y(lMCO3)
      r(310) = dbrk(310)*y(lXRCH)*y(lRCO3)
      r(311) = dbrk(311)*y(lXRCH)*y(lBZC3)
      r(312) = dbrk(312)*y(lXRCH)*y(lMAC3)
      r(313) = dbrk(313)*y(lXACE)*y(lNO)
      r(314) = dbrk(314)*y(lXACE)*y(lHO2)
      r(315) = dbrk(315)*y(lXACE)*y(lNO3)
      r(316) = dbrk(316)*y(lXACE)*y(lMEO2)
      r(317) = dbrk(317)*y(lXACE)*y(lRO2C)
      r(318) = dbrk(318)*y(lXACE)*y(lRO2X)
      r(319) = dbrk(319)*y(lXACE)*y(lMCO3)
      r(320) = dbrk(320)*y(lXACE)*y(lRCO3)
      r(321) = dbrk(321)*y(lXACE)*y(lBZC3)
      r(322) = dbrk(322)*y(lXACE)*y(lMAC3)
      r(323) = dbrk(323)*y(lXMEK)*y(lNO)
      r(324) = dbrk(324)*y(lXMEK)*y(lHO2)
      r(325) = dbrk(325)*y(lXMEK)*y(lNO3)
      r(326) = dbrk(326)*y(lXMEK)*y(lMEO2)
      r(327) = dbrk(327)*y(lXMEK)*y(lRO2C)
      r(328) = dbrk(328)*y(lXMEK)*y(lRO2X)
      r(329) = dbrk(329)*y(lXMEK)*y(lMCO3)
      r(330) = dbrk(330)*y(lXMEK)*y(lRCO3)
      r(331) = dbrk(331)*y(lXMEK)*y(lBZC3)
      r(332) = dbrk(332)*y(lXMEK)*y(lMAC3)
      r(333) = dbrk(333)*y(lXPD2)*y(lNO)
      r(334) = dbrk(334)*y(lXPD2)*y(lHO2)
      r(335) = dbrk(335)*y(lXPD2)*y(lNO3)
      r(336) = dbrk(336)*y(lXPD2)*y(lMEO2)
      r(337) = dbrk(337)*y(lXPD2)*y(lRO2C)
      r(338) = dbrk(338)*y(lXPD2)*y(lRO2X)
      r(339) = dbrk(339)*y(lXPD2)*y(lMCO3)
      r(340) = dbrk(340)*y(lXPD2)*y(lRCO3)
      r(341) = dbrk(341)*y(lXPD2)*y(lBZC3)
      r(342) = dbrk(342)*y(lXPD2)*y(lMAC3)
      r(343) = dbrk(343)*y(lXGLY)*y(lNO)
      r(344) = dbrk(344)*y(lXGLY)*y(lHO2)
      r(345) = dbrk(345)*y(lXGLY)*y(lNO3)
      r(346) = dbrk(346)*y(lXGLY)*y(lMEO2)
      r(347) = dbrk(347)*y(lXGLY)*y(lRO2C)
      r(348) = dbrk(348)*y(lXGLY)*y(lRO2X)
      r(349) = dbrk(349)*y(lXGLY)*y(lMCO3)
      r(350) = dbrk(350)*y(lXGLY)*y(lRCO3)
      r(351) = dbrk(351)*y(lXGLY)*y(lBZC3)
      r(352) = dbrk(352)*y(lXGLY)*y(lMAC3)
      r(353) = dbrk(353)*y(lXMGL)*y(lNO)
      r(354) = dbrk(354)*y(lXMGL)*y(lHO2)
      r(355) = dbrk(355)*y(lXMGL)*y(lNO3)
      r(356) = dbrk(356)*y(lXMGL)*y(lMEO2)
      r(357) = dbrk(357)*y(lXMGL)*y(lRO2C)
      r(358) = dbrk(358)*y(lXMGL)*y(lRO2X)
      r(359) = dbrk(359)*y(lXMGL)*y(lMCO3)
      r(360) = dbrk(360)*y(lXMGL)*y(lRCO3)
      r(361) = dbrk(361)*y(lXMGL)*y(lBZC3)
      r(362) = dbrk(362)*y(lXMGL)*y(lMAC3)
      r(363) = dbrk(363)*y(lXBAC)*y(lNO)
      r(364) = dbrk(364)*y(lXBAC)*y(lHO2)
      r(365) = dbrk(365)*y(lXBAC)*y(lNO3)
      r(366) = dbrk(366)*y(lXBAC)*y(lMEO2)
      r(367) = dbrk(367)*y(lXBAC)*y(lRO2C)
      r(368) = dbrk(368)*y(lXBAC)*y(lRO2X)
      r(369) = dbrk(369)*y(lXBAC)*y(lMCO3)
      r(370) = dbrk(370)*y(lXBAC)*y(lRCO3)
      r(371) = dbrk(371)*y(lXBAC)*y(lBZC3)
      r(372) = dbrk(372)*y(lXBAC)*y(lMAC3)
      r(373) = dbrk(373)*y(lXBAL)*y(lNO)
      r(374) = dbrk(374)*y(lXBAL)*y(lHO2)
      r(375) = dbrk(375)*y(lXBAL)*y(lNO3)
      r(376) = dbrk(376)*y(lXBAL)*y(lMEO2)
      r(377) = dbrk(377)*y(lXBAL)*y(lRO2C)
      r(378) = dbrk(378)*y(lXBAL)*y(lRO2X)
      r(379) = dbrk(379)*y(lXBAL)*y(lMCO3)
      r(380) = dbrk(380)*y(lXBAL)*y(lRCO3)
      r(381) = dbrk(381)*y(lXBAL)*y(lBZC3)
      r(382) = dbrk(382)*y(lXBAL)*y(lMAC3)
      r(383) = dbrk(383)*y(lXAF1)*y(lNO)
      r(384) = dbrk(384)*y(lXAF1)*y(lHO2)
      r(385) = dbrk(385)*y(lXAF1)*y(lNO3)
      r(386) = dbrk(386)*y(lXAF1)*y(lMEO2)
      r(387) = dbrk(387)*y(lXAF1)*y(lRO2C)
      r(388) = dbrk(388)*y(lXAF1)*y(lRO2X)
      r(389) = dbrk(389)*y(lXAF1)*y(lMCO3)
      r(390) = dbrk(390)*y(lXAF1)*y(lRCO3)
      r(391) = dbrk(391)*y(lXAF1)*y(lBZC3)
      r(392) = dbrk(392)*y(lXAF1)*y(lMAC3)
      r(393) = dbrk(393)*y(lXAF2)*y(lNO)
      r(394) = dbrk(394)*y(lXAF2)*y(lHO2)
      r(395) = dbrk(395)*y(lXAF2)*y(lNO3)
      r(396) = dbrk(396)*y(lXAF2)*y(lMEO2)
      r(397) = dbrk(397)*y(lXAF2)*y(lRO2C)
      r(398) = dbrk(398)*y(lXAF2)*y(lRO2X)
      r(399) = dbrk(399)*y(lXAF2)*y(lMCO3)
      r(400) = dbrk(400)*y(lXAF2)*y(lRCO3)
      r(401) = dbrk(401)*y(lXAF2)*y(lBZC3)
      r(402) = dbrk(402)*y(lXAF2)*y(lMAC3)
      r(403) = dbrk(403)*y(lXAF3)*y(lNO)
      r(404) = dbrk(404)*y(lXAF3)*y(lHO2)
      r(405) = dbrk(405)*y(lXAF3)*y(lNO3)
      r(406) = dbrk(406)*y(lXAF3)*y(lMEO2)
      r(407) = dbrk(407)*y(lXAF3)*y(lRO2C)
      r(408) = dbrk(408)*y(lXAF3)*y(lRO2X)
      r(409) = dbrk(409)*y(lXAF3)*y(lMCO3)
      r(410) = dbrk(410)*y(lXAF3)*y(lRCO3)
      r(411) = dbrk(411)*y(lXAF3)*y(lBZC3)
      r(412) = dbrk(412)*y(lXAF3)*y(lMAC3)
      r(413) = dbrk(413)*y(lXMAC)*y(lNO)
      r(414) = dbrk(414)*y(lXMAC)*y(lHO2)
      r(415) = dbrk(415)*y(lXMAC)*y(lNO3)
      r(416) = dbrk(416)*y(lXMAC)*y(lMEO2)
      r(417) = dbrk(417)*y(lXMAC)*y(lRO2C)
      r(418) = dbrk(418)*y(lXMAC)*y(lRO2X)
      r(419) = dbrk(419)*y(lXMAC)*y(lMCO3)
      r(420) = dbrk(420)*y(lXMAC)*y(lRCO3)
      r(421) = dbrk(421)*y(lXMAC)*y(lBZC3)
      r(422) = dbrk(422)*y(lXMAC)*y(lMAC3)
      r(423) = dbrk(423)*y(lXMVK)*y(lNO)
      r(424) = dbrk(424)*y(lXMVK)*y(lHO2)
      r(425) = dbrk(425)*y(lXMVK)*y(lNO3)
      r(426) = dbrk(426)*y(lXMVK)*y(lMEO2)
      r(427) = dbrk(427)*y(lXMVK)*y(lRO2C)
      r(428) = dbrk(428)*y(lXMVK)*y(lRO2X)
      r(429) = dbrk(429)*y(lXMVK)*y(lMCO3)
      r(430) = dbrk(430)*y(lXMVK)*y(lRCO3)
      r(431) = dbrk(431)*y(lXMVK)*y(lBZC3)
      r(432) = dbrk(432)*y(lXMVK)*y(lMAC3)
      r(433) = dbrk(433)*y(lXIPR)*y(lNO)
      r(434) = dbrk(434)*y(lXIPR)*y(lHO2)
      r(435) = dbrk(435)*y(lXIPR)*y(lNO3)
      r(436) = dbrk(436)*y(lXIPR)*y(lMEO2)
      r(437) = dbrk(437)*y(lXIPR)*y(lRO2C)
      r(438) = dbrk(438)*y(lXIPR)*y(lRO2X)
      r(439) = dbrk(439)*y(lXIPR)*y(lMCO3)
      r(440) = dbrk(440)*y(lXIPR)*y(lRCO3)
      r(441) = dbrk(441)*y(lXIPR)*y(lBZC3)
      r(442) = dbrk(442)*y(lXIPR)*y(lMAC3)
      r(443) = dbrk(443)*y(lXRN3)*y(lNO)
      r(444) = dbrk(444)*y(lXRN3)*y(lHO2)
      r(445) = dbrk(445)*y(lXRN3)*y(lNO3)
      r(446) = dbrk(446)*y(lXRN3)*y(lMEO2)
      r(447) = dbrk(447)*y(lXRN3)*y(lRO2C)
      r(448) = dbrk(448)*y(lXRN3)*y(lRO2X)
      r(449) = dbrk(449)*y(lXRN3)*y(lMCO3)
      r(450) = dbrk(450)*y(lXRN3)*y(lRCO3)
      r(451) = dbrk(451)*y(lXRN3)*y(lBZC3)
      r(452) = dbrk(452)*y(lXRN3)*y(lMAC3)
      r(453) = dbrk(453)*y(lYRPX)*y(lNO)
      r(454) = dbrk(454)*y(lYRPX)*y(lHO2)
      r(455) = dbrk(455)*y(lYRPX)*y(lNO3)
      r(456) = dbrk(456)*y(lYRPX)*y(lMEO2)
      r(457) = dbrk(457)*y(lYRPX)*y(lRO2C)
      r(458) = dbrk(458)*y(lYRPX)*y(lRO2X)
      r(459) = dbrk(459)*y(lYRPX)*y(lMCO3)
      r(460) = dbrk(460)*y(lYRPX)*y(lRCO3)
      r(461) = dbrk(461)*y(lYRPX)*y(lBZC3)
      r(462) = dbrk(462)*y(lYRPX)*y(lMAC3)
      r(463) = dbrk(463)*y(lY6PX)*y(lNO)
      r(464) = dbrk(464)*y(lY6PX)*y(lHO2)
      r(465) = dbrk(465)*y(lY6PX)*y(lNO3)
      r(466) = dbrk(466)*y(lY6PX)*y(lMEO2)
      r(467) = dbrk(467)*y(lY6PX)*y(lRO2C)
      r(468) = dbrk(468)*y(lY6PX)*y(lRO2X)
      r(469) = dbrk(469)*y(lY6PX)*y(lMCO3)
      r(470) = dbrk(470)*y(lY6PX)*y(lRCO3)
      r(471) = dbrk(471)*y(lY6PX)*y(lBZC3)
      r(472) = dbrk(472)*y(lY6PX)*y(lMAC3)
      r(473) = dbrk(473)*y(lYAPX)*y(lNO)
      r(474) = dbrk(474)*y(lYAPX)*y(lHO2)
      r(475) = dbrk(475)*y(lYAPX)*y(lNO3)
      r(476) = dbrk(476)*y(lYAPX)*y(lMEO2)
      r(477) = dbrk(477)*y(lYAPX)*y(lRO2C)
      r(478) = dbrk(478)*y(lYAPX)*y(lRO2X)
      r(479) = dbrk(479)*y(lYAPX)*y(lMCO3)
      r(480) = dbrk(480)*y(lYAPX)*y(lRCO3)
      r(481) = dbrk(481)*y(lYAPX)*y(lBZC3)
      r(482) = dbrk(482)*y(lYAPX)*y(lMAC3)
      r(483) = dbrk(483)*y(lZRN3)*y(lNO)
      r(484) = dbrk(484)*y(lZRN3)*y(lHO2)
      r(485) = dbrk(485)*y(lZRN3)*y(lNO3)
      r(486) = dbrk(486)*y(lZRN3)*y(lMEO2)
      r(487) = dbrk(487)*y(lZRN3)*y(lRO2C)
      r(488) = dbrk(488)*y(lZRN3)*y(lRO2X)
      r(489) = dbrk(489)*y(lZRN3)*y(lMCO3)
      r(490) = dbrk(490)*y(lZRN3)*y(lRCO3)
      r(491) = dbrk(491)*y(lZRN3)*y(lBZC3)
      r(492) = dbrk(492)*y(lZRN3)*y(lMAC3)
      r(493) = dbrk(493)*y(lXGLD)*y(lNO)
      r(494) = dbrk(494)*y(lXGLD)*y(lHO2)
      r(495) = dbrk(495)*y(lXGLD)*y(lNO3)
      r(496) = dbrk(496)*y(lXGLD)*y(lMEO2)
      r(497) = dbrk(497)*y(lXGLD)*y(lRO2C)
      r(498) = dbrk(498)*y(lXGLD)*y(lRO2X)
      r(499) = dbrk(499)*y(lXGLD)*y(lMCO3)
      r(500) = dbrk(500)*y(lXGLD)*y(lRCO3)
      r(501) = dbrk(501)*y(lXGLD)*y(lBZC3)
      r(502) = dbrk(502)*y(lXGLD)*y(lMAC3)
      r(503) = dbrk(503)*y(lXACR)*y(lNO)
      r(504) = dbrk(504)*y(lXACR)*y(lHO2)
      r(505) = dbrk(505)*y(lXACR)*y(lNO3)
      r(506) = dbrk(506)*y(lXACR)*y(lMEO2)
      r(507) = dbrk(507)*y(lXACR)*y(lRO2C)
      r(508) = dbrk(508)*y(lXACR)*y(lRO2X)
      r(509) = dbrk(509)*y(lXACR)*y(lMCO3)
      r(510) = dbrk(510)*y(lXACR)*y(lRCO3)
      r(511) = dbrk(511)*y(lXACR)*y(lBZC3)
      r(512) = dbrk(512)*y(lXACR)*y(lMAC3)
      r(513) = dbrk(513)*CH4*y(lOH)
      r(514) = dbrk(514)*y(lETHE)*y(lOH)
      r(515) = dbrk(515)*y(lETHE)*y(lO3)
      r(516) = dbrk(516)*y(lETHE)*y(lNO3)
      r(517) = dbrk(517)*y(lETHE)*y(lO3P)
      r(518) = dbrk(518)*y(lPRPE)*y(lOH)
      r(519) = dbrk(519)*y(lPRPE)*y(lO3)
      r(520) = dbrk(520)*y(lPRPE)*y(lNO3)
      r(521) = dbrk(521)*y(lPRPE)*y(lO3P)
      r(522) = dbrk(522)*y(lBD13)*y(lOH)
      r(523) = dbrk(523)*y(lBD13)*y(lO3)
      r(524) = dbrk(524)*y(lBD13)*y(lNO3)
      r(525) = dbrk(525)*y(lBD13)*y(lO3P)
      r(526) = dbrk(526)*y(lISOP)*y(lOH)
      r(527) = dbrk(527)*y(lISOP)*y(lO3)
      r(528) = dbrk(528)*y(lISOP)*y(lNO3)
      r(529) = dbrk(529)*y(lISOP)*y(lO3P)
      r(530) = dbrk(530)*y(lAPIN)*y(lOH)
      r(531) = dbrk(531)*y(lAPIN)*y(lO3)
      r(532) = dbrk(532)*y(lAPIN)*y(lNO3)
      r(533) = dbrk(533)*y(lAPIN)*y(lO3P)
      r(534) = dbrk(534)*y(lACYE)*y(lOH)
      r(535) = dbrk(535)*y(lACYE)*y(lO3)
      r(536) = dbrk(536)*y(lBENZ)*y(lOH)
      r(537) = dbrk(537)*y(lTOLU)*y(lOH)
      r(538) = dbrk(538)*y(lMXYL)*y(lOH)
      r(539) = dbrk(539)*y(lOXYL)*y(lOH)
      r(540) = dbrk(540)*y(lPXYL)*y(lOH)
      r(541) = dbrk(541)*y(lB124)*y(lOH)
      r(542) = dbrk(542)*y(lETOH)*y(lOH)
      r(543) = dbrk(543)*y(lALK1)*y(lOH)
      r(544) = dbrk(544)*y(lALK2)*y(lOH)
      r(545) = dbrk(545)*y(lALK3)*y(lOH)
      r(546) = dbrk(546)*y(lALK4)*y(lOH)
      r(547) = dbrk(547)*y(lALK5)*y(lOH)
      r(548) = dbrk(548)*y(lOLE1)*y(lOH)
      r(549) = dbrk(549)*y(lOLE1)*y(lO3)
      r(550) = dbrk(550)*y(lOLE1)*y(lNO3)
      r(551) = dbrk(551)*y(lOLE1)*y(lO3P)
      r(552) = dbrk(552)*y(lOLE2)*y(lOH)
      r(553) = dbrk(553)*y(lOLE2)*y(lO3)
      r(554) = dbrk(554)*y(lOLE2)*y(lNO3)
      r(555) = dbrk(555)*y(lOLE2)*y(lO3P)
      r(556) = dbrk(556)*y(lARO1)*y(lOH)
      r(557) = dbrk(557)*y(lARO2)*y(lOH)
      r(558) = dbrk(558)*y(lTERP)*y(lOH)
      r(559) = dbrk(559)*y(lTERP)*y(lO3)
      r(560) = dbrk(560)*y(lTERP)*y(lNO3)
      r(561) = dbrk(561)*y(lTERP)*y(lO3P)
      r(562) = dbrk(562)*y(lSESQ)*y(lOH)
      r(563) = dbrk(563)*y(lSESQ)*y(lO3)
      r(564) = dbrk(564)*y(lSESQ)*y(lNO3)
      r(565) = dbrk(565)*y(lSESQ)*y(lO3P)
c
      return
      end
