subroutine ebirate4(ny,nr,r,gain,loss) implicit none c c----CAMx v7Beta6 190902 c c EBIRATE4 computes species production and loss c for the EBI solver c c Copyright 1996 - 2019 c Ramboll c Created by the CMC version 5.2.6 c c --- Subroutines Called: c none c c --- Called by: c EBISOLV c c --- Argument definitions: c ny - dimension of gain and loss c nr - dimension of r c r - reaction rates (hr-1) c gain - species production (ppm/hr) c loss - species destruction (ppm/hr) c c --- Includes: include "camx.prm" include "chmdat.inc" include "ddmchm.inc" c c --- Arguments: integer ny, nr real loss(ny+1), gain(ny+1) real r(nr) c c --- Entry Point: c c c-----Calculate the species rates c c c H2O2 HNO3 SO2 SULF MEO2 RO2 PACD AACD CXO3 ALD2 c XO2H PANX FORM MEPX MEOH ROOH XO2 XO2N NTR1 NTR2 c FACD CO HCO3 ALDX GLYD GLY MGLY ETHA ETOH KET c PAR ACET PRPA XPRP XPAR ROR ETHY ETH OLE IOLE c ISOP ISO2 INTR ISPD ISPX HPLD OPO3 EPOX EPX2 TERP c BENZ CRES BZO2 OPEN TOL TO2 XOPN XYL XLO2 CRO c CAT1 CRON OPAN ECH4 I2 HOI I2O2 INO3 HIO3 IXOY c DMS c Loss(lH2O2 )= + r( 21)+ r( 22)+ r( 23) c Gain(lH2O2 )= + r( 17)+ r( 19)+ r( 20) & +( 0.040)*r(138)+( 0.080)*r(141) c Loss(lHNO3 )= + r( 46)+ r( 47) c Gain(lHNO3 )= +( 2.000)*r( 39)+ r( 45)+ r( 99) & + r(105)+ r(108)+ r(112) & + r(115)+ r(117)+( 0.717)*r(153) & + r(157)+ r(183)+ r(187) & + r(196)+ r(198)+ r(207) & + r(224)+ r(229)+ r(233) c Loss(lSO2 )= + r( 52)+ r(230) c Gain(lSO2 )= + r(231)+ r(233) c Loss(lSULF )= 0.0 c Gain(lSULF )= + r( 52)+ r(230)+ r(232) c Loss(lMEO2 )= + r( 71)+ r( 72)+ r( 73) & + r( 74) c Gain(lMEO2 )= + r( 53)+( 0.400)*r( 56)+( 0.440)*r( 57) & +( 2.000)*r( 59)+ r( 60)+( 0.900)*r( 73) & +( 0.800)*r( 77)+( 0.800)*r( 81)+( 0.800)*r( 85) & +( 0.600)*r( 87)+ r( 88)+ r( 94) & + r(106)+ r(121)+( 0.500)*r(125) & +( 1.380)*r(126)+( 0.800)*r(146)+( 0.115)*r(151) & +( 0.340)*r(154)+( 0.800)*r(161)+ r(169) & + r(174)+ r(180)+ r(203) & + r(208)+ r(231)+ r(232) Gain(lMEO2 ) = Gain(lMEO2 ) & + r(233) c Loss(lRO2 )= + r( 58)+ r( 66)+ r( 68) & + r( 69)+( 2.000)*r( 70)+ r( 74) & + r( 78)+ r( 82)+ r( 86) & + r(147)+ r(162)+ r(171) & + r(176)+ r(181)+ r(204) c Gain(lRO2 )= + r( 53)+( 0.400)*r( 56)+( 0.440)*r( 57) & +( 2.000)*r( 59)+( 2.000)*r( 60)+ r( 61) & +( 0.400)*r( 64)+( 0.440)*r( 65)+( 0.800)*r( 66) & +( 2.000)*r( 67)+( 0.900)*r( 73)+ r( 74) & +( 0.800)*r( 77)+ r( 78)+( 0.800)*r( 81) & + r( 82)+( 0.800)*r( 85)+ r( 86) & +( 0.600)*r( 87)+ r( 88)+( 0.600)*r( 89) & + r( 94)+ r(106)+ r(109) & +( 0.110)*r(111)+( 0.200)*r(113)+( 0.500)*r(115) Gain(lRO2 ) = Gain(lRO2 ) & + r(117)+ r(121)+ r(122) & +( 0.100)*r(124)+ r(125)+( 1.380)*r(126) & + r(127)+( 0.980)*r(130)+ r(134) & + r(136)+( 1.195)*r(137)+( 0.150)*r(138) & + r(139)+ r(140)+( 0.300)*r(141) & + r(142)+ r(143)+( 0.082)*r(144) & +( 0.872)*r(146)+( 1.072)*r(147)+( 0.200)*r(149) & + r(150)+( 0.658)*r(151)+( 0.284)*r(153) Gain(lRO2 ) = Gain(lRO2 ) & +( 0.840)*r(154)+( 0.067)*r(155)+ r(158) & +( 0.800)*r(161)+ r(162)+ r(163) & +( 1.500)*r(164)+( 0.940)*r(165)+( 1.280)*r(166) & +( 0.352)*r(167)+ r(169)+ r(171) & +( 0.720)*r(172)+ r(174)+ r(176) & +( 0.602)*r(177)+ r(180)+ r(181) & +( 0.020)*r(182)+( 0.700)*r(183)+( 2.000)*r(190) & +( 0.300)*r(191)+ r(192)+( 0.400)*r(194) Gain(lRO2 ) = Gain(lRO2 ) & +( 0.440)*r(202)+( 2.000)*r(203)+( 1.800)*r(204) & + r(208)+ r(225)+ r(226) & + r(227)+ r(228) c Loss(lPACD )= + r( 95) c Gain(lPACD )= +( 0.410)*r( 57)+( 0.410)*r( 65)+( 0.410)*r(202) c Loss(lAACD )= + r( 94) c Gain(lAACD )= +( 0.150)*r( 57)+( 0.150)*r( 65)+( 0.100)*r( 73) & +( 0.200)*r( 77)+( 0.200)*r( 81)+( 0.200)*r( 85) & +( 0.130)*r(138)+( 0.080)*r(141)+( 0.200)*r(146) & +( 0.200)*r(161)+( 0.150)*r(202)+( 0.200)*r(204) c Loss(lCXO3 )= + r( 60)+ r( 61)+ r( 62) & + r( 65)+ r( 66)+( 2.000)*r( 67) c Gain(lCXO3 )= + r( 63)+( 0.600)*r( 64)+ r(107) & + r(108)+( 0.500)*r(125)+( 0.200)*r(149) & +( 0.717)*r(153)+( 0.390)*r(165)+( 0.200)*r(199) c Loss(lALD2 )= + r(104)+ r(105)+ r(106) c Gain(lALD2 )= + r( 60)+ r( 61)+( 0.400)*r( 64) & +( 0.440)*r( 65)+( 0.800)*r( 66)+( 2.000)*r( 67) & + r(109)+( 0.991)*r(122)+( 0.950)*r(124) & +( 0.500)*r(125)+( 0.740)*r(130)+( 0.488)*r(137) & +( 0.295)*r(138)+( 0.250)*r(139)+( 1.300)*r(140) & +( 0.732)*r(141)+( 0.500)*r(142)+( 0.040)*r(152) & +( 0.100)*r(191)+( 0.020)*r(195)+ r(206) c Loss(lXO2H )= + r( 75)+ r( 76)+ r( 77) & + r( 78) c Gain(lXO2H )= + r( 60)+ r( 61)+( 0.400)*r( 64) & +( 0.440)*r( 65)+( 0.800)*r( 66)+( 2.000)*r( 67) & +( 0.540)*r( 89)+ r(109)+( 0.110)*r(111) & +( 0.991)*r(122)+( 0.100)*r(124)+( 0.500)*r(125) & +( 0.940)*r(130)+ r(134)+( 0.500)*r(136) & +( 0.976)*r(137)+( 0.150)*r(138)+( 0.480)*r(139) & + r(140)+( 0.300)*r(141)+( 0.480)*r(142) & +( 0.082)*r(144)+( 0.072)*r(146)+( 0.072)*r(147) & +( 0.640)*r(150)+( 0.142)*r(153)+( 0.340)*r(154) Gain(lXO2H ) = Gain(lXO2H ) & +( 0.370)*r(163)+( 0.750)*r(164)+( 0.070)*r(165) & +( 0.280)*r(166)+( 0.070)*r(172)+( 0.058)*r(177) & +( 0.120)*r(183)+ r(189)+( 2.000)*r(190) & +( 0.300)*r(191)+( 0.450)*r(192)+( 0.400)*r(194) & +( 0.440)*r(202)+( 0.800)*r(204)+ r(226) & +( 0.126)*r(228) c Loss(lPANX )= + r( 63)+ r( 64)+ r(206) c Gain(lPANX )= + r( 62) c Loss(lFORM )= + r( 96)+ r( 97)+ r( 98) & + r( 99)+ r(100) c Gain(lFORM )= + r( 71)+( 0.100)*r( 72)+ r( 73) & +( 0.685)*r( 74)+( 0.400)*r( 87)+ r(101) & +( 0.740)*r(111)+ r(123)+( 0.078)*r(124) & + r(127)+( 1.560)*r(134)+ r(135) & +( 1.125)*r(136)+( 0.781)*r(137)+( 0.555)*r(138) & +( 0.500)*r(139)+( 0.128)*r(141)+( 0.673)*r(144) & +( 0.120)*r(145)+( 0.598)*r(146)+( 0.598)*r(147) & +( 0.600)*r(149)+( 0.350)*r(150)+( 0.231)*r(152) & +( 0.260)*r(154)+( 0.375)*r(159)+( 0.375)*r(160) Gain(lFORM ) = Gain(lFORM ) & +( 0.300)*r(161)+( 0.375)*r(162)+( 0.592)*r(163) & +( 0.280)*r(164)+( 0.240)*r(165)+ r(188) & +( 0.080)*r(195)+( 0.140)*r(197)+ r(231) & + r(233) c Loss(lMEPX )= + r( 87)+ r( 88) c Gain(lMEPX )= +( 0.900)*r( 72)+( 0.500)*r(103) c Loss(lMEOH )= + r(123) c Gain(lMEOH )= +( 0.315)*r( 74)+( 0.150)*r(111) c Loss(lROOH )= + r( 89)+ r( 90) c Gain(lROOH )= + r( 76)+ r( 80)+ r( 84) c Loss(lXO2 )= + r( 79)+ r( 80)+ r( 81) & + r( 82) c Gain(lXO2 )= +( 0.200)*r(113)+( 0.500)*r(115)+ r(117) & + r(127)+( 0.500)*r(136)+( 0.195)*r(137) & +( 0.480)*r(139)+( 0.480)*r(142)+( 0.200)*r(149) & +( 0.330)*r(150)+( 0.521)*r(151)+( 0.142)*r(153) & +( 0.160)*r(154)+( 0.630)*r(163)+( 0.500)*r(164) & +( 0.690)*r(165)+( 0.750)*r(166)+( 0.480)*r(183) & +( 0.450)*r(192)+ r(203)+( 0.874)*r(228) c Loss(lXO2N )= + r( 83)+ r( 84)+ r( 85) & + r( 86) c Gain(lXO2N )= +( 0.060)*r( 89)+( 0.009)*r(122)+( 0.040)*r(130) & +( 0.024)*r(137)+( 0.040)*r(139)+( 0.040)*r(142) & +( 0.030)*r(150)+( 0.022)*r(151)+( 0.250)*r(164) & +( 0.180)*r(165)+( 0.250)*r(166)+( 0.020)*r(182) & +( 0.100)*r(183)+( 0.100)*r(192)+ r(225) & + r(227) c Loss(lNTR1 )= + r( 91)+ r( 92) c Gain(lNTR1 )= +( 0.500)*r( 83)+ r(132)+( 0.500)*r(136) & +( 0.500)*r(139)+( 0.500)*r(142) c Loss(lNTR2 )= + r(207) c Gain(lNTR2 )= +( 0.500)*r( 83)+ r( 91)+( 0.650)*r(150) & +( 0.142)*r(153)+( 0.266)*r(163)+( 0.530)*r(166) & +( 0.082)*r(168)+( 0.140)*r(173)+( 0.140)*r(178) & + r(186)+ r(187)+( 0.500)*r(192) & +( 0.500)*r(205) c Loss(lFACD )= + r( 93) c Gain(lFACD )= + r(102)+( 0.500)*r(103)+( 0.300)*r(133) & +( 0.370)*r(135)+( 0.090)*r(138)+( 0.150)*r(152) & +( 0.074)*r(159)+( 0.185)*r(163) c Loss(lCO )= + r(120) c Gain(lCO )= + r( 96)+ r( 97)+ r( 98) & + r( 99)+ r(106)+ r(109) & +( 0.890)*r(111)+( 1.800)*r(113)+( 2.000)*r(114) & +( 1.500)*r(115)+ r(116)+ r(118) & +( 0.380)*r(126)+( 0.300)*r(133)+( 0.510)*r(135) & +( 0.378)*r(138)+( 0.245)*r(141)+( 0.066)*r(149) & +( 0.137)*r(151)+( 0.543)*r(152)+( 0.251)*r(159) & +( 0.251)*r(160)+( 0.200)*r(161)+( 0.251)*r(162) & +( 0.001)*r(165)+( 0.700)*r(189)+( 0.500)*r(191) Gain(lCO ) = Gain(lCO ) & + r(193)+( 1.980)*r(195)+( 0.500)*r(199) & + r(205) c Loss(lHCO3 )= + r(101)+ r(102)+ r(103) c Gain(lHCO3 )= + r(100) c Loss(lALDX )= + r(107)+ r(108)+ r(109) c Gain(lALDX )= +( 0.370)*r(130)+( 0.488)*r(137)+( 0.270)*r(138) & +( 0.375)*r(139)+( 0.700)*r(140)+( 0.442)*r(141) & +( 0.625)*r(142)+( 0.150)*r(149)+( 0.029)*r(155) & +( 0.078)*r(163)+( 0.470)*r(164)+( 0.210)*r(165) & +( 0.470)*r(166)+( 0.440)*r(202)+ r(203) & +( 0.800)*r(204)+( 0.268)*r(226)+( 0.126)*r(228) c Loss(lGLYD )= + r(110)+ r(111)+ r(112) c Gain(lGLYD )= +( 0.011)*r(124)+( 0.220)*r(134)+( 0.269)*r(151) & +( 0.113)*r(153)+( 0.128)*r(154)+( 0.275)*r(159) & +( 0.275)*r(160)+( 0.220)*r(161)+( 0.275)*r(162) & +( 0.331)*r(163) c Loss(lGLY )= + r(113)+ r(114)+ r(115) c Gain(lGLY )= +( 0.200)*r(110)+( 0.110)*r(111)+( 0.700)*r(133) & +( 0.075)*r(138)+( 0.240)*r(141)+( 0.170)*r(152) & +( 0.275)*r(159)+( 0.275)*r(160)+( 0.220)*r(161) & +( 0.275)*r(162)+( 0.918)*r(168)+ r(169) & + r(171)+( 0.417)*r(173)+( 0.480)*r(174) & +( 0.480)*r(176)+( 0.221)*r(178)+( 0.260)*r(180) & +( 0.260)*r(181)+( 0.025)*r(182)+( 0.240)*r(183) & +( 0.400)*r(189)+( 0.400)*r(190)+( 0.400)*r(194) & +( 1.400)*r(195)+( 0.500)*r(199)+( 0.500)*r(205) c Loss(lMGLY )= + r(116)+ r(117)+ r(118) c Gain(lMGLY )= +( 0.075)*r(138)+( 0.060)*r(141)+( 0.115)*r(151) & +( 0.531)*r(152)+( 0.113)*r(153)+( 0.275)*r(159) & +( 0.275)*r(160)+( 0.220)*r(161)+( 0.275)*r(162) & +( 0.443)*r(173)+( 0.520)*r(174)+( 0.520)*r(176) & +( 0.675)*r(178)+( 0.770)*r(180)+( 0.770)*r(181) & +( 0.240)*r(183)+ r(190)+( 1.200)*r(191) & +( 0.250)*r(192)+( 0.240)*r(195) c Loss(lETHA )= + r(122) c Gain(lETHA )= 0.0 c Loss(lETOH )= + r(124) c Gain(lETOH )= 0.0 c Loss(lKET )= + r(125) c Gain(lKET )= +( 0.200)*r(130)+ r(131) c Loss(lPAR )= + r(129) c Gain(lPAR )= +(-2.500)*r(125)+(-2.700)*r(130)+(-0.730)*r(137) & +(-0.790)*r(138)+(-1.000)*r(139)+( 0.290)*r(141) & + r(142)+( 0.350)*r(149)+( 0.117)*r(151) & +( 0.717)*r(153)+( 0.240)*r(154)+( 2.175)*r(159) & +( 2.175)*r(160)+( 1.740)*r(161)+( 2.175)*r(162) & +( 2.700)*r(163)+( 1.660)*r(164)+( 7.000)*r(165) & +( 0.268)*r(226)+(-0.126)*r(228) c Loss(lACET )= + r(126)+ r(127) c Gain(lACET )= +( 0.420)*r(130)+( 0.137)*r(151)+( 0.170)*r(152) & +( 0.170)*r(154)+( 0.732)*r(226) c Loss(lPRPA )= + r(128) c Gain(lPRPA )= 0.0 c Loss(lXPRP )= + r(225)+ r(226) c Gain(lXPRP )= + r(128) c Loss(lXPAR )= + r(227)+ r(228) c Gain(lXPAR )= + r(129) c Loss(lROR )= + r(130)+ r(131)+ r(132) c Gain(lROR )= +( 0.020)*r(130)+( 0.874)*r(228) c Loss(lETHY )= + r(133) c Gain(lETHY )= 0.0 c Loss(lETH )= + r(134)+ r(135)+ r(136) c Gain(lETH )= 0.0 c Loss(lOLE )= + r(137)+ r(138)+ r(139) c Gain(lOLE )= +( 0.240)*r(154)+( 0.098)*r(163) c Loss(lIOLE )= + r(140)+ r(141)+ r(142) c Gain(lIOLE )= +( 0.029)*r(155) c Loss(lISOP )= + r(143)+ r(149)+ r(150) c Gain(lISOP )= 0.0 c Loss(lISO2 )= + r(144)+ r(145)+ r(146) & + r(147)+ r(148) c Gain(lISO2 )= + r(143)+( 0.067)*r(155) c Loss(lINTR )= + r(163)+ r(229) c Gain(lINTR )= +( 0.100)*r(144)+( 0.104)*r(163) c Loss(lISPD )= + r(151)+ r(152)+ r(153) & + r(154) c Gain(lISPD )= +( 0.900)*r(144)+( 0.120)*r(145)+ r(146) & + r(147)+( 0.650)*r(149)+( 0.350)*r(150) & + r(156)+ r(157) c Loss(lISPX )= + r(155) c Gain(lISPX )= +( 0.880)*r(145) c Loss(lHPLD )= + r(156)+ r(157) c Gain(lHPLD )= + r(148) c Loss(lOPO3 )= + r(199)+ r(200)+ r(202) & + r(203)+ r(204) c Gain(lOPO3 )= +( 0.457)*r(151)+( 0.480)*r(183)+ r(193) & +( 0.600)*r(194)+ r(196)+ r(201) c Loss(lEPOX )= + r(158) c Gain(lEPOX )= +( 0.904)*r(155) c Loss(lEPX2 )= + r(159)+ r(160)+ r(161) & + r(162) c Gain(lEPX2 )= + r(158) c Loss(lTERP )= + r(164)+ r(165)+ r(166) c Gain(lTERP )= 0.0 c Loss(lBENZ )= + r(167) c Gain(lBENZ )= 0.0 c Loss(lCRES )= + r(182)+ r(183) c Gain(lCRES )= +( 0.530)*r(167)+( 0.180)*r(172)+( 0.155)*r(177) & + r(185) c Loss(lBZO2 )= + r(168)+ r(169)+ r(170) & + r(171) c Gain(lBZO2 )= +( 0.352)*r(167) c Loss(lOPEN )= + r(193)+ r(194)+ r(195) & + r(196) c Gain(lOPEN )= +( 0.118)*r(167)+( 0.918)*r(168)+ r(169) & + r(171)+( 0.100)*r(172)+( 0.660)*r(173) & +( 0.770)*r(174)+( 0.770)*r(176)+( 0.300)*r(178) & +( 0.350)*r(180)+( 0.350)*r(181)+( 0.025)*r(182) & + r(188)+( 0.250)*r(192) c Loss(lTOL )= + r(172) c Gain(lTOL )= 0.0 c Loss(lTO2 )= + r(173)+ r(174)+ r(175) & + r(176) c Gain(lTO2 )= +( 0.650)*r(172) c Loss(lXOPN )= + r(189)+ r(190)+ r(191) & + r(192) c Gain(lXOPN )= +( 0.200)*r(173)+( 0.230)*r(174)+( 0.230)*r(176) & +( 0.244)*r(177)+( 0.560)*r(178)+( 0.650)*r(180) & +( 0.650)*r(181) c Loss(lXYL )= + r(177) c Gain(lXYL )= 0.0 c Loss(lXLO2 )= + r(178)+ r(179)+ r(180) & + r(181) c Gain(lXLO2 )= +( 0.544)*r(177) c Loss(lCRO )= + r(184)+ r(185) c Gain(lCRO )= +( 0.200)*r(182)+( 0.300)*r(183)+( 0.500)*r(186) & +( 0.500)*r(187)+( 0.500)*r(197)+ r(198) c Loss(lCAT1 )= + r(197)+ r(198) c Gain(lCAT1 )= +( 0.732)*r(182) c Loss(lCRON )= + r(186)+ r(187)+ r(188) c Gain(lCRON )= + r(184) c Loss(lOPAN )= + r(201)+ r(205) c Gain(lOPAN )= + r(200) c Loss(lECH4 )= + r(208) c Gain(lECH4 )= 0.0 c Loss(lI2 )= + r(209) c Gain(lI2 )= 0.0 c Loss(lHOI )= + r(210) c Gain(lHOI )= + r(214)+ r(224) c Loss(lI2O2 )= + r(221)+ r(222) c Gain(lI2O2 )= +( 0.600)*r(213) c Loss(lINO3 )= + r(223)+ r(224) c Gain(lINO3 )= + r(216) c Loss(lHIO3 )= 0.0 c Gain(lHIO3 )= + r(218) c Loss(lIXOY )= 0.0 c Gain(lIXOY )= + r(219)+ r(222) c Loss(lDMS )= + r(231)+ r(232)+ r(233) c Gain(lDMS )= 0.0 c c return end