* ************************************************************** ; * gas phase: RACM ; * aqueous phase: CAPRAM 2.4; * urban conditions; * ************************************************************** ; * Parameter definitions; * R: Perfect gas constant; * pi: pi; * TEMP: temperature; * NAVOG: Avogadro's number; * ZCL: Height of the boundary layer; * DG: Gas phase diffusion coefficient; * L: water content; * For the chemical species the following nomenclature is applied: XXXg: gas phase XXXa: aqueous phase (neutral or radical) XXXm: anion XXXp: cation * KG Gas phase reaction; * KA Aqueous phase reaction; * J Photolysis rate constant; * Chosen units; * Length: cm; * Volume: cm3; * Pressure: atm; * ; * *********************************************** ; * ; * Gas phase species used in RACM; * ; * *********************************************** ; VARIABLE O3g H2O2g NOg NO2g NO3g N2O5g HNO2g HNO3g HNO4g SO2g SULFg COg CO2g N2g O2g H2Og H2g O3Pg O1Dg OHg HO2g CH4g ETHg HC3g HC5g HC8g ETEg OLTg OLIg DIENg ISOg APIg LIMg TOLg XYLg CSLg HCHOg ALDg KETg GLYg MGLYg DCBg MACRg UDDg HKETg ONITg PANg TPANg OP1g OP2g PAAg ORA1g ORA2g MO2g HC3P HC5P HC8P ETEPg OLTPg OLIPg ISOPg APIPg LIMPg PHOg ADDTg ADDXg ADDCg TOLPg XYLPg CSLPg ACO3g TCO3g KETPg OLNNg OLNDg XO2g HC3Pg HC5Pg HC8Pg KHETg XOg CH3OHg OL2Pg ETOHg ETHPg DMSg NH3g N2O4g NOClg HClg ClNO3g HBrg HOBrg BrNO3g Br2g BrClg BrNO2g ClNO2g HOClg Cl2g SO3g BrOg ClOg ; * ; * ************************************************ ; * ; * Chemical species used in CAPRAM ; * ; * ************************************************ ; * ; VARIABLE FeOHpp FeOH2p CH3COOm ACO3a Br2a Br2m Bra BrCla Brm BrNO2a BrOHm COH22a HOA1a C2O4m C2O4mm CHOH2a CHOH22a HOA2a HOA2m CH2OH2a CH2COOHa CH2OHa HOA3m CH3a CH3COOOm CH3COH2a HOP1a ALCOOL1a CH3CH2Oa CH3CHOa CH3CHOHa CH3COa CH3COOHa CH3O2a CH3Oa CH3OHa CH3O2Ha CHOHSO3m CHOSO3m Cl2m Clm ClNO2a ClOHm CO2a CO2Ha CO3m CO3mm COOHa Cup Cupp ETHPa ETOHa FeMO2pp FeClpp FeOpp Fepp Feppp H2C2O4a H2O2a H2Oa HC2O4m HCHOa HCO3m HCOOHa HCOOm HMSm HNO2a HO2a HOBra HOCla HNO3a Hp HSO3m HSO5m MnOHpp MnBr2p MnCl2p MnHSO3p MnO2p MnOpp Mnpp Mnppp Mnpppp N2O5a NO2a NO2m NO2p NO3a NO3m O2a HOP2a HOA3a HOA4m O2CH2OHa HOP3a O2CHOa HS2m O2m O2NO2m O3a O3m OHa OHm S2O8mm SO2a SO3a SO3m SO3mm SO4m SO4mm SO5m SO5O2mm OHCH2SO3a CH2COOm HCla NH3a Cl2a HSO4m PAAa GLYa HNO4a HO3a HFe4p OHCH2SO3a H2SO4a NH4p H2CO3a SO5O2Hm FeSO4p Cla MnOH2p CuOHpp CHOSO3mm FeC2O4p NOa FeC2O42m FeORG3 HOONOa OONOm O3Pa CO2m hours; * Renamed species: ; * COH2CHOH2a - COH22a COH2COOHa - HOA1a CHOH2CHOH2a - CHOH22a CHOH2COOHa - HOA2a CHOH2COOm - HOA2m CH2OHCOOm - HOA3m CH3COH2O2a - HOP1a CH3CHOH2a - ALCOOL1a FeCH3O2pp - FeMO2pp O2COH2CHOH2a - HOP2a O2COH2COOHa - HOA3a O2CH2COOm - HOA4m O2CH3CHOHa - HOP3a O2CHOHSO3m - HS2m FeOH2Fepppp - Hfe4p FeC2O43mmm - FeORG3 ; * ; * *************************************************** ; * ; PARAMETER KCAP1 KCAP2 KCAP3 KCAP4 KCAP5 KCAP6 KCAP7 KCAP8 KCAP9 KCAP10 KCAP11 KCAP12 KCAP13 KCAP14 KCAP15 KCAP16 KCAP17 KCAP18 KCAP19 KCAP20 KCAP21 KCAP22 KCAP23 KCAP24 KCAP25 KCAP26 KCAP27 KCAP28 KCAP29 KCAP30 KCAP31 KCAP32 KCAP33 KCAP34 KCAP35 KCAP36 KCAP37 KCAP38 KCAP39 KCAP40 KCAP41 KCAP42 KCAP43 KCAP44 KCAP45 KCAP46 KCAP47 KCAP48 KCAP49 KCAP50 KCAP51 KCAP52 KCAP53 KCAP54 KCAP55 KCAP56 KCAP57 KCAP58 KCAP59 KCAP60 KCAP61 KCAP62 KCAP63 KCAP64 KCAP65 KCAP66 KCAP67 KCAP68 KCAP69 KCAP70 KCAP71 KCAP72 KCAP73 KCAP74 KCAP75 KCAP76 KCAP77 KCAP78 KCAP79 KCAP80 KCAP81 KCAP82 KCAP83 KCAP84 KCAP85 KCAP86 KCAP87 KCAP88 KCAP89 KCAP90 KCAP91 KCAP92 KCAP93 KCAP94 KCAP95 KCAP96 KCAP97 KCAP98 KCAP99 KCAP100 KCAP101 KCAP102 KCAP103 KCAP104 KCAP105 KCAP106 KCAP107 KCAP108 KCAP109 KCAP110 KCAP111 KCAP112 KCAP113 KCAP114 KCAP115 KCAP116 KCAP117 KCAP118 KCAP119 KCAP120 KCAP121 KCAP122 KCAP123 KCAP124 KCAP125 KCAP126 KCAP127 KCAP128 KCAP129 KCAP130 KCAP131 KCAP132 KCAP133 KCAP134 KCAP135 KCAP136 KCAP137 KCAP138 KCAP139 KCAP140 KCAP141 KCAP142 KCAP143 KCAP144 KCAP145 KCAP146 KCAP147 KCAP148 KCAP149 KCAP150 KCAP151 KCAP152 KCAP153 KCAP154 KCAP155 KCAP156 KCAP157 KCAP158 KCAP159 KCAP160 KCAP161 KCAP162 KCAP163 KCAP164 KCAP165 KCAP166 KCAP167 KCAP168 KCAP169 KCAP170 KCAP171 KCAP172 KCAP173 KCAP174 KCAP175 KCAP176 KCAP177 KCAP178 KCAP179 KCAP180 KCAP181 KCAP182 KCAP183 KCAP184 KCAP185 KCAP186 KCAP187 KCAP188 KCAP189 KCAP190 KCAP191 KCAP192 KCAP193 KCAP194 KCAP195 KCAP196 KCAP197 KCAP198 KCAP199 KCAP200 KCAP201 KCAP202 KCAP203 KCAP204 KCAP205 KCAP206 KCAP207 KCAP208 KCAP209 KCAP210 KCAP211 KCAP212 KCAP213 KCAP214 KCAP215 KCAP216 KCAP217 KCAP218 KCAP219 KCAP220 KCAP221 KCAP222 KCAP223 KCAP224 KCAP225 KCAP226 KCAP227 KCAP228 KCAP229 KCAP230 KCAP231 KCAP232 KCAP233 KCAP234 KCAP235 KCAP236 KCAP237 KCAP238 KCAP239 KCAP240 KCAP241 KCAP242 KCAP243 KCAP244 KCAP245 KCAP246 KCAP247 KCAP248 KCAP249 KCAP250 KCAP251 KCAP252 KCAP253 KCAP254 KCAP255 KCAP256 KCAP257 KCAP258 KCAP259 KCAP260 KCAP261 KCAP262 KCAP263 KCAP264 KCAP265 KCAP266 KCAP267 KCAP268 KCAP269 KCAP270 KCAP271 KCAP272 KCAP273 KCAP274 KCAP275 KCAP276 KCAP277 KCAP278 KCAP279 KCAP280 KCAP281 KCAP282 KCAP283 KCAP284 KCAP285 KCAP286 KCAP287 KCAP288 KCAP289 KCAP290 KCAP291 KCAP292 KCAP293 KCAP294 KCAP295 KCAP296 KCAP297 KCAP298 KCAP299 KCAP300 KCAP301 KCAP302 KCAP303 KCAP304 KCAP305 KCAP306 KCAP307 KCAP308 KCAP309 KCAP310 KCAP311 KCAP312 KCAP313 E1FWD E2FWD E3FWD E4FWD E5FWD E6FWD E7FWD E8FWD E9FWD E10FWD E11FWD E12FWD E13FWD E14FWD E15FWD E16FWD E17FWD E18FWD E19FWD E20FWD E21FWD E22FWD E23FWD E24FWD E25FWD E26FWD E27FWD E28FWD E29FWD E30FWD E31FWD E32FWD E33FWD E34FWD E35FWD E36FWD E37FWD E38FWD E39FWD E40FWD E41FWD E42FWD E43FWD E44FWD E45FWD E46FWD E47FWD E48FWD E49FWD E50FWD E51FWD E52FWD E53FWD E54FWD E55FWD E56FWD E57FWD E58FWD E1BCK E2BCK E3BCK E4BCK E5BCK E6BCK E7BCK E8BCK E9BCK E10BCK E11BCK E12BCK E13BCK E14BCK E15BCK E16BCK E17BCK E18BCK E19BCK E20BCK E21BCK E22BCK E23BCK E24BCK E25BCK E26BCK E27BCK E28BCK E29BCK E30BCK E31BCK E32BCK E33BCK E34BCK E35BCK E36BCK E37BCK E38BCK E39BCK E40BCK E41BCK E42BCK E43BCK E44BCK E45BCK E46BCK E47BCK E48BCK E49BCK E50BCK E51BCK E52BCK E53BCK E54BCK E55BCK E56BCK E57BCK E58BCK P1 P2 P3 P4 P5 P6 P7 P8 P9 P10 P11 J1 J2 J3 J4 J5 J6 J7 J8 J9 J10 J11 J12 J13 J14 J15 J16 J17 J18 J19 J20 J21 J22 J23 K24 K25 K26 K27 K28 K29 K30 K31 K32 K33 K34 K35 K36 K37 K38 K39 K40 K41 K42 K43 K44 K45 K46 K47 K48 K49 K50 K51 K52 K53 K54 K55 K56 K57 K58 K59 K60 K61 K62 K63 K64 K65 K66 K67 K68 K69 K70 K71 K72 K73 K74 K75 K76 K77 K78 K79 K80 K81 K82 K83 K84 K85 K86 K87 K88 K89 K90 K91 K92 K93 K94 K95 K96 K97 K98 K99 K100 K101 K102 K103 K104 K105 K106 K107 K108 K109 K110 K111 K112 K113 K114 K115 K116 K117 K118 K119 K120 K121 K122 K123 K124 K125 K126 K127 K128 K129 K130 K131 K132 K133 K134 K135 K136 K137 K138 K139 K140 K141 K142 K143 K144 K145 K146 K147 K148 K149 K150 K151 K152 K153 K154 K155 K156 K157 K158 K159 K160 K161 K162 K163 K164 K165 K166 K167 K168 K169 K170 K171 K172 K173 K174 K175 K176 K177 K178 K179 K180 K181 K182 K183 K184 K185 K186 K187 K188 K189 K190 K191 K192 K193 K194 K195 K196 K197 K198 K199 K200 K201 K202 K203 K204 K205 K206 K207 K208 K209 K210 K211 K212 K213 K214 K215 K216 K217 K218 K219 K220 K221 K222 K223 K224 K225 K226 K227 K228 K229 K230 K231 K232 K233 K234 K235 K236 K237 ; * *************************************************************** ; * ************************************************************ ; PARAMETER K11ZK10 FACT 1000; PARAMETER Relhum Psat Rgas Time0 OneDay InitSpr Latitu Cancer SolarDecli SINPSI DayReal Delai sza DAY0 start z hourangle Longcor lat eqt teta; * ; * ************************************************************** ; * Depositions, emissions, phase exchange; * ****************************************************************; * ; PARAMETER ACO3 ACO2 ACOH ACHO2 ACH2O2 ACNO3 ACN2O5 ACHNO3 ACNH3 ACMO2 ACOP1 ACHCHO ACORA1 ACCO2 ACHCL ACHOCL ACCLNO3 ACCL2 ACHBR ACHOBR ACBRNO3 ACBR2 ACBRCL ACSO2 ACSO3 ACHNO2 ACCLNO2 ACBRNO2 ACNOCL ACN2O4 ACNO2 ACHI ACHOI ACINO3 ACI2 ACICL ACINO2 ACACO3 ACALD ACBRO ACCH3OH ACCLO ACETH ACETHP ACETOH ACOL2 ACOL2P ACOP2 ACORA2 ACPAA ACPAN ACSULF ACGLY ACHNO4 ACNO HCO3 HCO2 HCOH HCHO2 HCH2O2 HCNO3 HCN2O5 HCHNO3 HCNH3 HCMO2 HCOP1 HCHCHO HCORA1 HCCO2 HCHCL HCHOCL HCCLNO3 HCCL2 HCHBR HCHOBR HCBRNO3 HCBR2 HCBRCL HCSO2 HCSO3 HCHNO2 HCCLNO2 HCBRNO2 HCN2O4 HCNO2 HCHI HCHOI HCINO3 HCI2 HCICL HCINO2 HCACO3 HCALD HCBRO HCCH3OH HCCLO HCETH HCETHP HCETOH HCOL2 HCOL2P HCOP2 HCORA2 HCPAA HCPAN HCSULF HCHNO4 HCCH4 HCGLY HCNO MMO3 MMO2 MMOH MMHO2 MMH2O2 MMNO3 MMN2O5 MMHNO3 MMNH3 MMMO2 MMOP1 MMHCHO MMORA1 MMCO2 MMHCL MMHOCL MMCLNO3 MMCL2 MMHBR MMHOBR MMBRNO3 MMBR2 MMBRCL MMSO2 MMSO3 MMHNO2 MMCLNO2 MMBRNO2 MMNOCL MMN2O4 MMNO2 MMHI MMHOI MMINO3 MMI2 MMINO2 MMICL MMACO3 MMALD MMBRO MMCH3OH MMCLO MMETH MMETHP MMETOH MMOL2 MMOL2P MMOP2 MMORA2 MMPAA MMPAN MMSULF MMGLY MMHNO4 MMNO DGCO2 DGHCl DGNH3 DGO3 DGHO2 DGOH DGH2O2 DGHNO3 DGNO3 DGN2O5 DGNO2 DGHNO2 DGHNO4 DGSO2 DGHCHO DGOP1 DGPAA DGCH3OH DGETOH DGALD DGORA1 DGORA2 DGMO2 DGETHP DGCl2 DGBr2 DGSULF DGACO3 DGGLY DGO2 DGClNO2 DGBrNO2 DGBrCl DGNO VO3 VO2 VOH VHO2 VH2O2 VNO3 VN2O5 VHNO3 VNH3 VMO2 VOP1 VHCHO VORA1 VCO2 VHCL VHOCL VCLNO3 VCL2 VHBR VHOBR VBRNO3 VBR2 VBRCL VSO2 VSO3 VHNO2 VCLNO2 VBRNO2 VNOCL VNO2 VN2O4 VHI VHOI VINO3 VI2 VICL VINO2 VACO3 VALD VBRO VCH3OH VCLO VETH VETHP VETOH VOL2 VOL2P VOP2 VORA2 VPAA VPAN VSULF VGLY VHNO4 VNO VDCO VDHOBr VDHOCl VDNO2 VDO3 VDETOH VDOP1 VDCH3OH VDH2O2 VDHBr VDHCHO VDHCl VDNH3 VDORA1 VDSO2 VDHNO3 VDN2O5 VDSULF KMTO3 KMTO2 KMTOH KMTHO2 KMTH2O2 KMTNO3 KMTN2O5 KMTHNO3 KMTNH3 KMTMO2 KMTOP1 KMTHCHO KMTORA1 KMTCO2 KMTHCL KMTHOCL KMTCLNO3 KMTCL2 KMTHBR KMTHOBR KMTBRNO3 KMTBR2 KMTBRCL KMTSO2 KMTSO3 KMTHNO2 KMTCLNO2 KMTBRNO2 KMTNOCL KMTN2O4 KMTNO2 KMTHI KMTHOI KMTINO3 KMTI2 KMTICL KMTINO2 KMTACO3 KMTALD KMTBRO KMTCH3OH KMTCLO KMTETH KMTETHP KMTETOH KMTOL2 KMTOL2P KMTOP2 KMTORA2 KMTPAA KMTPAN KMTSULF KMTGLY KMTHNO4 KMTNO EMICH3OH EMIETOH EMIALC EMIHC3 EMIXYL EMIKET EMIHC5 EMIETH EMIOLT EMIALD EMIHCHO EMITOL EMIDIEN EMICSL EMIETE EMIHC8 EMICO EMITERPENE EMIORA2 EMISO2 EMINH3 EMINO EMIAPI EMILIM EMIDMS EMIISO ; PARAMETER <100> WORK; * ; * *****************************************************************; * ; * Definition of general parameters; * ; * *****************************************************************; * ; PARAMETER TEMP TCORR R PI NAVOG RAD VOLUME ZCL CM DG L LOCEAN BETA KSS ON 0 HET 0 k0t kinft DUSK PHOTON CONV1 CONV2 CONV3 CONVJ TMP0 TMP1 TMP2 TMP3 P PNTCOLL KEQRET 1E10 PRES 1 Mg SZA; * ***************************************************; * ; VARIABLE KCAP1t KCAP2t KCAP3t KCAP4t KCAP5t KCAP6t KCAP7t KCAP8t KCAP9t KCAP10t KCAP11t KCAP12t KCAP13t KCAP14t KCAP15t KCAP16t KCAP17t KCAP18t KCAP19t KCAP20t KCAP21t KCAP22t KCAP23t KCAP24t KCAP25t KCAP26t KCAP27t KCAP28t KCAP29t KCAP30t KCAP31t KCAP32t KCAP33t KCAP34t KCAP35t KCAP36t KCAP37t KCAP38t KCAP39t KCAP40t KCAP41t KCAP42t KCAP43t KCAP44t KCAP45t KCAP46t KCAP47t KCAP48t KCAP49t KCAP50t KCAP51t KCAP52t KCAP53t KCAP54t KCAP55t KCAP56t KCAP57t KCAP58t KCAP59t KCAP60t KCAP61t KCAP62t KCAP63t KCAP64t KCAP65t KCAP66t KCAP67t KCAP68t KCAP69t KCAP70t KCAP71t KCAP72t KCAP73t KCAP74t KCAP75t KCAP76t KCAP77t KCAP78t KCAP79t KCAP80t KCAP81t KCAP82t KCAP83t KCAP84t KCAP85t KCAP86t KCAP87t KCAP88t KCAP89t KCAP90t KCAP91t KCAP92t KCAP93t KCAP94t KCAP95t KCAP96t KCAP97t KCAP98t KCAP99t KCAP100t KCAP101t KCAP102t KCAP103t KCAP104t KCAP105t KCAP106t KCAP107t KCAP108t KCAP109t KCAP110t KCAP111t KCAP112t KCAP113t KCAP114t KCAP115t KCAP116t KCAP117t KCAP118t KCAP119t KCAP120t KCAP121t KCAP122t KCAP123t KCAP124t KCAP125t KCAP126t KCAP127t KCAP128t KCAP129t KCAP130t KCAP131t KCAP132t KCAP133t KCAP134t KCAP135t KCAP136t KCAP137t KCAP138t KCAP139t KCAP140t KCAP141t KCAP142t KCAP143t KCAP144t KCAP145t KCAP146t KCAP147t KCAP148t KCAP149t KCAP150t KCAP151t KCAP152t KCAP153t KCAP154t KCAP155t KCAP156t KCAP157t KCAP158t KCAP159t KCAP160t KCAP161t KCAP162t KCAP163t KCAP164t KCAP165t KCAP166t KCAP167t KCAP168t KCAP169t KCAP170t KCAP171t KCAP172t KCAP173t KCAP174t KCAP175t KCAP176t KCAP177t KCAP178t KCAP179t KCAP180t KCAP181t KCAP182t KCAP183t KCAP184t KCAP185t KCAP186t KCAP187t KCAP188t KCAP189t KCAP190t KCAP191t KCAP192t KCAP193t KCAP194t KCAP195t KCAP196t KCAP197t KCAP198t KCAP199t KCAP200t KCAP201t KCAP202t KCAP203t KCAP204t KCAP205t KCAP206t KCAP207t KCAP208t KCAP209t KCAP210t KCAP211t KCAP212t KCAP213t KCAP214t KCAP215t KCAP216t KCAP217t KCAP218t KCAP219t KCAP220t KCAP221t KCAP222t KCAP223t KCAP224t KCAP225t KCAP226t KCAP227t KCAP228t KCAP229t KCAP230t KCAP231t KCAP232t KCAP233t KCAP234t KCAP235t KCAP236t KCAP237t KCAP238t KCAP239t KCAP240t KCAP241t KCAP242t KCAP243t KCAP244t KCAP245t KCAP246t KCAP247t KCAP248t KCAP249t KCAP250t KCAP251t KCAP252t KCAP253t KCAP254t KCAP255t KCAP256t KCAP257t KCAP258t KCAP259t KCAP260t KCAP261t KCAP262t KCAP263t KCAP264t KCAP265t KCAP266t KCAP267t KCAP268t KCAP269t KCAP270t KCAP271t KCAP272t KCAP273t KCAP274t KCAP275t KCAP276t KCAP277t KCAP278t KCAP279t KCAP280t KCAP281t KCAP282t KCAP283t KCAP284t KCAP285t KCAP286t KCAP287t KCAP288t KCAP289t KCAP290t KCAP291t KCAP292t KCAP293t KCAP294t KCAP295t KCAP296t KCAP297t KCAP298t KCAP299t KCAP300t KCAP301t KCAP302t KCAP303t KCAP304t KCAP305t KCAP306t KCAP307t KCAP308t KCAP309t KCAP310t KCAP311t KCAP312t KCAP313t P1t P2t P3t P4t P5t P6t P7t P8t P9t P10t P11t ; * ***************************************************; * ; PARAMETER KCAP1r KCAP2r KCAP3r KCAP4r KCAP5r KCAP6r KCAP7r KCAP8r KCAP9r KCAP10r KCAP11r KCAP12r KCAP13r KCAP14r KCAP15r KCAP16r KCAP17r KCAP18r KCAP19r KCAP20r KCAP21r KCAP22r KCAP23r KCAP24r KCAP25r KCAP26r KCAP27r KCAP28r KCAP29r KCAP30r KCAP31r KCAP32r KCAP33r KCAP34r KCAP35r KCAP36r KCAP37r KCAP38r KCAP39r KCAP40r KCAP41r KCAP42r KCAP43r KCAP44r KCAP45r KCAP46r KCAP47r KCAP48r KCAP49r KCAP50r KCAP51r KCAP52r KCAP53r KCAP54r KCAP55r KCAP56r KCAP57r KCAP58r KCAP59r KCAP60r KCAP61r KCAP62r KCAP63r KCAP64r KCAP65r KCAP66r KCAP67r KCAP68r KCAP69r KCAP70r KCAP71r KCAP72r KCAP73r KCAP74r KCAP75r KCAP76r KCAP77r KCAP78r KCAP79r KCAP80r KCAP81r KCAP82r KCAP83r KCAP84r KCAP85r KCAP86r KCAP87r KCAP88r KCAP89r KCAP90r KCAP91r KCAP92r KCAP93r KCAP94r KCAP95r KCAP96r KCAP97r KCAP98r KCAP99r KCAP100r KCAP101r KCAP102r KCAP103r KCAP104r KCAP105r KCAP106r KCAP107r KCAP108r KCAP109r KCAP110r KCAP111r KCAP112r KCAP113r KCAP114r KCAP115r KCAP116r KCAP117r KCAP118r KCAP119r KCAP120r KCAP121r KCAP122r KCAP123r KCAP124r KCAP125r KCAP126r KCAP127r KCAP128r KCAP129r KCAP130r KCAP131r KCAP132r KCAP133r KCAP134r KCAP135r KCAP136r KCAP137r KCAP138r KCAP139r KCAP140r KCAP141r KCAP142r KCAP143r KCAP144r KCAP145r KCAP146r KCAP147r KCAP148r KCAP149r KCAP150r KCAP151r KCAP152r KCAP153r KCAP154r KCAP155r KCAP156r KCAP157r KCAP158r KCAP159r KCAP160r KCAP161r KCAP162r KCAP163r KCAP164r KCAP165r KCAP166r KCAP167r KCAP168r KCAP169r KCAP170r KCAP171r KCAP172r KCAP173r KCAP174r KCAP175r KCAP176r KCAP177r KCAP178r KCAP179r KCAP180r KCAP181r KCAP182r KCAP183r KCAP184r KCAP185r KCAP186r KCAP187r KCAP188r KCAP189r KCAP190r KCAP191r KCAP192r KCAP193r KCAP194r KCAP195r KCAP196r KCAP197r KCAP198r KCAP199r KCAP200r KCAP201r KCAP202r KCAP203r KCAP204r KCAP205r KCAP206r KCAP207r KCAP208r KCAP209r KCAP210r KCAP211r KCAP212r KCAP213r KCAP214r KCAP215r KCAP216r KCAP217r KCAP218r KCAP219r KCAP220r KCAP221r KCAP222r KCAP223r KCAP224r KCAP225r KCAP226r KCAP227r KCAP228r KCAP229r KCAP230r KCAP231r KCAP232r KCAP233r KCAP234r KCAP235r KCAP236r KCAP237r KCAP238r KCAP239r KCAP240r KCAP241r KCAP242r KCAP243r KCAP244r KCAP245r KCAP246r KCAP247r KCAP248r KCAP249r KCAP250r KCAP251r KCAP252r KCAP253r KCAP254r KCAP255r KCAP256r KCAP257r KCAP258r KCAP259r KCAP260r KCAP261r KCAP262r KCAP263r KCAP264r KCAP265r KCAP266r KCAP267r KCAP268r KCAP269r KCAP270r KCAP271r KCAP272r KCAP273r KCAP274r KCAP275r KCAP276r KCAP277r KCAP278r KCAP279r KCAP280r KCAP281r KCAP282r KCAP283r KCAP284r KCAP285r KCAP286r KCAP287r KCAP288r KCAP289r KCAP290r KCAP291r KCAP292r KCAP293r KCAP294r KCAP295r KCAP296r KCAP297r KCAP298r KCAP299r KCAP300r KCAP301r KCAP302r KCAP303r KCAP304r KCAP305r KCAP306r KCAP307r KCAP308r KCAP309r KCAP310r KCAP311r KCAP312r KCAP313r P1r P2r P3r P4r P5r P6r P7r P8r P9r P10r P11r ; * ***************************************************; * ; * ; * ********************************************************************; * Otput definition; * ********************************************************************; * ; COMPILE INSTANT; OPEN 1 "output1.OUT"; OPEN 2 "output2.OUT"; OPEN 3 "output3.OUT"; OPEN 4 "output4.OUT"; **; COMPILE OUTPUT; WRITE 1=1,"Sim. without depositions for 30 days"; ** ; * ****************************************************; * ; * Parameter initialization ; * ; * ****************************************************; * ; COMPILE INITIAL; * ; * General parameters; * ; PI=3.141592654; R=82.056; Perfect gas constant; TEMP=288; temperature; NAVOG=6.022E23; Avogadro's number; L=3.0E-7; Liquid water content; ZCL=1.0E5; Boundary layer hight cm; RAD=1.0E-3; droplet RAD cm; BETA=0.5; heterogene chemistry branching ratio; KSS=1/(48*3600); renewal speed of aerosols; DG=0.1; gas phase diffusion; TCORR = 1./TEMP-1./298.0; temperature correction term [1/K]; RELHUM = 0.762; relative humidity; PSAT = (6.982616467E+5 - 1.888612677E+4*TEMP + 2.132971127E+2*TEMP**2 - 1.288405237E+0*TEMP**3 + 4.393186046E-3*TEMP**4 - 8.023554873E-6*TEMP**5 + 6.136820929E-9*TEMP**6) /101325.000; * = saturation partial pressure of gaseous H2O [atm]; Mg = PRES / ( R * TEMP ); * = concentration of 'air molecules' [mol/cm3]; CM=Mg*NAVOG; * H2Og = Mg * RELHUM * PSAT / PRES; Latitu=51.*PI/180; Day0=171; 21st June LONGCOR=0. ; * CONVn converts a unit from SI to something else; CONV1 = 1.0 / Mg; [mol/m3(air)] to mixing ratio [1]; CONV2 = NAVOG ; CONV3 = 1000 ; * CONVJ conversion for J values; * PHOTON/COS(LATITU) reaches about 1 on noon of first day of spring; * PI/2. converts from 12h-mean to noon maximum; CONVJ = 0.5*PI/COS(LATITU); * ********************************************************************; * Constants relating to gases; * ********************************************************************; * ; * ********************************************************************; * Henry constants (mol/cm3/atm); * ********************************************************************; *; HCCO2 = 0.001 * 0.031 * exp( 20.14 * 1000 * TCORR /8.314) ; HCHCl = 0.001 * 1.10 * exp( 16.8 * 1000 * TCORR /8.314) ; HCNH3 = 0.001 * 60.7 * exp( 32.6 * 1000 * TCORR /8.314) ; HCO3 = 0.001 * 1.14E-2 * exp( 19.1 * 1000 * TCORR /8.314) ; HCHO2 = 0.001 * 9E3 * exp(0 * 1000 * TCORR /8.314) ; HCOH = 0.001 * 25 * exp( 43.9 * 1000 * TCORR /8.314) ; HCH2O2 = 0.001 * 1.02E5 * exp( 52.7 * 1000 * TCORR /8.314) ; HCHNO3 = 0.001 * 2.1E5 * exp( 72.3 * 1000 * TCORR /8.314) ; HCNO3 = 0.001 * 0.6 * exp( 0 * 1000 * TCORR /8.314) ; HCN2O5 = 0.001 * 1.4 * exp( 0 * 1000 * TCORR /8.314) ; HCNO2 = 0.001 * 1.2E-2 * exp( 10.5 * 1000 * TCORR /8.314) ; HCHNO2 = 0.001 * 49 * exp( 40.6 * 1000 * TCORR /8.314) ; HCHNO4 = 0.001 * 3.0E4 * exp( 0 * 1000 * TCORR /8.314) ; HCSO2 = 0.001 * 1.24 * exp( 27 * 1000 * TCORR /8.314) ; HCHCHO = 0.001 * 2.5 * exp( 60 * 1000 * TCORR /8.314) ; HCOP1 = 0.001 * 310 * exp( 46.6 * 1000 * TCORR /8.314) ; HCPAA = 0.001 * 6.69E2 * exp( 49.0 * 1000 * TCORR /8.314) ; HCCH3OH = 0.001 * 2.2E2 * exp( 44.8 * 1000 * TCORR /8.314) ; HCETOH = 0.001 * 1.9E2 * exp( 52.3 * 1000 * TCORR /8.314) ; HCALD = 0.001 * 4.8 * exp( 52 * 1000 * TCORR /8.314) ; HCORA1 = 0.001 * 5.53E3 * exp( 46.8 * 1000 * TCORR /8.314) ; HCORA2 = 0.001 * 5.50E3 * exp( 49.0 * 1000 * TCORR /8.314) ; HCMO2 = 0.001 * 310 * exp( 46.9 * 1000 * TCORR /8.314) ; HCETHP = 0.001 * 340.0 * exp( 87.0 * TCORR ) ; HCCl2 = 0.001 * 9.15E-2 * exp( 20.7 * 1000 * TCORR /8.314) ; HCBr2 = 0.001 * 0.758 * exp( 31.6 * 1000 * TCORR /8.314) ; HCSULF = 0.001 * 8.7E11 * exp( 0 * 1000 * TCORR /8.314) ; HCACO3 = 0.001 * 669 * exp( 49 * 1000 * TCORR /8.314) ; HCGLY = 0.001 * 1.4 * exp( 0 * 1000 * TCORR /8.314) ; * HCO2 = 0.001 * 1.3E-3 * exp( 14.1 * 1000 * TCORR /8.314) ; HCO2 = 0.001 * 1.3E-3 * exp( 1700 * TCORR ) ; HCClNO2 = 0.001 * 0.024 * exp( 0 * 1000 * TCORR /8.314) ; HCBrNO2 = 0.001 * 0.3 * exp( 0 * 1000 * TCORR /8.314) ; HCBrCl = 0.001 * 0.94 * exp( 0 * 1000 * TCORR /8.314) ; HCNO= 0.001 * 1.9E-3 ; * ; * ********************************************************************; * Accomodation coefficients; * ********************************************************************; * ; ACCO2 = 1.5E-4 ; ACHCl = 0.1158 ; ACNH3 = 0.091 ; ACO3 = 0.1 ; ACHO2 = 0.01 ; ACOH = 0.05 ; ACH2O2 = 0.1532 ; ACHNO3 = 0.0868 ; ACNO3 = 4E-3 ; ACN2O5 = 0.018 ; ACNO2 = 0.0015 ; ACHNO2 = 0.05 ; ACHNO4 = 0.1 ; ACSO2 = 0.11 ; ACHCHO = 0.02 ; ACOP1 = 6.7581E-3 ; ACPAA = 0.019 ; ACCH3OH = 0.0271 ; ACETOH = 0.0176 ; ACALD = 0.03 ; ACORA1 = 0.0229 ; ACORA2 = 0.0322 ; ACMO2 = 6.7581E-3 ; ACETHP = 8.2E-3 ; ACCl2 = 0.08 ; ACBr2 = 0.08 ; ACSULF = 0.12 ; ACACO3 = 0.019 ; ACGLY = 0.023 ; ACO2 = 0.01 ; ACClNO2 = 0.01 ; ACBrNO2 = 0.01 ; ACBrCl = 0.33 ; ACNO = 0.001 ; * ; * *******************************************************************; * Molar mass initializations; * *******************************************************************; * ; MMHNO2=47.0; MMO3=48.0; MMO2=32.0; MMOH=17.0; MMHO2=33.0; MMH2O2=34.0; MMNO3=62.0; MMN2O5=108.0; MMHNO3=63.0; MMNH3=17.0; MMMO2=47.0; MMOP1=48.0; MMHCHO=30.0; MMORA1=46.0; MMCO2=44.0; MMHCL=36.0; MMCL2=70.0; MMBR2=160.0; MMBRCL=115.0; MMSO2=64.0; MMCLNO2=81.0; MMBRNO2=125.909; MMNO2=46.0; MMACO3=75.0 ; MMALD=44.0 ; MMCH3OH=32.0 ; MMETHP=61.0 ; CH3CH2OO MMETOH=46.0 ; MMOP2=61.0 ; CH3CH2OOH MMORA2=60.0 ; MMPAA=76.0 ; CH3C(O)OOH MMSULF=98.0 ; MMGLY=58.0; MMHNO4=79.; MMNO=30. ; * ; * ******************************************************************; * Diffusion coefficients; * ******************************************************************; * ; DGCO2 = 0.1550 ; DGHCl = 0.1890 ; DGNH3 = 0.2300 ; DGO3 = 0.1480 ; DGHO2 = 0.1040 ; DGOH = 0.1530 ; DGH2O2 = 0.1460 ; DGHNO3 = 0.132 ; DGNO3 = 0.1000 ; DGN2O5 = 0.1100 ; DGNO2 = 0.1920 ; DGHNO2 = 0.1300 ; DGHNO4 = 0.1300 ; DGSO2 = 0.1280 ; DGHCHO = 0.1640 ; DGOP1 = 0.1310 ; DGPAA = 0.1020 ; DGCH3OH = 0.116 ; DGETOH = 0.095 ; DGALD = 0.1220 ; DGORA1 = 0.1530 ; DGORA2 = 0.1240 ; DGMO2 = 0.1350 ; DGETHP = 0.1080 ; DGCl2 = 0.1280 ; DGBr2 = 0.1000 ; DGSULF = 0.1300 ; DGACO3 = 0.1000 ; DGGLY = 0.1150 ; DGO2 = 0.112 ; DGClNO2 = 0.1440 ; DGBrNO2 = 0.1440 ; DGBrCl = 0.1190 ; DGNO = 0.224 ; * ; * ******************************************************************; * Molecular speeds; * ******************************************************************; * ; VO3=100*SQRT(8*8.314*TEMP*1000/(MMO3*PI)); VO2=100*SQRT(8*8.314*TEMP*1000/(MMO2*PI)); VOH=100*SQRT(8*8.314*TEMP*1000/(MMOH*PI)); VHO2=100*SQRT(8*8.314*TEMP*1000/(MMHO2*PI)); VH2O2=100*SQRT(8*8.314*TEMP*1000/(MMH2O2*PI)); VNO3=100*SQRT(8*8.314*TEMP*1000/(MMNO3*PI)); VN2O5=100*SQRT(8*8.314*TEMP*1000/(MMN2O5*PI)); VHNO3=100*SQRT(8*8.314*TEMP*1000/(MMHNO3*PI)); VHNO4=100*SQRT(8*8.314*TEMP*1000/(MMHNO4*PI)); VHNO2=100*SQRT(8*8.314*TEMP*1000/(MMHNO2*PI)); VNH3=100*SQRT(8*8.314*TEMP*1000/(MMNH3*PI)); VMO2=100*SQRT(8*8.314*TEMP*1000/(MMMO2*PI)); VOP1=100*SQRT(8*8.314*TEMP*1000/(MMOP1*PI)); VHCHO=100*SQRT(8*8.314*TEMP*1000/(MMHCHO*PI)); VORA1=100*SQRT(8*8.314*TEMP*1000/(MMORA1*PI)); VCO2=100*SQRT(8*8.314*TEMP*1000/(MMCO2*PI)); VHCL=100*SQRT(8*8.314*TEMP*1000/(MMHCL*PI)); VCL2=100*SQRT(8*8.314*TEMP*1000/(MMCL2*PI)); VBR2=100*SQRT(8*8.314*TEMP*1000/(MMBR2*PI)); VBRCL=100*SQRT(8*8.314*TEMP*1000/(MMBRCL*PI)); VSO2=100*SQRT(8*8.314*TEMP*1000/(MMSO2*PI)); VBRNO2=100*SQRT(8*8.314*TEMP*1000/(MMBRNO2*PI)); VCLNO2=100*SQRT(8*8.314*TEMP*1000/(MMCLNO2*PI)); VNO2=100*SQRT(8*8.314*TEMP*1000/(MMNO2*PI)); VACO3=100*SQRT(8*8.314*TEMP*1000/(MMACO3*PI)); VALD=100*SQRT(8*8.314*TEMP*1000/(MMALD*PI)); VCH3OH=100*SQRT(8*8.314*TEMP*1000/(MMCH3OH*PI)); VETHP=100*SQRT(8*8.314*TEMP*1000/(MMETHP*PI)); VETOH=100*SQRT(8*8.314*TEMP*1000/(MMETOH*PI)); VOP2=100*SQRT(8*8.314*TEMP*1000/(MMOP2*PI)); VORA2=100*SQRT(8*8.314*TEMP*1000/(MMORA2*PI)); VPAA=100*SQRT(8*8.314*TEMP*1000/(MMPAA*PI)); VSULF=100*SQRT(8*8.314*TEMP*1000/(MMSULF*PI)); VGLY=100*SQRT(8*8.314*TEMP*1000/(MMGLY*PI)); VNO=100*SQRT(8*8.314*TEMP*1000/(MMNO*PI)); * *********; * *******************************************************************; * Emissions; * *******************************************************************; * ; * ; EMICH3OH = 1.9204e-13 /ZCL ; EMIETOH = 6.6928e-13 /ZCL ; EMIHC3 = 4.5870e-13 /ZCL ; EMIXYL = 1.6411e-13 /ZCL ; EMIKET = 1.6447e-13 /ZCL ; EMIHC5 = 8.1436e-13 /ZCL ; EMIETH = 2.5648e-13 /ZCL ; EMIOLT = 8.1948e-14 /ZCL ; EMIALD = 9.8507e-14 /ZCL ; EMIHCHO = 4.2862e-14 /ZCL ; EMITOL = 2.8301e-13 /ZCL ; EMIDIEN = 2.0664e-13 /ZCL ; EMICSL = 3.0284e-13 /ZCL ; EMIETE = 4.3388e-13 /ZCL ; EMIHC8 = 5.3026e-13 /ZCL ; EMICO = 1.4922e-11 /ZCL ; EMIORA2 = 1.4019e-14 /ZCL ; EMISO2 = 5.4241e-13 /ZCL ; EMINH3 = 5.0232e-13 /ZCL ; EMINO = 1.6744e-12 /ZCL ; EMIAPI = 3.26e-15 /ZCL ; EMILIM = 3.26e-15 /ZCL ; EMIISO = 2.55814e-14 /ZCL ; * ; * ***************************************************************; * Deposition velocities; * Urban scenario; * **************************************************************; * ; VDCO = 0.1 /ZCL ; VDHOBr = 0.2 /ZCL ; VDHOCl = 0.2 /ZCL ; VDNO2 = 0.4 /ZCL ; VDO3 = 0.4 /ZCL ; VDETOH = 0.5 /ZCL ; VDOP1 = 0.5 /ZCL ; VDCH3OH = 1 /ZCL ; VDH2O2 = 1 /ZCL ; VDHBr = 1 /ZCL ; VDHCHO = 1 /ZCL ; VDHCl = 1 /ZCL ; VDNH3 = 1 /ZCL ; VDORA1 = 1 /ZCL ; VDSO2 = 1 /ZCL ; VDHNO3 = 2 /ZCL ; VDN2O5 = 2 /ZCL ; VDSULF = 2 /ZCL ; * ; * *******************************************************************; * Initial concentrations; * *******************************************************************; * ; O2g=2.0E8*(1E-9)/(R*TEMP); N2g=7.8E8*(1E-9)/(R*TEMP); H2Og=2.0E7*(1E-9)/(R*TEMP); * *****************; * **************************************************************; NO2g = 4.5/(R*TEMP*1E9); HNO3g = 1/(R*TEMP*1E9) ; CH4g = 1700/(R*TEMP*1E9) ; H2O2g = 1/(R*TEMP*1E9) ; H2g = 500/(R*TEMP*1E9) ; COg = 300/(R*TEMP*1E9) ; O3g = 90/(R*TEMP*1E9) ; HClg = 0.2/(R*TEMP*1E9) ; NH3g = 5/(R*TEMP*1E9) ; CO2g = 3.57e5/(R*TEMP*1E9) ; SO2g = 5/(R*TEMP*1E9) ; HCHOg = 1/(R*TEMP*1E9) ; ETHg = 2/(R*TEMP*1E9) ; HC3g = 2/(R*TEMP*1E9) ; HC5g = 1/(R*TEMP*1E9) ; HC8g = 0.1/(R*TEMP*1E9) ; OLTg = 0.1/(R*TEMP*1E9) ; Isog = 0/(R*TEMP*1E9) ; TOLg = 0.1/(R*TEMP*1E9) ; CSLg = 0.001/(R*TEMP*1E9) ; XYLg = 0.1/(R*TEMP*1E9) ; ALDg = 0.1/(R*TEMP*1E9) ; Ketg = 1/(R*TEMP*1E9) ; Glyg = 0.1/(R*TEMP*1E9) ; Mglyg = 0.1/(R*TEMP*1E9) ; PANg = 0.5/(R*TEMP*1E9) ; OP1g = 1/(R*TEMP*1E9) ; OP2g = 0.1/(R*TEMP*1E9) ; PAAg = 0.001/(R*TEMP*1E9) ; CH3OHg = 5/(R*TEMP*1E9) ; EtOHg = 1/(R*TEMP*1E9) ; APIg = 0/(R*TEMP*1E9) ; LIMg = 0/(R*TEMP*1E9) ; * Reaction rates ; K24 = ( ( CM*6.00E-34*(TEMP/300)@(-2.3))) * CONV2 ; K25 = ( ( 8.00E-12*EXP( -2060.0/TEMP))) * CONV2 ; K26 = ( ( 1.80E-11*EXP( 110.0/TEMP)))* CONV2 ; K27 = ( ( 3.20E-11*EXP( 70.0/TEMP))) * CONV2 ; K28 = ( ( 2.20E-10)) * CONV2 ; K29 = ( ( 1.60E-12*EXP( -940.0/TEMP))) * CONV2 ; K30 = ( ( 1.10E-14*EXP( -500.0/TEMP))) * CONV2 ; K31 = ( ( 4.80E-11*EXP( 250.0/TEMP))) * CONV2 ; K32 = ( ( 2.90E-12*EXP( -160.0/TEMP))) * CONV2 ; K33 = ( ( 2.3E-13*EXP(600./TEMP) +1.7E-33 * CM * EXP(1000./TEMP))) * CONV2 ; K34 = ( ( 3.22E-34 * EXP(2800./TEMP) +2.38E-54*CM*EXP(3200./TEMP))) * CONV2* CONV2 ; k0t = 9.00E-32*(TEMP/300.0)@(-1.5); kinft = 3.00E-11; K35 = ( ( (k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2)))) * CONV2 ; K36 = ( ( 6.50E-12*EXP( 120.0/TEMP))) * CONV2 ; k0t = 9.00E-32*(TEMP/300.0)@(-2.0); kinft = 2.20E-11; K37 = ( ( (k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2)))) * CONV2 ; k0t = 7.00E-31*(TEMP/300.0)@(-2.6); kinft = 1.50E-11*(TEMP/300.0)@(-0.5); K38 = ( ( (k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2)))) * CONV2 ; k0t = 2.60E-30*(TEMP/300.0)@(-3.2); kinft = 2.40E-11*(TEMP/300.0)@(-1.3); K39 = ( ( (k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2)))) * CONV2 ; K40 = ( ( 2.20E-11)) * CONV2 ; K41 = ( ( 3.70E-12*EXP( 250.0/TEMP))) * CONV2 ; k0t = 1.80E-31*(TEMP/300.0)@(-3.2); kinft = 4.70E-12*(TEMP/300.0)@(-1.4); K42 = ( ( (k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2)))) * CONV2 ; k0t = 1.80E-31*(TEMP/300.0)@(-3.2); kinft = 4.70E-12*(TEMP/300.0)@(-1.4); K43 = ( 4.76E+26*EXP( -10900.0/TEMP)*(k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2))) ; K44 = ( ( 3.50E-12)) * CONV2 ; K45 = ( ( 1.80E-11*EXP( -390.0/TEMP))) * CONV2 ; TMP0 = 7.2E-15*EXP(785/TEMP); TMP2 = 4.1E-16*EXP(1440/TEMP); TMP3 = 1.9E-33*EXP(725/TEMP)*CM; K46 = ( ( TMP0 + TMP3/(1+TMP3/TMP2))) * CONV2 ; K47 = ( ( 1.30E-12*EXP( 380.0/TEMP))) * CONV2 ; K48 = ( ( 2.00E-12*EXP( -1400.0/TEMP))) * CONV2 ; K49 = ( ( 1.20E-13*EXP( -2450.0/TEMP))) * CONV2 ; K50 = ( ( 3.30E-39*EXP( 530.0/TEMP))) * CONV2* CONV2 ; K51 = ( ( 1.50E-11*EXP( 170.0/TEMP))) * CONV2 ; K52 = ( ( 4.50E-14*EXP( -1260.0/TEMP))) * CONV2 ; k0t = 2.20E-30*(TEMP/300.0)@(-3.9); kinft = 1.50E-12*(TEMP/300.0)@(-0.7); K53 = ( ( (k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2)))) * CONV2 ; k0t = 2.20E-30*(TEMP/300.0)@(-3.9); kinft = 1.50E-12*(TEMP/300.0)@(-0.7); K54 = ( 3.70E+26*EXP( -11000.0/TEMP)*(k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2))); K55 = ( ( 8.50E-13*EXP( -2450.0/TEMP))) * CONV2 ; K56 = ( ( 5.50E-12*EXP( -2000.0/TEMP))) * CONV2 ; k0t = 3.00E-31*(TEMP/300.0)@(-3.3); kinft = 1.50E-12; K57 = ( ( (k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2)))) * CONV2 ; K58 = ( ( 1.5E-13 * (1. + 2.439E-20 * CM))) * CONV2 ; K59 = ( ( 6.00E-11)) * CONV2 ; K60 = ( ( 1.59E-11*EXP( 13.0/TEMP))) * CONV2 ; K61 = ( ( 7.44E-18 *(TEMP**2)*EXP( -1361.0/TEMP))) * CONV2 ; K62 = ( ( 1.51E-17 *(TEMP**2)*EXP( -492.0/TEMP))) * CONV2 ; K63 = ( ( 5.26E-12*EXP( -260.0/TEMP))) * CONV2 ; K64 = ( ( 8.02E-12*EXP( -155.0/TEMP))) * CONV2 ; K65 = ( ( 1.64E-11*EXP( -125.0/TEMP))) * CONV2 ; K66 = ( ( 1.96E-12*EXP( 438.0/TEMP))) * CONV2 ; K67 = ( ( 5.72E-12*EXP( 500.0/TEMP))) * CONV2 ; K68 = ( ( 1.33E-11*EXP( 500.0/TEMP))) * CONV2 ; K69 = ( ( 1.48E-11*EXP( 448.0/TEMP))) * CONV2 ; K70 = ( ( 2.54E-11*EXP( 410.0/TEMP))) * CONV2 ; K71 = ( ( 1.21E-11*EXP( 444.0/TEMP))) * CONV2 ; K72 = ( ( 1.71E-10)) * CONV2 ; K73 = ( ( 1.81E-12*EXP( 355.0/TEMP))) * CONV2 ; K74 = ( ( 7.30E-12*EXP( 355.0/TEMP))) * CONV2 ; K75 = ( ( 6.00E-11)) * CONV2 ; K76 = ( ( 1.00E-11)) * CONV2 ; K77 = ( ( 5.55E-12*EXP( 331.0/TEMP))) * CONV2 ; K78 = ( ( 5.68E-18 *(TEMP**2)*EXP( 92.0/TEMP))) * CONV2 ; K79 = ( ( 3.00E-12)) * CONV2 ; K80 = ( ( 1.14E-11)) * CONV2 ; K81 = ( ( 1.72E-11)) * CONV2 ; K82 = ( ( 1.86E-11*EXP( 175.0/TEMP))) * CONV2 ; K83 = ( ( 2.80E-11*EXP( 175.0/TEMP))) * CONV2 ; K84 = ( ( 2.70E-10)) * CONV2 ; K85 = ( ( 2.93E-12*EXP( 190.0/TEMP))) * CONV2 ; K86 = ( ( 3.40E-12*EXP( 190.0/TEMP))) * CONV2 ; K87 = ( ( 2.93E-12*EXP( 190.0/TEMP))) * CONV2 ; K88 = ( ( 4.00E-14)) * CONV2 ; K89 = ( ( 3.25E-13*EXP( 500.0/TEMP))) * CONV2 ; K90 = ( ( 5.31E-12*EXP( -260.0/TEMP))) * CONV2 ; K91 = ( ( 3.40E-13*EXP( -1900.0/TEMP))) * CONV2 ; K92 = ( ( 1.40E-12*EXP( -1900.0/TEMP))) * CONV2 ; K93 = ( ( 2.90E-12*EXP( -1900.0/TEMP))) * CONV2 ; K94 = ( ( 1.40E-12*EXP( -1900.0/TEMP))) * CONV2 ; K95 = ( ( 8.27E-15*EXP( -150.0/TEMP))) * CONV2 ; K96 = ( ( 2.87E-13*EXP( -1000.0/TEMP))) * CONV2 ; K97 = ( ( 2.20E-11)) * CONV2 ; K98 = ( ( 4.88E-18 *(TEMP**2)*EXP( -2282.0/TEMP))) * CONV2; K99 = ( ( 1.79E-13*EXP( -450.0/TEMP))) * CONV2 ; K100 = ( ( 8.64E-13*EXP( 450.0/TEMP))) * CONV2 ; K101 = ( ( 1.00E-13)) * CONV2 ; K102 = ( ( 4.00E-12*EXP( -446.0/TEMP))) * CONV2 ; K103 = ( ( 1.19E-12*EXP( 490.0/TEMP))) * CONV2 ; K104 = ( ( 1.22E-11)) * CONV2 ; K105 = ( ( 2.20E-14*EXP( -500.0/TEMP))) * CONV2 ; K106 = ( ( 9.14E-15*EXP( -2580.0/TEMP))) * CONV2 ; K107 = ( ( 4.33E-15*EXP( -1800.0/TEMP))) * CONV2 ; K108 = ( ( 4.40E-15*EXP( -845.0/TEMP))) * CONV2 ; K109 = ( ( 1.34E-14*EXP( -2283.0/TEMP))) * CONV2 ; K110 = ( ( 7.86E-15*EXP( -1913.0/TEMP))) * CONV2 ; K111 = ( ( 1.01E-15*EXP( -736.0/TEMP))) * CONV2 ; K112 = ( ( 2.00E-16)) * CONV2 ; K113 = ( ( 1.36E-15*EXP( -2112.0/TEMP))) * CONV2 ; K114 = ( ( 2.00E-18)) * CONV2 ; K115 = ( ( 2.46E-15*EXP( -1700.0/TEMP))) * CONV2 ; K116 = ( ( 2.00E-11)) * CONV2 ; K117 = ( ( 1.00E-11)) * CONV2 ; K118 = ( ( 3.60E-11)) * CONV2 ; K119 = ( ( 1.66E-17*EXP( 1044.0/TEMP))) * CONV2 ; K120 = ( ( 5.00E-11)) * CONV2 ; K121 = ( ( 3.60E-11)) * CONV2 ; K122 = ( ( 1.66E-17*EXP( 1044.0/TEMP))) * CONV2 ; K123 = ( ( 1.00E-11)) * CONV2 ; K124 = ( ( 3.60E-11)) * CONV2 ; K125 = ( ( 1.66E-17*EXP( 1044.0/TEMP))) * CONV2 ; K126 = ( ( 5.00E-11)) * CONV2 ; k0t = 9.70E-29*(TEMP/300.0)@(-5.6); kinft = 9.30E-12*(TEMP/300.0)@(-1.5); K127 = ( ( ((k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2))))) * CONV2 ; k0t = 9.70E-29*(TEMP/300.0)@(-5.6); kinft = 9.30E-12*(TEMP/300.0)@(-1.5); K128 = ( (1.16E+28*EXP( -13954.0/TEMP)*(k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2)))) ; k0t = 9.70E-29*(TEMP/300.0)@(-5.6); kinft = 9.30E-12*(TEMP/300.0)@(-1.5); K129 = ( ( (k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2)))) * CONV2 ; k0t = 9.70E-29*(TEMP/300.0)@(-5.6); kinft = 9.30E-12*(TEMP/300.0)@(-1.5); K130 = ( 1.16E+28*EXP( -13954.0/TEMP)*(k0t*CM/(1+k0t*CM/kinft)) *0.6@(1./(1.+(log10(k0t*CM/kinft))**2))) ; K131 = ( ( 4.20E-12*EXP( 180.0/TEMP))) * CONV2 ; K132 = ( ( 8.70E-12)) * CONV2 ; K133 = ( ( 4.00E-12)) * CONV2 ; K134 = ( ( 4.00E-12)) * CONV2 ; K135 = ( ( 4.00E-12)) * CONV2 ; K136 = ( ( 9.00E-12)) * CONV2 ; K137 = ( ( 4.00E-12)) * CONV2 ; K138 = ( ( 4.00E-12)) * CONV2 ; K139 = ( ( 4.00E-12)) * CONV2 ; K140 = ( ( 4.00E-12)) * CONV2 ; K141 = ( ( 4.00E-12)) * CONV2 ; K142 = ( ( 4.00E-12)) * CONV2 ; K143 = ( ( 4.00E-12)) * CONV2 ; K144 = ( ( 4.00E-12)) * CONV2 ; K145 = ( ( 2.00E-11)) * CONV2 ; K146 = ( ( 2.00E-11)) * CONV2 ; K147 = ( ( 4.00E-12)) * CONV2 ; K148 = ( ( 4.00E-12)) * CONV2 ; K149 = ( ( 4.00E-12)) * CONV2 ; K150 = ( ( 3.80E-13*EXP( 800.0/TEMP))) * CONV2 ; K151 = ( ( 7.50E-13*EXP( 700.0/TEMP))) * CONV2 ; K152 = ( ( 1.66E-13*EXP( 1300.0/TEMP))) * CONV2 ; K153 = ( ( 1.66E-13*EXP( 1300.0/TEMP))) * CONV2 ; K154 = ( ( 1.66E-13*EXP( 1300.0/TEMP))) * CONV2 ; K155 = ( ( 1.90E-13*EXP( 1300.0/TEMP))) * CONV2 ; K156 = ( ( 1.66E-13*EXP( 1300.0/TEMP))) * CONV2 ; K157 = ( ( 1.66E-13*EXP( 1300.0/TEMP))) * CONV2 ; K158 = ( ( 1.28E-13*EXP( 1300.0/TEMP))) * CONV2 ; K159 = ( ( 1.50E-11)) * CONV2 ; K160 = ( ( 1.50E-11)) * CONV2 ; K161 = ( ( 3.75E-13*EXP( 980.0/TEMP))) * CONV2 ; K162 = ( ( 3.75E-13*EXP( 980.0/TEMP))) * CONV2 ; K163 = ( ( 3.75E-13*EXP( 980.0/TEMP))) * CONV2 ; K164 = ( ( 1.15E-12*EXP( 550.0/TEMP))) * CONV2 ; K165 = ( ( 3.86E-16*EXP( 2640.0/TEMP))) * CONV2 ; K166 = ( ( 1.15E-12*EXP( 550.0/TEMP))) * CONV2 ; K167 = ( ( 3.86E-16*EXP( 2640.0/TEMP))) * CONV2 ; K168 = ( ( 1.15E-13*EXP( 1300.0/TEMP))) * CONV2 ; K169 = ( ( 1.66E-13*EXP( 1300.0/TEMP))) * CONV2 ; K170 = ( ( 1.66E-13*EXP( 1300.0/TEMP))) * CONV2 ; K171 = ( ( 9.10E-14*EXP( 416.0/TEMP))) * CONV2 ; K172 = ( ( 1.18E-13*EXP( 158.0/TEMP))) * CONV2 ; K173 = ( ( 9.46E-14*EXP( 431.0/TEMP))) * CONV2 ; K174 = ( ( 1.00E-13*EXP( 467.0/TEMP))) * CONV2 ; K175 = ( ( 4.34E-14*EXP( 633.0/TEMP))) * CONV2 ; K176 = ( ( 1.71E-13*EXP( 708.0/TEMP))) * CONV2 ; K177 = ( ( 1.46E-13*EXP( 708.0/TEMP))) * CONV2 ; K178 = ( ( 9.18E-14*EXP( 708.0/TEMP))) * CONV2 ; K179 = ( ( 1.36E-13*EXP( 708.0/TEMP))) * CONV2 ; K180 = ( ( 3.56E-14*EXP( 708.0/TEMP))) * CONV2 ; K181 = ( ( 3.56E-14*EXP( 708.0/TEMP))) * CONV2 ; K182 = ( ( 3.56E-14*EXP( 708.0/TEMP))) * CONV2 ; K183 = ( ( 3.56E-14*EXP( 708.0/TEMP))) * CONV2 ; K184 = ( ( 3.56E-14*EXP( 708.0/TEMP))) * CONV2 ; K185 = ( ( 3.21E-11*EXP( -440.0/TEMP))) * CONV2 ; K186 = ( ( 2.68E-16*EXP( 2510.0/TEMP))) * CONV2 ; K187 = ( ( 3.21E-11*EXP( -440.0/TEMP))) * CONV2 ; K188 = ( ( 2.68E-16*EXP( 2510.0/TEMP))) * CONV2 ; K189 = ( ( 6.91E-13*EXP( 508.0/TEMP))) * CONV2 ; K190 = ( ( 1.60E-13*EXP( 708.0/TEMP))) * CONV2 ; K191 = ( ( 9.68E-14*EXP( 708.0/TEMP))) * CONV2 ; K192 = ( ( 1.03E-12*EXP( 211.0/TEMP))) * CONV2 ; K193 = ( ( 6.90E-13*EXP( 460.0/TEMP))) * CONV2 ; K194 = ( ( 5.59E-13*EXP( 522.0/TEMP))) * CONV2 ; K195 = ( ( 2.47E-13*EXP( 683.0/TEMP))) * CONV2 ; K196 = ( ( 9.48E-13*EXP( 765.0/TEMP))) * CONV2 ; K197 = ( ( 8.11E-13*EXP( 765.0/TEMP))) * CONV2 ; K198 = ( ( 5.09E-13*EXP( 765.0/TEMP))) * CONV2 ; K199 = ( ( 7.60E-13*EXP( 765.0/TEMP))) * CONV2 ; K200 = ( ( 7.40E-13*EXP( 765.0/TEMP))) * CONV2 ; K201 = ( ( 7.40E-13*EXP( 765.0/TEMP))) * CONV2 ; K202 = ( ( 7.40E-13*EXP( 765.0/TEMP))) * CONV2 ; K203 = ( ( 7.40E-13*EXP( 765.0/TEMP))) * CONV2 ; K204 = ( ( 7.40E-13*EXP( 765.0/TEMP))) * CONV2 ; K205 = ( ( 2.80E-12*EXP( 530.0/TEMP))) * CONV2 ; K206 = ( ( 2.80E-12*EXP( 530.0/TEMP))) * CONV2 ; K207 = ( ( 7.51E-13*EXP( 565.0/TEMP))) * CONV2 ; K208 = ( ( 8.85E-13*EXP( 765.0/TEMP))) * CONV2 ; K209 = ( ( 5.37E-13*EXP( 765.0/TEMP))) * CONV2 ; K210 = ( ( 7.00E-14*EXP( 1000.0/TEMP))) * CONV2 ; K211 = ( ( 4.25E-14*EXP( 1000.0/TEMP))) * CONV2 ; K212 = ( ( 2.96E-14*EXP( 1000.0/TEMP))) * CONV2 ; K213 = ( ( 1.20E-12)) * CONV2 ; K214 = ( ( 1.20E-12)) * CONV2 ; K215 = ( ( 1.20E-12)) * CONV2 ; K216 = ( ( 1.20E-12)) * CONV2 ; K217 = ( ( 1.20E-12)) * CONV2 ; K218 = ( ( 1.20E-12)) * CONV2 ; K219 = ( ( 1.20E-12)) * CONV2 ; K220 = ( ( 1.20E-12)) * CONV2 ; K221 = ( ( 1.20E-12)) * CONV2 ; K222 = ( ( 1.20E-12)) * CONV2 ; K223 = ( ( 1.20E-12)) * CONV2 ; K224 = ( ( 1.20E-12)) * CONV2 ; K225 = ( ( 1.20E-12)) * CONV2 ; K226 = ( ( 1.20E-12)) * CONV2 ; K227 = ( ( 4.00E-12)) * CONV2 ; K228 = ( ( 4.00E-12)) * CONV2 ; K229 = ( ( 1.20E-12)) * CONV2 ; K230 = ( ( 1.20E-12)) * CONV2 ; K231 = ( ( 1.20E-12)) * CONV2 ; K232 = ( ( 1.66E-13*EXP( 1300.0/TEMP))) * CONV2 ; K233 = ( ( 5.99E-15*EXP( 1510.0/TEMP))) * CONV2 ; K234 = ( ( 3.40E-14*EXP( 1560.0/TEMP))) * CONV2 ; K235 = ( ( 7.13E-17*EXP( 2950.0/TEMP))) * CONV2 ; K236 = ( ( 4.00E-12)) * CONV2 ; K237 = ( ( 1.20E-12)) * CONV2 ; H2Oa=0.0556; OHm=3.16e-10/CONV3; Hp=3.16e-05/CONV3; ** ; COMPILE HETERO; HET=1; KMTCO2 =HET/(((RAD**2)/(3* DGCO2 ))+(4*RAD/(3* VCO2 *ACCO2 ))); KMTHCl =HET/(((RAD**2)/(3* DGHCl ))+(4*RAD/(3* VHCl *ACHCl ))); KMTNH3 =HET/(((RAD**2)/(3* DGNH3 ))+(4*RAD/(3* VNH3 *ACNH3 ))); KMTO3 =HET/(((RAD**2)/(3* DGO3 ))+(4*RAD/(3* VO3 *ACO3 ))); KMTHO2 =HET/(((RAD**2)/(3* DGHO2 ))+(4*RAD/(3* VHO2 *ACHO2 ))); KMTOH =HET/(((RAD**2)/(3* DGOH ))+(4*RAD/(3* VOH *ACOH ))); KMTH2O2 =HET/(((RAD**2)/(3* DGH2O2 ))+(4*RAD/(3* VH2O2 *ACH2O2 ))); KMTHNO3 =HET/(((RAD**2)/(3* DGHNO3 ))+(4*RAD/(3* VHNO3 *ACHNO3 ))); KMTNO3 =HET/(((RAD**2)/(3* DGNO3 ))+(4*RAD/(3* VNO3 *ACNO3 ))); KMTN2O5 =HET/(((RAD**2)/(3* DGN2O5 ))+(4*RAD/(3* VN2O5 *ACN2O5 ))); KMTNO2 =HET/(((RAD**2)/(3* DGNO2 ))+(4*RAD/(3* VNO2 *ACNO2 ))); KMTHNO2 =HET/(((RAD**2)/(3* DGHNO2 ))+(4*RAD/(3* VHNO2 *ACHNO2 ))); KMTHNO4 =HET/(((RAD**2)/(3* DGHNO4 ))+(4*RAD/(3* VHNO4 *ACHNO4 ))); KMTSO2 =HET/(((RAD**2)/(3* DGSO2 ))+(4*RAD/(3* VSO2 *ACSO2 ))); KMTHCHO =HET/(((RAD**2)/(3* DGHCHO ))+(4*RAD/(3* VHCHO *ACHCHO ))); KMTOP1 =HET/(((RAD**2)/(3* DGOP1 ))+(4*RAD/(3* VOP1 *ACOP1 ))); KMTPAA =HET/(((RAD**2)/(3* DGPAA ))+(4*RAD/(3* VPAA *ACPAA ))); KMTCH3OH =HET/(((RAD**2)/(3*DGCH3OH))+(4*RAD/(3*VCH3OH*ACCH3OH))); KMTETOH =HET/(((RAD**2)/(3* DGETOH ))+(4*RAD/(3* VETOH *ACETOH ))); KMTALD =HET/(((RAD**2)/(3* DGALD ))+(4*RAD/(3* VALD *ACALD ))); KMTORA1 =HET/(((RAD**2)/(3* DGORA1 ))+(4*RAD/(3* VORA1 *ACORA1 ))); KMTORA2 =HET/(((RAD**2)/(3* DGORA2 ))+(4*RAD/(3* VORA2 *ACORA2 ))); KMTMO2 =HET/(((RAD**2)/(3* DGMO2 ))+(4*RAD/(3* VMO2 *ACMO2 ))); KMTETHP =HET/(((RAD**2)/(3* DGETHP ))+(4*RAD/(3* VETHP *ACETHP ))); KMTCl2 =HET/(((RAD**2)/(3* DGCl2 ))+(4*RAD/(3* VCl2 *ACCl2 ))); KMTBr2 =HET/(((RAD**2)/(3* DGBr2 ))+(4*RAD/(3* VBr2 *ACBr2 ))); KMTSULF =HET/(((RAD**2)/(3*DGSULF ))+(4*RAD/(3* VSULF *ACSULF ))); KMTACO3 =HET/(((RAD**2)/(3* DGACO3 ))+(4*RAD/(3* VACO3 *ACACO3 ))); KMTGLY =HET/(((RAD**2)/(3* DGGLY ))+(4*RAD/(3* VGLY *ACGLY ))); KMTO2 =HET/(((RAD**2)/(3* DGO2 ))+(4*RAD/(3* VO2 *ACO2 ))); KMTClNO2 =HET/(((RAD**2)/(3*DGClNO2))+(4*RAD/(3*VClNO2*ACClNO2))); KMTBrNO2 =HET/(((RAD**2)/(3*DGBrNO2))+(4*RAD/(3*VBrNO2*ACBrNO2))); KMTBrCl =HET/(((RAD**2)/(3* DGBrCl ))+(4*RAD/(3* VBrCl *ACBrCl ))); KMTNO =HET/(((RAD**2)/(3* DGNO ))+(4*RAD/(3* VNO *ACNO ))); * *********; * ; * Kinetics in aqueous solution Capram 2.4; * ; * ; KCAP1= HET * 50 *Exp( -TCORR * 0 ) * CONV3 ; KCAP2= HET * 7.3E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP3= HET * 2.8E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP4= HET * 1E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP5= HET * 7.0E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP6= HET * 1.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP7= HET * 1.3E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP8= HET * 1.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP9= HET * 1.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP10= HET * 1.0E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP11= HET * 1.2E6 *Exp( -TCORR * 5050 ) * CONV3 ; KCAP12= HET * 4.3E8 *Exp( -TCORR * 1100 ) * CONV3 ; KCAP13= HET * 1E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP14= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP15= HET * 1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP16= HET * 2.3E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP17= HET * 3E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP18= HET * 1E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP19= HET * 8E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP20= HET * 1.3E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP21= HET * 1.3E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP22= HET * 1.3E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP23= HET * 1.6E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP24= HET * 1.5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP25= HET * 330 *Exp( -TCORR * 0 ) ; KCAP26= HET * 4.15E5 *Exp( -TCORR * 2720 ) * CONV3 ; KCAP27= HET * 4.85E7 *Exp( -TCORR * 1060 ) * CONV3 ; KCAP28= HET * 1.0E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP29= HET * 1.1E10 *Exp( -TCORR * 2120 ) * CONV3 ; KCAP30= HET * 3.0E7 *Exp( -TCORR * 1680 ) * CONV3 ; KCAP31= HET * 2.4E7 *Exp( -TCORR * 1680 ) * CONV3 ; KCAP32= HET * 6.0E6 *Exp( -TCORR * 1680 ) * CONV3 ; KCAP33= HET * 2.7E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP34= HET * 4.6E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP35= HET * 6E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP36= HET * 4.6E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP37= HET * 8.2E5 *Exp( -TCORR * 4690.0 ) * CONV3 ; KCAP38= HET * 9.5E3 *Exp( -TCORR * 2766.0 ) * CONV3 ; KCAP39= HET * 2.34E-4 *Exp( -TCORR * 4089.0 ) * CONV3; KCAP40= HET * 2.0E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP41= HET * 1.0E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP42= HET * 1.1E4 *Exp( -TCORR * 4150 ) * CONV3 ; KCAP43= HET * 2.5E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP44= HET * 100 *Exp( -TCORR * 0 ) * CONV3 ; KCAP45= HET * 1.0E4 *Exp( -TCORR * 2562 ) * CONV3 ; KCAP46= HET * 160 *Exp( -TCORR * 2680 ) * CONV3 ; KCAP47= HET * 3E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP48= HET * 2.5E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP49= HET * 400 *Exp( -TCORR * 5352 ) * CONV3 ; KCAP50= HET * 1E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP51= HET * 7.2E4 *Exp( -TCORR * 842.0 ) * CONV3 ; KCAP52= HET * 1.8E4 *Exp( -TCORR * 5052.0 ) * CONV3 ; ; KCAP53= HET * 1.95 *Exp( -TCORR * 5653 ) * CONV3 ; KCAP54= HET * 0.49 *Exp( -TCORR * 8780.0 ) ; KCAP55= HET * 1E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP56= HET * 1.65E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP57= HET * 3E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP58= HET * 5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP59= HET * 1.0E9 *Exp( -TCORR * 0.0 ) ; KCAP60= HET * 8.9E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP61= HET * 9.4E7 *Exp( -TCORR * 2700 ) * CONV3 ; KCAP62= HET * 8.0E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP63= HET * 1.1E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP64= HET * 4.9E6 *Exp( -TCORR * 2000 ) * CONV3 ; KCAP65= HET * 4.9E6 *Exp( -TCORR * 2000 ) * CONV3 ; KCAP66= HET * 3.0E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP67= HET * 3E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP68= HET * 1.3E9 *Exp( -TCORR * 2000 ) * CONV3 ; KCAP69= HET * 3.0E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP70= HET * 2.6E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP71= HET * 1.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP72= HET * 1.2E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP73= HET * 1E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP74= HET * 8.4E7 *Exp( -TCORR * (-2900) ) * CONV3 ; KCAP75= HET * 4.5E-2 *Exp( -TCORR * 0 ) ; KCAP76= HET * 9.1E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP77= HET * 7.2E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP78= HET * 1.4E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP79= HET * 6E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP80= HET * 1.2E7 *Exp( -TCORR * 1720 ) * CONV3 ; KCAP81= HET * 6.6E5 *Exp( -TCORR * 850 ) * CONV3 ; KCAP82= HET * 5E5 *Exp( -TCORR * 6900 ) * CONV3 ; KCAP83= HET * 1.1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP84= HET * 3.3E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP85= HET * 1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP86= HET * 1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP87= HET * 5E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP88= HET * 2.55E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP89= HET * 10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP90= HET * 3E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP91= HET * 1.3E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP92= HET * 4.2E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP93= HET * 5.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP94= HET * 5E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP95= HET * 1E5 *Exp( -TCORR * 0 ) ; KCAP96= HET * 1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP97= HET * 2.6E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP98= HET * 1.7E4 *Exp( -TCORR * 0 ) ; KCAP99= HET * 7.0E3 *Exp( -TCORR * 0 ) ; KCAP100= HET * 44.32 *Exp( -TCORR * 0 ) * CONV3 ; KCAP101= HET * 1.6E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP102= HET * 7.2E7 *Exp( -TCORR * 4000 ) * CONV3 * CONV3 ; KCAP103= HET * 1.7E7 *Exp( -TCORR * 3800 ) * CONV3 * CONV3 ; KCAP104= HET * 5.6E7 *Exp( -TCORR * 3990 ) * CONV3 * CONV3 ; KCAP105= HET * 2.4E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP106= HET * 3.7E5 *Exp( -TCORR * 5530 ) * CONV3 ; KCAP107= HET * 1.5E9 *Exp( -TCORR * 5280 ) * CONV3 ; KCAP108= HET * 30 *Exp( -TCORR * 0 ) * CONV3 ; KCAP109= HET * 2.65E7 *Exp( -TCORR * 5809 ) * CONV3 ; KCAP110= HET * 3E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP111= HET * 1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP112= HET * 4.6E9 *Exp( -TCORR * (-2165) ) * CONV3 ; KCAP113= HET * 17 *Exp( -TCORR * 0 ) * CONV3 ; KCAP114= HET * 4.8E7 *Exp( -TCORR * 2600 ) * CONV3 ; KCAP115= HET * 2.2E8 *Exp( -TCORR * 2600 ) * CONV3 ; KCAP116= HET * 1.7E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP117= HET * 2.34E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP118= HET * 3.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP119= HET * 5.2E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP120= HET * 1200 *Exp( -TCORR * 0 ) ; KCAP121= HET * 2.5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP122= HET * 8.6E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP123= HET * 3.6E2 *Exp( -TCORR * 0 ) * CONV3 ; KCAP124= HET * 2.13E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP125= HET * 5.5E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP126= HET * 3.5E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP127= HET * 6.1E8 *Exp( -TCORR * 840 ) * CONV3 ; KCAP128= HET * 5.8E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP129= HET * 3.4E8 *Exp( -TCORR * 1200 ) * CONV3 ; KCAP130= HET * 1.4E7 *Exp( -TCORR * 4089 ) * CONV3 ; KCAP131= HET * 3E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP132= HET * 1.7E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP133= HET * 1.7E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP134= HET * 3.5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP135= HET * 3.5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP136= HET * 5.0E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP137= HET * 1.4E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP138= HET * 11 *Exp( -TCORR * 1110 ) * CONV3 ; KCAP139= HET * 7.14E6 *Exp( -TCORR * 0 ) * CONV3 * CONV3 ; KCAP140= HET * 7.14E6 *Exp( -TCORR * 0 ) * CONV3 * CONV3 ; KCAP141= HET * 1.7E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP142= HET * 1.0E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP143= HET * 1.3E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP144= HET * 1.0E9 *Exp( -TCORR * 580 ) * CONV3 ; KCAP145= HET * 9.0E6 *Exp( -TCORR * 2190 ) * CONV3 ; KCAP146= HET * 5.4E5 *Exp( -TCORR * 4300 ) * CONV3 ; KCAP147= HET * 5.1E4 *Exp( -TCORR * 5500 ) * CONV3 ; KCAP148= HET * 1000 *Exp( -TCORR * 0 ) * CONV3 ; KCAP149= HET * 2.6E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP150= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP151= HET * 1.65E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP152= HET * 1.05E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP153= HET * 1.9E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP154= HET * 4.1E7 *Exp( -TCORR * 1760 ) * CONV3 ; KCAP155= HET * 2.2E6 *Exp( -TCORR * 3300 ) * CONV3 ; KCAP156= HET * 1.2E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP157= HET * 3.8E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP158= HET * 1.5E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP159= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP160= HET * 52 *Exp( -TCORR * 7217 ) ; KCAP161= HET * 8E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP162= HET * 1.0E9 *Exp( -TCORR * 1020 ) * CONV3 ; KCAP163= HET * 1.4E7 *Exp( -TCORR * 1300 ) * CONV3 ; KCAP164= HET * 1.0E6 *Exp( -TCORR * 4500 ) * CONV3 ; KCAP165= HET * 3.6E4 *Exp( -TCORR * 4330 ) * CONV3 ; KCAP166= HET * 3E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP167= HET * 1.3E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP168= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP169= HET * 1.2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP170= HET * 3.6E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP171= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP172= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP173= HET * 1.9E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP174= HET * 1.9E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP175= HET * 4E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP176= HET * 4E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP177= HET * 4E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP178= HET * 4E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP179= HET * 1E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP180= HET * 1E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP181= HET * 1.3E8 *Exp( -TCORR * 1000 ) * CONV3 ; KCAP182= HET * 3.2E9 *Exp( -TCORR * 1000 ) * CONV3 ; KCAP183= HET * 2.5E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP184= HET * 2.1E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP185= HET * 3.8E5 *Exp( -TCORR * 3400 ) * CONV3 ; KCAP186= HET * 5.1E7 *Exp( -TCORR * 2200 ) * CONV3 ; KCAP187= HET * 8E4 *Exp( -TCORR * 4450 ) * CONV3 ; KCAP188= HET * 1.3E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP189= HET * 4E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP190= HET * 4.9E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP191= HET * 1.4E5 *Exp( -TCORR * 3300 ) * CONV3 ; KCAP192= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP193= HET * 1.5E7 *Exp( -TCORR * 1330 ) * CONV3 ; KCAP194= HET * 1.0E8 *Exp( -TCORR * 1800 ) * CONV3 ; KCAP195= HET * 2.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP196= HET * 2.8E7 *Exp( -TCORR * 1210 ) * CONV3 ; KCAP197= HET * 1.4E4 *Exp( -TCORR * 3800 ) * CONV3 ; KCAP198= HET * 2.9E6 *Exp( -TCORR * 3800 ) * CONV3 ; KCAP199= HET * 1500 *Exp( -TCORR * 4930 ) * CONV3 ; KCAP200= HET * 2.6E5 *Exp( -TCORR * 4800 ) * CONV3 ; KCAP201= HET * 10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP202= HET * 1.7E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP203= HET * 100 *Exp( -TCORR * 0 ) * CONV3 ; KCAP204= HET * 580 *Exp( -TCORR * 0 ) * CONV3 ; KCAP205= HET * 4.1E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP206= HET * 7.4E7 *Exp( -TCORR * 2200 ) * CONV3 ; KCAP207= HET * 3.6E7 *Exp( -TCORR * 2200 ) * CONV3 ; KCAP208= HET * 1.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP209= HET * 5E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP210= HET * 6E7 *Exp( -TCORR * (-750) ) * CONV3 ; KCAP211= HET * 1E8 *Exp( -TCORR * (-750) ) * CONV3 ; KCAP212= HET * 6E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP213= HET * 1E6 *Exp( -TCORR * 0 ) ; KCAP214= HET * 8.6E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP215= HET * 3.0E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP216= HET * 100 *Exp( -TCORR * 0 ) ; KCAP217= HET * 3.2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP218= HET * 6.8E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP219= HET * 1.3E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP220= HET * 3.7E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP221= HET * 3.35E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP222= HET * 5.3E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP223= HET * 2.2E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP224= HET * 4.0E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP225= HET * 1.1E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP226= HET * 1.05E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP227= HET * 2.0E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP228= HET * 1.1E9 *Exp( -TCORR * 1516 ) * CONV3 ; KCAP229= HET * 1.1E6 *Exp( -TCORR * 3368 ) * CONV3 ; KCAP230= HET * 4.0E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP231= HET * 500 *Exp( -TCORR * 0 ) * CONV3 ; KCAP232= HET * 2.35E7 *Exp( -TCORR * 1395 ) * CONV3 ; KCAP233= HET * 1.38E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP234= HET * 2E9 *Exp( -TCORR * 0 ) ; KCAP235= HET * 1.1E9 *Exp( -TCORR * 1516 ) * CONV3 ; KCAP236= HET * 1.1E6 *Exp( -TCORR * 3368 ) * CONV3 ; KCAP237= HET * 4.0E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP238= HET * 500 *Exp( -TCORR * 0 ) * CONV3 ; KCAP239= HET * 2.35E7 *Exp( -TCORR * 1395 ) * CONV3 ; KCAP240= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP241= HET * 2E9 *Exp( -TCORR * 0 ) ; KCAP242= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP243= HET * 1E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP244= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP245= HET * 1E5 *Exp( -TCORR * 0 ) ; KCAP246= HET * 1.2E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP247= HET * 1E6 *Exp( -TCORR * 0 ) ; KCAP248= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP249= HET * 1.88E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP250= HET * 1.88E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP251= HET * 7.5E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP252= HET * 8.7E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP253= HET * 1.0E7 *Exp( -TCORR * 3030 ) * CONV3 ; KCAP254= HET * 4E6 *Exp( -TCORR * 3490 ) * CONV3 ; KCAP255= HET * 2E7 *Exp( -TCORR * 4090 ) * CONV3 ; KCAP256= HET * 3E5 *Exp( -TCORR * 0 ) ; KCAP257= HET * 2.1E5 *Exp( -TCORR * 0 ) ; KCAP258= HET * 1E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP259= HET * 5E4 *Exp( -TCORR * 3340 ) * CONV3 ; KCAP260= HET * 5E4 *Exp( -TCORR * 3340 ) * CONV3 ; KCAP261= HET * 4.0E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP262= HET * 1.3E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP263= HET * 6E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP264= HET * 1.7E8 *Exp( -TCORR * 400 ) * CONV3 ; KCAP265= HET * 6.2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP266= HET * 0.4 *Exp( -TCORR * 7900 ) * CONV3 ; KCAP267= HET * 23.4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP268= HET * 2.1E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP269= HET * 3.8E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP270= HET * 1.7E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP271= HET * 3.6E6 *Exp( -TCORR * 3330 ) * CONV3 ; KCAP272= HET * 6.3E6 *Exp( -TCORR * 4330 ) * CONV3 ; KCAP273= HET * 3E5 *Exp( -TCORR * 0 ) ; KCAP274= HET * 2.2E5 *Exp( -TCORR * 0 ) ; KCAP275= HET * 3.6E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP276= HET * 1.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP277= HET * 1.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP278= HET * 1.1E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP279= HET * 6.5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP280= HET * 1.7E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP281= HET * 5.0E7 *Exp( -TCORR * 780 ) * CONV3 ; KCAP282= HET * 3.3E7 *Exp( -TCORR * 650 ) * CONV3 ; KCAP283= HET * 0.031 *Exp( -TCORR * 7500 ) * CONV3 ; KCAP284= HET * 4.2E6 *Exp( -TCORR * 0 ) ; KCAP285= HET * 1.7E7 *Exp( -TCORR * 1900 ) * CONV3 ; KCAP286= HET * 3.9E8 *Exp( -TCORR * 2840 ) * CONV3 ; KCAP287= HET * 4.1E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP288= HET * 2.8E6 *Exp( -TCORR * 2090 ) * CONV3 ; KCAP289= HET * 4.1E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP290= HET * 4.1E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP291= HET * 2.7E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP292= HET * 2.7E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP293= HET * 1.1E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP294= HET * 2.2E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP295= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP296= HET * 1.5E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP297= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP298= HET * 4.3E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP299= HET * 4.3E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP300= HET * 6.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP301= HET * 6.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP302= HET * 1E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP303= HET * 5.0E6 *Exp( -TCORR * 470 ) * CONV3 ; KCAP304= HET * 1E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP305= HET * 1E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP306= HET * 5.0E3* Exp( -TCORR * 0 ) * CONV3; KCAP307= HET * 1.0E4* Exp( -TCORR * 0 ) * CONV3; KCAP308= HET * 5.0E5* Exp( -TCORR * 0 ) * CONV3; KCAP309= HET * 2.0E9* Exp( -TCORR * 0 ) * CONV3; KCAP310= HET * 1.0E8* Exp( -TCORR * 0 ) * CONV3; KCAP311= HET * 1.7E4* Exp( -TCORR * 0 ) ; KCAP312= HET * 3.0E4 * Exp( -TCORR * 0 ) * CONV3; KCAP313= HET * 0.0126 * Exp( -TCORR * 0 ) * CONV3; E1FWD= HET * 2.34E-5 * exp(-TCORR* 6800 ); E2FWD= HET * 4.3E-2 * exp(-TCORR* 9250 ) * CONV3 ; E3FWD= HET * 1E7 * exp(-TCORR* 0 ) ; E4FWD= HET * 2.35 * exp(-TCORR* 1820 ) ; E5FWD= HET * 5E11 * exp(-TCORR* (-6890) ) ; E6FWD= HET * 6.02E5 * exp(-TCORR* 560 ) * CONV3 ; E7FWD= HET * 8.0E5 * exp(-TCORR* 0 ) ; E8FWD= HET * 1.1E12 * exp(-TCORR* (-1800) ) ; E9FWD= HET * 2.65E7 * exp(-TCORR* 1760 ) ; E10FWD= HET * 5E5 * exp(-TCORR* 0 ) ; E11FWD= HET * 1.0E7 * exp(-TCORR* 0 ) * CONV3 ; E12FWD= HET * 6.27E4 * exp(-TCORR* (-1940) ) * CONV3 ; E13FWD= HET * 3110 * exp(-TCORR* (-1960) ) ; E14FWD= HET * 5E13 * exp(-TCORR* 0 ) ; E15FWD= HET * 1.02E9 * exp(-TCORR* (-2700) ) ; E16FWD= HET * 8.85E6 * exp(-TCORR* (-12) ) ; E17FWD= HET * 8.75E5 * exp(-TCORR* (-46) ) ; E18FWD= HET * 4.7E4 * exp(-TCORR* 0 ) * CONV3 ; E19FWD= HET * 1.1E3 * exp(-TCORR* 0 ) * CONV3 ; E20FWD= HET * 3.2E3 * exp(-TCORR* 0 ) * CONV3 ; E21FWD= HET * 0.18 * exp(-TCORR* (-4030) ) * CONV3 ; E22FWD= HET * 1.4E-4 * exp(-TCORR* (-2500) ) * CONV3 ; E23FWD= HET * 790 * exp(-TCORR* 2500 ) * CONV3 ; E24FWD= HET * 2.5E7 * exp(-TCORR* 2450 ) * CONV3 ; E25FWD= HET * 8.5E9 * exp(-TCORR* 0 ) * CONV3; E26FWD= HET * 1.2E10 * exp(-TCORR* 0 ) * CONV3; E27FWD= HET * 4.3E9 * exp(-TCORR* 0 ) * CONV3; E28FWD= HET * 2.1E10 * exp(-TCORR* 0 ) * CONV3; E29FWD= HET * 1.0E4 * exp(-TCORR* 0 ) * CONV3; E30FWD= HET * 1.1E10 * exp(-TCORR* 0 ) * CONV3; E31FWD= HET * 4.4E10 * exp(-TCORR* 0 ) * CONV3; E32FWD= HET * 1.9E8 * exp(-TCORR* 0 ) * CONV3; E33FWD= HET * 1.86E10 * exp(-TCORR* 0 ) * CONV3 ; E34FWD= HET * 1.5E8 * exp(-TCORR* 0 ) * CONV3 ; E35FWD= HET * 1.1E6 * exp(-TCORR* 0 ) * CONV3 ; E36FWD= HET * 1E7 * exp(-TCORR* 0 ) * CONV3 ; E37FWD= HET * 2E5 * exp(-TCORR* 0 ) * CONV3 ; E38FWD= HET * 3.1E7 * exp(-TCORR* 0 ) * CONV3 ; E39FWD= HET * 3.5E8 * exp(-TCORR* 0 ) * CONV3 ; E40FWD= HET * 330 * exp(-TCORR* 0 ) ; E41FWD= HET * 5E4 * exp(-TCORR* 0 ) ; E42FWD= HET * 5.9E4 * exp(-TCORR* 0 ) ; E43FWD= HET * 7.5E5 * exp(-TCORR* 0 ) ; E44FWD= HET * 3.2E9 * exp(-TCORR* 0 ) ; E45FWD= HET * 2.6E6 * exp(-TCORR* 0 ) ; E46FWD= HET * 6.32E6 * exp(-TCORR* 0 ) ; E47FWD= HET * 21.5 * exp(-TCORR* 0 ) * CONV3 ; E48FWD= HET * 3E-3 * exp(-TCORR* 0 ) ; E49FWD= HET * 3E-3 * exp(-TCORR* 0 ) ; E50FWD= HET * 3E-3 * exp(-TCORR* 0 ) ; E51FWD= HET * 2.52E8 * exp(-TCORR* 0 ) * CONV3 ; E52FWD= HET * 3.4E8 * exp(-TCORR* 4300 ) * CONV3 ; E53FWD= HET * 3E3 * exp(-TCORR* 0 ) * CONV3 ; E54FWD= HET * 315 * exp(-TCORR* 0 ) ; E55FWD= HET * 1.1E7 * exp(-TCORR* 0 ) * CONV3 ; E56FWD= HET * 3.4E4 * exp(-TCORR* 0 ) * CONV3 * CONV3 ; E57FWD= HET * 8.75E5 * exp(-TCORR* (-46) ) ; E58FWD= HET * 8.75E5 * exp(-TCORR* (-46) ) ; E1BCK= HET * 1.3E11 * exp(-TCORR* 0 ) * CONV3 ; E2BCK= HET * 5.6E4 * exp(-TCORR* 8500 ) ; E3BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E4BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E5BCK= HET * 2.9E5 * exp(-TCORR* 0 ) * CONV3 ; E6BCK= HET * 3.4E10 * exp(-TCORR* 0 ) * CONV3 ; E7BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E8BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E9BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E10BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E11BCK= HET * 4.6E-3 * exp(-TCORR* 0 ) ; E12BCK= HET * 2.0E8 * exp(-TCORR* 0 ) * CONV3 ; E13BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E14BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E15BCK= HET * 1E11 * exp(-TCORR* 0 ) * CONV3 ; E16BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E17BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E18BCK= HET * 4.3E8 * exp(-TCORR* 0 ) * CONV3 ; E19BCK= HET * 8.0E9 * exp(-TCORR* 0 ) * CONV3 ; E20BCK= HET * 1.8E5 * exp(-TCORR* 0 ) ; E21BCK= HET * 5.1E-3 * exp(-TCORR* 0 ) ; E22BCK= HET * 5.69E-3 * exp(-TCORR* 0 ) ; E23BCK= HET * 3.95E-6 * exp(-TCORR* 0 ) *CONV3 ; E24BCK= HET * 3.95E-6 * exp(-TCORR* 5530 ) * CONV3 ; E25BCK= HET * 6.0E4 * exp(-TCORR* 0 ) ; E26BCK= HET * 1.9E4 * exp(-TCORR* 0 ) ; E27BCK= HET * 6.1E9 * exp(-TCORR* 0 ) ; E28BCK= HET * 4100 * exp(-TCORR* 0 ) * CONV3 ; E29BCK= HET * 4.5E7 * exp(-TCORR* 0 ) * CONV3 ; E30BCK= HET * 3.3E7 * exp(-TCORR* 0 ) ; E31BCK= HET * 2.45E-2 * exp(-TCORR* 0 ) * CONV3 ; E32BCK= HET * 2.7E6 * exp(-TCORR* 0 ) * CONV3 ; E33BCK= HET * 2.0E10 * exp(-TCORR* 0 ) * CONV3 ; E34BCK= HET * 6.5E3 * exp(-TCORR* 0 ) ; E35BCK= HET * 6.5E6 * exp(-TCORR* 0 ) * CONV3 ; E36BCK= HET * 1E10 * exp(-TCORR* 0 ) * CONV3 ; E37BCK= HET * 2E10 * exp(-TCORR* 0 ) * CONV3 ; E38BCK= HET * 1.033E3 * exp(-TCORR* 0 ) ; E39BCK= HET * 3E4 * exp(-TCORR* 0 ) ; E40BCK= HET * 5.2E10 * exp(-TCORR* 0 ) * CONV3 ; E41BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E42BCK= HET * 4.4E10 * exp(-TCORR* 0 ) * CONV3 ; E43BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E44BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E45BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E46BCK= HET * 2E10 * exp(-TCORR* 0 ) * CONV3 ; E47BCK= HET * 5.5E-3 * exp(-TCORR* 0 ) ; E48BCK= HET * 7.5E6 * exp(-TCORR* 0 ) * CONV3; E49BCK= HET * 1.89E4 * exp(-TCORR* 0 ) * CONV3; E50BCK= HET * 114 * exp(-TCORR* 0 ) * CONV3; E51BCK= HET * 2.1E8 * exp(-TCORR* 0 ) * CONV3; E52BCK= HET * 1E8 * exp(-TCORR* 0 ) * CONV3; E53BCK= HET * 2160 * exp(-TCORR* 0 ) ; E54BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3; E55BCK= HET * 3E4 * exp(-TCORR* 0 ) ; E56BCK= HET * 1.1E5 * exp(-TCORR* 0 ) * CONV3; E57BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3; E58BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3; * **********************************************************; *O2a = 3E10-4/1000; OHm = 3.16E-10/1000; Hp = 3.16E-5/1000; Clm = 1E-4/1000; Brm = 3E-6/1000; Feppp = 5E-6/1000; Mnppp = 2.5E-7/1000; Cupp = 2.5E-7/1000; HSO4m = 3E-7/1000; SO4mm = 5.97E-5/1000; C2O4mm = 5E-6/1000; H2Oa=0.0556; CH3COOm = 0.0 ; ACO3a = 0.0 ; Br2a = 0.0 ; Br2m = 0.0 ; Bra = 0.0 ; BrCla = 0.0 ; BrNO2a = 0.0 ; BrOHm = 0.0 ; C2O4m = 0.0 ; CH2COOHa = 0.0 ; CH2OH2a = 0.0 ; CH2OHa = 0.0 ; HOA3m = 0.0 ; CH3a = 0.0 ; CH3CH2Oa = 0.0 ; CH3CHOa = 0.0 ; ALCOOL1a = 0.0 ; CH3CHOHa = 0.0 ; CH3COa = 0.0 ; CH3COH2a = 0.0 ; HOP1a = 0.0 ; CH3COOHa = 0.0 ; CH3COOOm = 0.0 ; CH3O2a = 0.0 ; CH3O2Ha = 0.0 ; CH3Oa = 0.0 ; CH3OHa = 0.0 ; CHOH2a = 0.0 ; CHOH22a = 0.0 ; HOA2a = 0.0 ; HOA2m = 0.0 ; OHCH2SO3a = 0.0 ; CHOHSO3m = 0.0 ; CHOSO3m = 0.0 ; CHOSO3mm = 0.0 ; Cl2a = 0.0 ; Cl2a = 0.0 ; Cl2m = 0.0 ; Cla = 0.0 ; ClNO2a = 0.0 ; ClOHm = 0.0 ; CO2a = 0.0 ; CO2Ha = 0.0 ; CO2m = 0.0 ; CO3m = 0.0 ; CO3mm = 0.0 ; COH22a = 0.0 ; HOA1a = 0.0 ; COOHa = 0.0 ; CuOHpp = 0.0 ; Cup = 0.0 ; ETHPa = 0.0 ; ETOHa = 0.0 ; FeC2O42m = 0.0 ; FeORG3 = 0.0 ; FeC2O4p = 0.0 ; FeMO2pp = 0.0 ; FeClpp = 0.0 ; HFe4p = 0.0 ; FeOH2p = 0.0 ; FeOHpp = 0.0 ; FeOpp = 0.0 ; Fepp = 0.0 ; FeSO4p = 0.0 ; GLYa = 0.0 ; H2C2O4a = 0.0 ; H2CO3a = 0.0 ; H2O2a = 0.0 ; H2SO4a = 0.0 ; CH3COOHa = 0.0 ; HC2O4m = 0.0 ; HCHOa = 0.0 ; HCla = 0.0 ; HCO3m = 0.0 ; HCOOHa = 0.0 ; HCOOm = 0.0 ; HMSm = 0.0 ; HNO2a = 0.0 ; HNO3a = 0.0 ; HNO4a = 0.0 ; HO2a = 0.0 ; HO3a = 0.0 ; HOBra = 0.0 ; HOCla = 0.0 ; HSO3m = 0.0 ; HSO5m = 0.0 ; MnBr2p = 0.0 ; MnCl2p = 0.0 ; MnHSO3p = 0.0 ; MnO2p = 0.0 ; MnOH2p = 0.0 ; MnOHpp = 0.0 ; MnOpp = 0.0 ; Mnpp = 0.0 ; Mnpppp = 0.0 ; N2O5a = 0.0 ; NH3a = 0.0 ; NH4p = 0.0 ; NO2a = 0.0 ; NO2m = 0.0 ; NO2p = 0.0 ; NO3a = 0.0 ; NO3m = 0.0 ; NOa = 0.0 ; HOA4m = 0.0 ; O2CH2OHa = 0.0 ; HOP3a = 0.0 ; O2CHOa = 0.0 ; HS2m = 0.0 ; HOP2a = 0.0 ; HOA3a = 0.0 ; O2m = 0.0 ; O2NO2m = 0.0 ; O3a = 0.0 ; O3m = 0.0 ; O3Pa = 0.0 ; OHa = 0.0 ; PAAa = 0.0 ; PAAa = 0.0 ; S2O8mm = 0.0 ; SO2a = 0.0 ; SO3a = 0.0 ; SO3m = 0.0 ; SO3mm = 0.0 ; SO4m = 0.0 ; SO5m = 0.0 ; SO5O2Hm = 0.0 ; SO5O2mm = 0.0 ; ** ; COMPILE GAS; HET=0; KMTCO2 =HET/(((RAD**2)/(3* DGCO2 ))+(4*RAD/(3* VCO2 *ACCO2 ))); KMTHCl =HET/(((RAD**2)/(3* DGHCl ))+(4*RAD/(3* VHCl *ACHCl ))); KMTNH3 =HET/(((RAD**2)/(3* DGNH3 ))+(4*RAD/(3* VNH3 *ACNH3 ))); KMTO3 =HET/(((RAD**2)/(3* DGO3 ))+(4*RAD/(3* VO3 *ACO3 ))); KMTHO2 =HET/(((RAD**2)/(3* DGHO2 ))+(4*RAD/(3* VHO2 *ACHO2 ))); KMTOH =HET/(((RAD**2)/(3* DGOH ))+(4*RAD/(3* VOH *ACOH ))); KMTH2O2 =HET/(((RAD**2)/(3* DGH2O2 ))+(4*RAD/(3* VH2O2 *ACH2O2 ))); KMTHNO3 =HET/(((RAD**2)/(3* DGHNO3 ))+(4*RAD/(3* VHNO3 *ACHNO3 ))); KMTNO3 =HET/(((RAD**2)/(3* DGNO3 ))+(4*RAD/(3* VNO3 *ACNO3 ))); KMTN2O5 =HET/(((RAD**2)/(3* DGN2O5 ))+(4*RAD/(3* VN2O5 *ACN2O5 ))); KMTNO2 =HET/(((RAD**2)/(3* DGNO2 ))+(4*RAD/(3* VNO2 *ACNO2 ))); KMTHNO2 =HET/(((RAD**2)/(3* DGHNO2 ))+(4*RAD/(3* VHNO2 *ACHNO2 ))); KMTHNO4 =HET/(((RAD**2)/(3* DGHNO4 ))+(4*RAD/(3* VHNO4 *ACHNO4 ))); KMTSO2 =HET/(((RAD**2)/(3* DGSO2 ))+(4*RAD/(3* VSO2 *ACSO2 ))); KMTHCHO =HET/(((RAD**2)/(3* DGHCHO ))+(4*RAD/(3* VHCHO *ACHCHO ))); KMTOP1 =HET/(((RAD**2)/(3* DGOP1 ))+(4*RAD/(3* VOP1 *ACOP1 ))); KMTPAA =HET/(((RAD**2)/(3* DGPAA ))+(4*RAD/(3* VPAA *ACPAA ))); KMTCH3OH =HET/(((RAD**2)/(3*DGCH3OH))+(4*RAD/(3*VCH3OH*ACCH3OH))); KMTETOH =HET/(((RAD**2)/(3* DGETOH ))+(4*RAD/(3* VETOH *ACETOH ))); KMTALD =HET/(((RAD**2)/(3* DGALD ))+(4*RAD/(3* VALD *ACALD ))); KMTORA1 =HET/(((RAD**2)/(3* DGORA1 ))+(4*RAD/(3* VORA1 *ACORA1 ))); KMTORA2 =HET/(((RAD**2)/(3* DGORA2 ))+(4*RAD/(3* VORA2 *ACORA2 ))); KMTMO2 =HET/(((RAD**2)/(3* DGMO2 ))+(4*RAD/(3* VMO2 *ACMO2 ))); KMTETHP =HET/(((RAD**2)/(3* DGETHP ))+(4*RAD/(3* VETHP *ACETHP ))); KMTCl2 =HET/(((RAD**2)/(3* DGCl2 ))+(4*RAD/(3* VCl2 *ACCl2 ))); KMTBr2 =HET/(((RAD**2)/(3* DGBr2 ))+(4*RAD/(3* VBr2 *ACBr2 ))); KMTSULF =HET/(((RAD**2)/(3*DGSULF ))+(4*RAD/(3* VSULF *ACSULF ))); KMTACO3 =HET/(((RAD**2)/(3* DGACO3 ))+(4*RAD/(3* VACO3 *ACACO3 ))); KMTGLY =HET/(((RAD**2)/(3* DGGLY ))+(4*RAD/(3* VGLY *ACGLY ))); KMTO2 =HET/(((RAD**2)/(3* DGO2 ))+(4*RAD/(3* VO2 *ACO2 ))); KMTClNO2 =HET/(((RAD**2)/(3*DGClNO2))+(4*RAD/(3*VClNO2*ACClNO2))); KMTBrNO2 =HET/(((RAD**2)/(3*DGBrNO2))+(4*RAD/(3*VBrNO2*ACBrNO2))); KMTBrCl =HET/(((RAD**2)/(3* DGBrCl ))+(4*RAD/(3* VBrCl *ACBrCl ))); KMTNO =HET/(((RAD**2)/(3* DGNO ))+(4*RAD/(3* VNO *ACNO ))); * *********; * ; * Kinetics in aqueous solution Capram 2.4; * ; * ; KCAP1= HET * 50 *Exp( -TCORR * 0 ) * CONV3 ; KCAP2= HET * 7.3E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP3= HET * 2.8E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP4= HET * 1E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP5= HET * 7.0E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP6= HET * 1.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP7= HET * 1.3E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP8= HET * 1.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP9= HET * 1.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP10= HET * 1.0E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP11= HET * 1.2E6 *Exp( -TCORR * 5050 ) * CONV3 ; KCAP12= HET * 4.3E8 *Exp( -TCORR * 1100 ) * CONV3 ; KCAP13= HET * 1E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP14= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP15= HET * 1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP16= HET * 2.3E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP17= HET * 3E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP18= HET * 1E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP19= HET * 8E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP20= HET * 1.3E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP21= HET * 1.3E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP22= HET * 1.3E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP23= HET * 1.6E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP24= HET * 1.5E9 *Exp( -TCORR * 2200 ) * CONV3 ; KCAP25= HET * 330 *Exp( -TCORR * 4500 ) ; KCAP26= HET * 4.15E5 *Exp( -TCORR * 2720 ) * CONV3 ; KCAP27= HET * 4.85E7 *Exp( -TCORR * 1060 ) * CONV3 ; KCAP28= HET * 1.0E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP29= HET * 1.1E10 *Exp( -TCORR * 2120 ) * CONV3 ; KCAP30= HET * 3.0E7 *Exp( -TCORR * 1680 ) * CONV3 ; KCAP31= HET * 2.4E7 *Exp( -TCORR * 1680 ) * CONV3 ; KCAP32= HET * 6.0E6 *Exp( -TCORR * 1680 ) * CONV3 ; KCAP33= HET * 2.7E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP34= HET * 4.6E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP35= HET * 6E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP36= HET * 4.6E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP37= HET * 8.2E5 *Exp( -TCORR * 4690.0 ) * CONV3 ; KCAP38= HET * 9.5E3 *Exp( -TCORR * 2766.0 ) * CONV3 ; KCAP39= HET * 2.34E-4 *Exp( -TCORR * 4089.0 ) * CONV3; KCAP40= HET * 2.0E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP41= HET * 1.0E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP42= HET * 1.1E4 *Exp( -TCORR * 4150 ) * CONV3 ; KCAP43= HET * 2.5E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP44= HET * 100 *Exp( -TCORR * 0 ) * CONV3 ; KCAP45= HET * 1.0E4 *Exp( -TCORR * 2562 ) * CONV3 ; KCAP46= HET * 160 *Exp( -TCORR * 2680 ) * CONV3 ; KCAP47= HET * 3E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP48= HET * 2.5E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP49= HET * 400 *Exp( -TCORR * 5352 ) * CONV3 ; KCAP50= HET * 1E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP51= HET * 7.2E4 *Exp( -TCORR * 842.0 ) * CONV3 ; KCAP52= HET * 1.8E4 *Exp( -TCORR * 5052.0 ) * CONV3 ; ; KCAP53= HET * 1.95 *Exp( -TCORR * 5653 ) * CONV3 ; KCAP54= HET * 0.49 *Exp( -TCORR * 8780.0 ) ; KCAP55= HET * 1E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP56= HET * 1.65E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP57= HET * 3E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP58= HET * 5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP59= HET * 1.0E9 *Exp( -TCORR * 0.0 ) ; KCAP60= HET * 8.9E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP61= HET * 9.4E7 *Exp( -TCORR * 2700 ) * CONV3 ; KCAP62= HET * 8.0E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP63= HET * 1.1E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP64= HET * 4.9E6 *Exp( -TCORR * 2000 ) * CONV3 ; KCAP65= HET * 4.9E6 *Exp( -TCORR * 2000 ) * CONV3 ; KCAP66= HET * 3.0E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP67= HET * 3E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP68= HET * 1.3E9 *Exp( -TCORR * 2000 ) * CONV3 ; KCAP69= HET * 3.0E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP70= HET * 2.6E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP71= HET * 1.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP72= HET * 1.2E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP73= HET * 1E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP74= HET * 8.4E7 *Exp( -TCORR * (-2900) ) * CONV3 ; KCAP75= HET * 4.5E-2 *Exp( -TCORR * 0 ) ; KCAP76= HET * 9.1E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP77= HET * 7.2E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP78= HET * 1.4E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP79= HET * 6E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP80= HET * 1.2E7 *Exp( -TCORR * 1720 ) * CONV3 ; KCAP81= HET * 6.6E5 *Exp( -TCORR * 850 ) * CONV3 ; KCAP82= HET * 5E5 *Exp( -TCORR * 6900 ) * CONV3 ; KCAP83= HET * 1.1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP84= HET * 3.3E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP85= HET * 1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP86= HET * 1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP87= HET * 5E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP88= HET * 2.55E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP89= HET * 10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP90= HET * 3E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP91= HET * 1.3E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP92= HET * 4.2E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP93= HET * 5.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP94= HET * 5E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP95= HET * 1E5 *Exp( -TCORR * 0 ) ; KCAP96= HET * 1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP97= HET * 2.6E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP98= HET * 1.7E4 *Exp( -TCORR * 0 ) ; KCAP99= HET * 7.0E3 *Exp( -TCORR * 0 ) ; KCAP100= HET * 44.32 *Exp( -TCORR * 0 ) * CONV3 ; KCAP101= HET * 1.6E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP102= HET * 7.2E7 *Exp( -TCORR * 4000 ) * CONV3 * CONV3 ; KCAP103= HET * 1.7E7 *Exp( -TCORR * 3800 ) * CONV3 * CONV3 ; KCAP104= HET * 5.6E7 *Exp( -TCORR * 3990 ) * CONV3 * CONV3 ; KCAP105= HET * 2.4E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP106= HET * 3.7E5 *Exp( -TCORR * 5530 ) * CONV3 ; KCAP107= HET * 1.5E9 *Exp( -TCORR * 5280 ) * CONV3 ; KCAP108= HET * 30 *Exp( -TCORR * 0 ) * CONV3 ; KCAP109= HET * 2.65E7 *Exp( -TCORR * 5809 ) * CONV3 ; KCAP110= HET * 3E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP111= HET * 1E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP112= HET * 4.6E9 *Exp( -TCORR * (-2165) ) * CONV3 ; KCAP113= HET * 17 *Exp( -TCORR * 0 ) * CONV3 ; KCAP114= HET * 4.8E7 *Exp( -TCORR * 2600 ) * CONV3 ; KCAP115= HET * 2.2E8 *Exp( -TCORR * 2600 ) * CONV3 ; KCAP116= HET * 1.7E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP117= HET * 2.34E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP118= HET * 3.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP119= HET * 5.2E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP120= HET * 1200 *Exp( -TCORR * 0 ) ; KCAP121= HET * 2.5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP122= HET * 8.6E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP123= HET * 3.6E2 *Exp( -TCORR * 0 ) * CONV3 ; KCAP124= HET * 2.13E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP125= HET * 5.5E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP126= HET * 3.5E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP127= HET * 6.1E8 *Exp( -TCORR * 840 ) * CONV3 ; KCAP128= HET * 5.8E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP129= HET * 3.4E8 *Exp( -TCORR * 1200 ) * CONV3 ; KCAP130= HET * 1.4E7 *Exp( -TCORR * 4089 ) * CONV3 ; KCAP131= HET * 3E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP132= HET * 1.7E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP133= HET * 1.7E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP134= HET * 3.5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP135= HET * 3.5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP136= HET * 5.0E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP137= HET * 1.4E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP138= HET * 11 *Exp( -TCORR * 1110 ) * CONV3 ; KCAP139= HET * 7.14E6 *Exp( -TCORR * 0 ) * CONV3 * CONV3 ; KCAP140= HET * 7.14E6 *Exp( -TCORR * 0 ) * CONV3 * CONV3 ; KCAP141= HET * 1.7E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP142= HET * 1.0E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP143= HET * 1.3E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP144= HET * 1.0E9 *Exp( -TCORR * 580 ) * CONV3 ; KCAP145= HET * 9.0E6 *Exp( -TCORR * 2190 ) * CONV3 ; KCAP146= HET * 5.4E5 *Exp( -TCORR * 4300 ) * CONV3 ; KCAP147= HET * 5.1E4 *Exp( -TCORR * 5500 ) * CONV3 ; KCAP148= HET * 1000 *Exp( -TCORR * 0 ) * CONV3 ; KCAP149= HET * 2.6E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP150= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP151= HET * 1.65E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP152= HET * 1.05E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP153= HET * 1.9E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP154= HET * 4.1E7 *Exp( -TCORR * 1760 ) * CONV3 ; KCAP155= HET * 2.2E6 *Exp( -TCORR * 3300 ) * CONV3 ; KCAP156= HET * 1.2E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP157= HET * 3.8E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP158= HET * 1.5E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP159= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP160= HET * 52 *Exp( -TCORR * 7217 ) ; KCAP161= HET * 8E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP162= HET * 1.0E9 *Exp( -TCORR * 1020 ) * CONV3 ; KCAP163= HET * 1.4E7 *Exp( -TCORR * 1300 ) * CONV3 ; KCAP164= HET * 1.0E6 *Exp( -TCORR * 4500 ) * CONV3 ; KCAP165= HET * 3.6E4 *Exp( -TCORR * 4330 ) * CONV3 ; KCAP166= HET * 3E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP167= HET * 1.3E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP168= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP169= HET * 1.2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP170= HET * 3.6E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP171= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP172= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP173= HET * 1.9E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP174= HET * 1.9E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP175= HET * 4E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP176= HET * 4E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP177= HET * 4E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP178= HET * 4E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP179= HET * 1E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP180= HET * 1E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP181= HET * 1.3E8 *Exp( -TCORR * 1000 ) * CONV3 ; KCAP182= HET * 3.2E9 *Exp( -TCORR * 1000 ) * CONV3 ; KCAP183= HET * 2.5E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP184= HET * 2.1E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP185= HET * 3.8E5 *Exp( -TCORR * 3400 ) * CONV3 ; KCAP186= HET * 5.1E7 *Exp( -TCORR * 2200 ) * CONV3 ; KCAP187= HET * 8E4 *Exp( -TCORR * 4450 ) * CONV3 ; KCAP188= HET * 1.3E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP189= HET * 4E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP190= HET * 4.9E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP191= HET * 1.4E5 *Exp( -TCORR * 3300 ) * CONV3 ; KCAP192= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP193= HET * 1.5E7 *Exp( -TCORR * 1330 ) * CONV3 ; KCAP194= HET * 1.0E8 *Exp( -TCORR * 1800 ) * CONV3 ; KCAP195= HET * 2.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP196= HET * 2.8E7 *Exp( -TCORR * 1210 ) * CONV3 ; KCAP197= HET * 1.4E4 *Exp( -TCORR * 3800 ) * CONV3 ; KCAP198= HET * 2.9E6 *Exp( -TCORR * 3800 ) * CONV3 ; KCAP199= HET * 1500 *Exp( -TCORR * 4930 ) * CONV3 ; KCAP200= HET * 2.6E5 *Exp( -TCORR * 4800 ) * CONV3 ; KCAP201= HET * 10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP202= HET * 1.7E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP203= HET * 100 *Exp( -TCORR * 0 ) * CONV3 ; KCAP204= HET * 580 *Exp( -TCORR * 0 ) * CONV3 ; KCAP205= HET * 4.1E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP206= HET * 7.4E7 *Exp( -TCORR * 2200 ) * CONV3 ; KCAP207= HET * 3.6E7 *Exp( -TCORR * 2200 ) * CONV3 ; KCAP208= HET * 1.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP209= HET * 5E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP210= HET * 6E7 *Exp( -TCORR * (-750) ) * CONV3 ; KCAP211= HET * 1E8 *Exp( -TCORR * (-750) ) * CONV3 ; KCAP212= HET * 6E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP213= HET * 1E6 *Exp( -TCORR * 0 ) ; KCAP214= HET * 8.6E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP215= HET * 3.0E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP216= HET * 100 *Exp( -TCORR * 0 ) ; KCAP217= HET * 3.2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP218= HET * 6.8E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP219= HET * 1.3E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP220= HET * 3.7E3 *Exp( -TCORR * 0 ) * CONV3 ; KCAP221= HET * 3.35E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP222= HET * 5.3E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP223= HET * 2.2E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP224= HET * 4.0E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP225= HET * 1.1E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP226= HET * 1.05E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP227= HET * 2.0E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP228= HET * 1.1E9 *Exp( -TCORR * 1516 ) * CONV3 ; KCAP229= HET * 1.1E6 *Exp( -TCORR * 3368 ) * CONV3 ; KCAP230= HET * 4.0E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP231= HET * 500 *Exp( -TCORR * 0 ) * CONV3 ; KCAP232= HET * 2.35E7 *Exp( -TCORR * 1395 ) * CONV3 ; KCAP233= HET * 1.38E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP234= HET * 2E9 *Exp( -TCORR * 0 ) ; KCAP235= HET * 1.1E9 *Exp( -TCORR * 1516 ) * CONV3 ; KCAP236= HET * 1.1E6 *Exp( -TCORR * 3368 ) * CONV3 ; KCAP237= HET * 4.0E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP238= HET * 500 *Exp( -TCORR * 0 ) * CONV3 ; KCAP239= HET * 2.35E7 *Exp( -TCORR * 1395 ) * CONV3 ; KCAP240= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP241= HET * 2E9 *Exp( -TCORR * 0 ) ; KCAP242= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP243= HET * 1E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP244= HET * 2E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP245= HET * 1E5 *Exp( -TCORR * 0 ) ; KCAP246= HET * 1.2E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP247= HET * 1E6 *Exp( -TCORR * 0 ) ; KCAP248= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP249= HET * 1.88E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP250= HET * 1.88E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP251= HET * 7.5E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP252= HET * 8.7E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP253= HET * 1.0E7 *Exp( -TCORR * 3030 ) * CONV3 ; KCAP254= HET * 4E6 *Exp( -TCORR * 3490 ) * CONV3 ; KCAP255= HET * 2E7 *Exp( -TCORR * 4090 ) * CONV3 ; KCAP256= HET * 3E5 *Exp( -TCORR * 0 ) ; KCAP257= HET * 2.1E5 *Exp( -TCORR * 0 ) ; KCAP258= HET * 1E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP259= HET * 5E4 *Exp( -TCORR * 3340 ) * CONV3 ; KCAP260= HET * 5E4 *Exp( -TCORR * 3340 ) * CONV3 ; KCAP261= HET * 4.0E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP262= HET * 1.3E10 *Exp( -TCORR * 0 ) * CONV3 ; KCAP263= HET * 6E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP264= HET * 1.7E8 *Exp( -TCORR * 400 ) * CONV3 ; KCAP265= HET * 6.2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP266= HET * 0.4 *Exp( -TCORR * 7900 ) * CONV3 ; KCAP267= HET * 23.4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP268= HET * 2.1E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP269= HET * 3.8E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP270= HET * 1.7E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP271= HET * 3.6E6 *Exp( -TCORR * 3330 ) * CONV3 ; KCAP272= HET * 6.3E6 *Exp( -TCORR * 4330 ) * CONV3 ; KCAP273= HET * 3E5 *Exp( -TCORR * 0 ) ; KCAP274= HET * 2.2E5 *Exp( -TCORR * 0 ) ; KCAP275= HET * 3.6E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP276= HET * 1.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP277= HET * 1.0E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP278= HET * 1.1E4 *Exp( -TCORR * 0 ) * CONV3 ; KCAP279= HET * 6.5E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP280= HET * 1.7E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP281= HET * 5.0E7 *Exp( -TCORR * 780 ) * CONV3 ; KCAP282= HET * 3.3E7 *Exp( -TCORR * 650 ) * CONV3 ; KCAP283= HET * 0.031 *Exp( -TCORR * 7500 ) * CONV3 ; KCAP284= HET * 4.2E6 *Exp( -TCORR * 0 ) ; KCAP285= HET * 1.7E7 *Exp( -TCORR * 1900 ) * CONV3 ; KCAP286= HET * 3.9E8 *Exp( -TCORR * 2840 ) * CONV3 ; KCAP287= HET * 4.1E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP288= HET * 2.8E6 *Exp( -TCORR * 2090 ) * CONV3 ; KCAP289= HET * 4.1E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP290= HET * 4.1E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP291= HET * 2.7E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP292= HET * 2.7E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP293= HET * 1.1E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP294= HET * 2.2E6 *Exp( -TCORR * 0 ) * CONV3 ; KCAP295= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP296= HET * 1.5E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP297= HET * 2E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP298= HET * 4.3E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP299= HET * 4.3E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP300= HET * 6.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP301= HET * 6.5E8 *Exp( -TCORR * 0 ) * CONV3 ; KCAP302= HET * 1E7 *Exp( -TCORR * 0 ) * CONV3 ; KCAP303= HET * 5.0E6 *Exp( -TCORR * 470 ) * CONV3 ; KCAP304= HET * 1E9 *Exp( -TCORR * 0 ) * CONV3 ; KCAP305= HET * 1E5 *Exp( -TCORR * 0 ) * CONV3 ; KCAP306= HET * 5.0E3* Exp( -TCORR * 0 ) * CONV3; KCAP307= HET * 1.0E4* Exp( -TCORR * 0 ) * CONV3; KCAP308= HET * 5.0E5* Exp( -TCORR * 0 ) * CONV3; KCAP309= HET * 2.0E9* Exp( -TCORR * 0 ) * CONV3; KCAP310= HET * 1.0E8* Exp( -TCORR * 0 ) * CONV3; KCAP311= HET * 1.7E4* Exp( -TCORR * 0 ) ; KCAP312= HET * 3.0E4 * Exp( -TCORR * 0 ) * CONV3; KCAP313= HET * 0.0126 * Exp( -TCORR * 0 ) * CONV3; E1FWD= HET * 2.34E-5 * exp(-TCORR* 6800 ); E2FWD= HET * 4.3E-2 * exp(-TCORR* 9250 ) * CONV3 ; E3FWD= HET * 1E7 * exp(-TCORR* 0 ) ; E4FWD= HET * 2.35 * exp(-TCORR* 1820 ) ; E5FWD= HET * 5E11 * exp(-TCORR* (-6890) ) ; E6FWD= HET * 6.02E5 * exp(-TCORR* 560 ) * CONV3 ; E7FWD= HET * 8.0E5 * exp(-TCORR* 0 ) ; E8FWD= HET * 1.1E12 * exp(-TCORR* (-1800) ) ; E9FWD= HET * 2.65E7 * exp(-TCORR* 1760 ) ; E10FWD= HET * 5E5 * exp(-TCORR* 0 ) ; E11FWD= HET * 1.0E7 * exp(-TCORR* 0 ) * CONV3 ; E12FWD= HET * 6.27E4 * exp(-TCORR* (-1940) ) * CONV3 ; E13FWD= HET * 3110 * exp(-TCORR* (-1960) ) ; E14FWD= HET * 5E13 * exp(-TCORR* 0 ) ; E15FWD= HET * 1.02E9 * exp(-TCORR* (-2700) ) ; E16FWD= HET * 8.85E6 * exp(-TCORR* (-12) ) ; E17FWD= HET * 8.75E5 * exp(-TCORR* (-46) ) ; E18FWD= HET * 4.7E4 * exp(-TCORR* 0 ) * CONV3 ; E19FWD= HET * 1.1E3 * exp(-TCORR* 0 ) * CONV3 ; E20FWD= HET * 3.2E3 * exp(-TCORR* 0 ) * CONV3 ; E21FWD= HET * 0.18 * exp(-TCORR* (-4030) ) * CONV3 ; E22FWD= HET * 1.4E-4 * exp(-TCORR* (-2500) ) * CONV3 ; E23FWD= HET * 790 * exp(-TCORR* 2500 ) * CONV3 ; E24FWD= HET * 2.5E7 * exp(-TCORR* 2450 ) * CONV3 ; E25FWD= HET * 8.5E9 * exp(-TCORR* 0 ) * CONV3; E26FWD= HET * 1.2E10 * exp(-TCORR* 0 ) * CONV3; E27FWD= HET * 4.3E9 * exp(-TCORR* 0 ) * CONV3; E28FWD= HET * 2.1E10 * exp(-TCORR* 0 ) * CONV3; E29FWD= HET * 1.0E4 * exp(-TCORR* 0 ) * CONV3; E30FWD= HET * 1.1E10 * exp(-TCORR* 0 ) * CONV3; E31FWD= HET * 4.4E10 * exp(-TCORR* 0 ) * CONV3; E32FWD= HET * 1.9E8 * exp(-TCORR* 0 ) * CONV3; E33FWD= HET * 1.86E10 * exp(-TCORR* 0 ) * CONV3 ; E34FWD= HET * 1.5E8 * exp(-TCORR* 0 ) * CONV3 ; E35FWD= HET * 1.1E6 * exp(-TCORR* 0 ) * CONV3 ; E36FWD= HET * 1E7 * exp(-TCORR* 0 ) * CONV3 ; E37FWD= HET * 2E5 * exp(-TCORR* 0 ) * CONV3 ; E38FWD= HET * 3.1E7 * exp(-TCORR* 0 ) * CONV3 ; E39FWD= HET * 3.5E8 * exp(-TCORR* 0 ) * CONV3 ; E40FWD= HET * 330 * exp(-TCORR* 0 ) ; E41FWD= HET * 5E4 * exp(-TCORR* 0 ) ; E42FWD= HET * 5.9E4 * exp(-TCORR* 0 ) ; E43FWD= HET * 7.5E5 * exp(-TCORR* 0 ) ; E44FWD= HET * 3.2E9 * exp(-TCORR* 0 ) ; E45FWD= HET * 2.6E6 * exp(-TCORR* 0 ) ; E46FWD= HET * 6.32E6 * exp(-TCORR* 0 ) ; E47FWD= HET * 21.5 * exp(-TCORR* 0 ) * CONV3 ; E48FWD= HET * 3E-3 * exp(-TCORR* 0 ) ; E49FWD= HET * 3E-3 * exp(-TCORR* 0 ) ; E50FWD= HET * 3E-3 * exp(-TCORR* 0 ) ; E51FWD= HET * 2.52E8 * exp(-TCORR* 0 ) * CONV3 ; E52FWD= HET * 3.4E8 * exp(-TCORR* 4300 ) * CONV3 ; E53FWD= HET * 3E3 * exp(-TCORR* 0 ) * CONV3 ; E54FWD= HET * 315 * exp(-TCORR* 0 ) ; E55FWD= HET * 1.1E7 * exp(-TCORR* 0 ) * CONV3 ; E56FWD= HET * 3.4E4 * exp(-TCORR* 0 ) * CONV3 * CONV3 ; E57FWD= HET * 8.75E5 * exp(-TCORR* (-46) ) ; E58FWD= HET * 8.75E5 * exp(-TCORR* (-46) ) ; E1BCK= HET * 1.3E11 * exp(-TCORR* 0 ) * CONV3 ; E2BCK= HET * 5.6E4 * exp(-TCORR* 8500 ) ; E3BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E4BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E5BCK= HET * 2.9E5 * exp(-TCORR* 0 ) * CONV3 ; E6BCK= HET * 3.4E10 * exp(-TCORR* 0 ) * CONV3 ; E7BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E8BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E9BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E10BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E11BCK= HET * 4.6E-3 * exp(-TCORR* 0 ) ; E12BCK= HET * 2.0E8 * exp(-TCORR* 0 ) * CONV3 ; E13BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E14BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E15BCK= HET * 1E11 * exp(-TCORR* 0 ) * CONV3 ; E16BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E17BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E18BCK= HET * 4.3E8 * exp(-TCORR* 0 ) * CONV3 ; E19BCK= HET * 8.0E9 * exp(-TCORR* 0 ) * CONV3 ; E20BCK= HET * 1.8E5 * exp(-TCORR* 0 ) ; E21BCK= HET * 5.1E-3 * exp(-TCORR* 0 ) ; E22BCK= HET * 5.69E-3 * exp(-TCORR* 0 ) ; E23BCK= HET * 3.95E-6 * exp(-TCORR* 0 ) *CONV3 ; E24BCK= HET * 3.95E-6 * exp(-TCORR* 5530 ) * CONV3 ; E25BCK= HET * 6.0E4 * exp(-TCORR* 0 ) ; E26BCK= HET * 1.9E4 * exp(-TCORR* 0 ) ; E27BCK= HET * 6.1E9 * exp(-TCORR* 0 ) ; E28BCK= HET * 4100 * exp(-TCORR* 0 ) * CONV3 ; E29BCK= HET * 4.5E7 * exp(-TCORR* 0 ) * CONV3 ; E30BCK= HET * 3.3E7 * exp(-TCORR* 0 ) ; E31BCK= HET * 2.45E-2 * exp(-TCORR* 0 ) * CONV3 ; E32BCK= HET * 2.7E6 * exp(-TCORR* 0 ) * CONV3 ; E33BCK= HET * 2.0E10 * exp(-TCORR* 0 ) * CONV3 ; E34BCK= HET * 6.5E3 * exp(-TCORR* 0 ) ; E35BCK= HET * 6.5E6 * exp(-TCORR* 0 ) * CONV3 ; E36BCK= HET * 1E10 * exp(-TCORR* 0 ) * CONV3 ; E37BCK= HET * 2E10 * exp(-TCORR* 0 ) * CONV3 ; E38BCK= HET * 1.033E3 * exp(-TCORR* 0 ) ; E39BCK= HET * 3E4 * exp(-TCORR* 0 ) ; E40BCK= HET * 5.2E10 * exp(-TCORR* 0 ) * CONV3 ; E41BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E42BCK= HET * 4.4E10 * exp(-TCORR* 0 ) * CONV3 ; E43BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E44BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E45BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3 ; E46BCK= HET * 2E10 * exp(-TCORR* 0 ) * CONV3 ; E47BCK= HET * 5.5E-3 * exp(-TCORR* 0 ) ; E48BCK= HET * 7.5E6 * exp(-TCORR* 0 ) * CONV3; E49BCK= HET * 1.89E4 * exp(-TCORR* 0 ) * CONV3; E50BCK= HET * 114 * exp(-TCORR* 0 ) * CONV3; E51BCK= HET * 2.1E8 * exp(-TCORR* 0 ) * CONV3; E52BCK= HET * 1E8 * exp(-TCORR* 0 ) * CONV3; E53BCK= HET * 2160 * exp(-TCORR* 0 ) ; E54BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3; E55BCK= HET * 3E4 * exp(-TCORR* 0 ) ; E56BCK= HET * 1.1E5 * exp(-TCORR* 0 ) * CONV3; E57BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3; E58BCK= HET * 5E10 * exp(-TCORR* 0 ) * CONV3; ** ; COMPILE INTEGRATE; KCAP102t=0 ; KCAP103t=0 ; KCAP104t=0 ; KCAP106t=0 ; KCAP107t=0 ; KCAP112t=0 ; KCAP113t=0 ; KCAP118t=0 ; KCAP118t=0 ; KCAP119t=0 ; KCAP119t=0 ; KCAP123t=0 ; KCAP125t=0 ; KCAP126t=0 ; KCAP128t=0 ; KCAP129t=0 ; KCAP12t=0 ; KCAP130t=0 ; KCAP131t=0 ; KCAP132t=0 ; KCAP133t=0 ; KCAP134t=0 ; KCAP135t=0 ; KCAP136t=0 ; KCAP136t=0 ; KCAP137t=0 ; KCAP137t=0 ; KCAP138t=0 ; KCAP138t=0 ; KCAP139t=0 ; KCAP139t=0 ; KCAP140t=0 ; KCAP140t=0 ; KCAP141t=0 ; KCAP142t=0 ; KCAP144t=0 ; KCAP145t=0 ; KCAP146t=0 ; KCAP14t=0 ; KCAP153t=0 ; KCAP154t=0 ; KCAP155t=0 ; KCAP162t=0 ; KCAP163t=0 ; KCAP164t=0 ; KCAP169t=0 ; KCAP170t=0 ; KCAP171t=0 ; KCAP172t=0 ; KCAP173t=0 ; KCAP174t=0 ; KCAP17t=0 ; KCAP181t=0 ; KCAP182t=0 ; KCAP183t=0 ; KCAP184t=0 ; KCAP185t=0 ; KCAP186t=0 ; KCAP193t=0 ; KCAP194t=0 ; KCAP195t=0 ; KCAP196t=0 ; KCAP197t=0 ; KCAP198t=0 ; KCAP1t=0 ; KCAP217t=0 ; KCAP218t=0 ; KCAP221t=0 ; KCAP222t=0 ; KCAP223t=0 ; KCAP226t=0 ; KCAP228t=0 ; KCAP229t=0 ; KCAP232t=0 ; KCAP235t=0 ; KCAP236t=0 ; KCAP239t=0 ; KCAP25t=0 ; KCAP261t=0 ; KCAP267t=0 ; KCAP268t=0 ; KCAP269t=0 ; KCAP278t=0 ; KCAP285t=0 ; KCAP286t=0 ; KCAP287t=0 ; KCAP288t=0 ; KCAP289t=0 ; KCAP28t=0 ; KCAP290t=0 ; KCAP29t=0 ; KCAP30t=0 ; KCAP310t=0 ; KCAP31t=0 ; KCAP32t=0 ; KCAP32t=0 ; KCAP33t=0 ; KCAP34t=0 ; KCAP39t=0 ; KCAP41t=0 ; KCAP57t=0 ; KCAP5t=0 ; KCAP61t=0 ; KCAP61t=0 ; KCAP62t=0 ; KCAP63t=0 ; KCAP64t=0 ; KCAP65t=0 ; KCAP66t=0 ; KCAP67t=0 ; KCAP68t=0 ; KCAP69t=0 ; KCAP70t=0 ; KCAP71t=0 ; KCAP71t=0 ; KCAP72t=0 ; KCAP76t=0 ; KCAP77t=0 ; KCAP78t=0 ; KCAP83t=0 ; KCAP84t=0 ; KCAP90t=0 ; KCAP91t=0 ; KCAP92t=0 ; KCAP96t=0 ; **; COMPILE DAY; ON=1; **; COMPILE NIGHT; ON=0; **; COMPILE GENERAL; DayReal=Day0+(TIME/(24*3600)); TETA=2*PI*DayReal/365.0 ; Solardecli=0.006918-0.399912*cos(TETA)+0.070257*sin(TETA) -0.0067518*cos(2*TETA)+0.000907*sin(2*TETA) -0.002697*cos(3*TETA)+0.00148*sin(3*TETA) ; EQT=0.000075+0.001868*cos(TETA)-0.030277*sin(TETA) -0.01465*cos(2*TETA)-0.040849*sin(2*TETA) ; LAT=(TIME/3600)+LONGCOR+EQT ; Local apparent time Hourangle=(LAT-12)*(15)*PI/180 ; z=(sin(LATITU)*sin(solardecli)+cos(LATITU)*cos(solarDecli) *cos(Hourangle)); sza=artan(sqrt(1-z**2)/z); IF (sza) *,*,1; sza=sza+PI; LABEL 1; photon=sza*180/PI; start=TIME/3600; IF (sza*0.78-1.570796327) *,10,10; J1 = ON* 7.67e-3 * exp( 1.773179 * (1-(1/cos( 0.77233 * sza)))) ; J2 = ON* 2.55e-5 * exp( 4.11188 * (1-(1/cos( 0.77335 * sza)))) ; J4 = ON* 1.72e-3 * exp( 1.86895 * (1-(1/cos( 0.76528 * sza)))) ; J5 = ON* 4.67e-7 * exp( 3.00003 * (1-(1/cos( 0.77346 * sza)))) ; J6 = ON* 2.93e-6 * exp( 2.74497 * (1-(1/cos( 0.77084 * sza)))) ; J9 = ON* 7.65e-6 * exp( 2.35766 * (1-(1/cos( 0.76470 * sza)))) ; J10 = ON* 3.3e-5 * exp( 2.16815 * (1-(1/cos( 0.76255 * sza)))) ; J11 = ON* 2.85e-5 * exp( 2.49092 * (1-(1/cos( 0.77402 * sza)))) ; J12 = ON* 4.45e-6 * exp( 2.88929 * (1-(1/cos( 0.77445 * sza)))) ; J13 = ON* 4.38e-6 * exp( 2.32912 * (1-(1/cos( 0.76201 * sza)))) ; J14 = ON* 4.38e-6 * exp( 2.32912 * (1-(1/cos( 0.76201 * sza)))) ; J15 = ON* 4.58e-9 * exp( 9.3050 * (1-(1/cos( 0.71754 * sza)))) ; J16 = ON* 9.87e-6 * exp( 2.90556 * (1-(1/cos( 0.77465 * sza)))) ; J20 = ON* 4.48e-4 * exp( 2.13016 * (1-(1/cos( 0.76047 * sza)))) ; J21 = ON* 8.26e-7 * exp( 2.83477 * (1-(1/cos( 0.77775 * sza)))) ; J22 = ON* 4.46e-4 * exp( 2.05326 * (1-(1/cos( 0.76199 * sza)))) ; J23 = ON* 9.87e-6 * exp( 2.90556 * (1-(1/cos( 0.77465 * sza)))) ; P1 = HET*ON*7.64E-6* EXP( 2.46425 * (1-(1/cos( 0.76355 * sza)))) ; P2 = HET*ON*4.76E-3* EXP( 2.19894 * (1-(1/cos( 0.76087 * sza)))) ; P3 = HET*ON*6.10E-3* EXP( 2.21811 * (1-(1/cos( 0.76118 * sza)))) ; P4 = HET*ON*6.80E-3* EXP( 2.24656 * (1-(1/cos( 0.76118 * sza)))) ; P5 = HET*ON*2.70E-5* EXP( 2.01536 * (1-(1/cos( 0.75915 * sza)))) ; P6 = HET*ON*4.57E-7* EXP( 2.59408 * (1-(1/cos( 0.77213 * sza)))) ; P7 = HET*ON*2.47E-2* EXP( 1.95825 * (1-(1/cos( 0.76782 * sza)))) ; P8 = HET*ON*1.55E-2* EXP( 1.98840 * (1-(1/cos( 0.76679 * sza)))) ; P9 = HET*ON*7.64E-6* EXP( 2.46425 * (1-(1/cos( 0.76355 * sza)))) ; LABEL 10 ; IF (sza*0.81-1.570796327) *,20,20; J17 = ON* 2.32e-3 * exp( 1.32625 * (1-(1/cos( 0.80629 * sza)))) ; J18 = ON* 2.32e-3 * exp( 1.32625 * (1-(1/cos( 0.80629 * sza)))) ; J19 = ON* 2.5e-3 * exp( 1.38450 * (1-(1/cos( 0.80103 * sza)))) ; LABEL 20 ; IF (sza*0.89103-1.570796327) *,30,30; * gas phase photolysis rates ; J3 = ON* 4.38e-4 * exp( 0.74524 * (1-(1/cos( 0.87973 * sza)))) ; J7 = ON* 2.32e-2 * exp( 0.66942 * (1-(1/cos( 0.89103 * sza)))) ; J8 = ON* 2.01e-1 * exp( 0.79342 * (1-(1/cos( 0.87147 * sza)))) ; P10 = HET * J7; P11 = HET * J8; LABEL 30 ; ** ; COMPILE EQUATIONS; * *****************************************************************; * Phase transfer ; * *****************************************************************; * ; * ******************************************; 'O3a='O3a+KMTO3*O3g-(KMTO3/(HCO3*R*TEMP))*O3a; 'O3g='O3g-(KMTO3*L)*O3g+(KMTO3*L/(HCO3*R*TEMP))*O3a; 'O2a='O2a+KMTO2*O2g-(KMTO2/(HCO2*R*TEMP))*O2a; 'O2g='O2g-(KMTO2*L)*O2g+(KMTO2*L/(HCO2*R*TEMP))*O2a; 'OHa='OHa+KMTOH*OHg-(KMTOH/(HCOH*R*TEMP))*OHa; 'OHg='OHg-(KMTOH*L)*OHg+(KMTOH*L/(HCOH*R*TEMP))*OHa; 'HO2a='HO2a+KMTHO2*HO2g-(KMTHO2/(HCHO2*R*TEMP))*HO2a; 'HO2g='HO2g-(KMTHO2*L)*HO2g+(KMTHO2*L/(HCHO2*R*TEMP))*HO2a; 'H2O2a='H2O2a+KMTH2O2*H2O2g-(KMTH2O2/(HCH2O2*R*TEMP))*H2O2a; 'H2O2g='H2O2g-(KMTH2O2*L)*H2O2g+(KMTH2O2*L/(HCH2O2*R*TEMP))*H2O2a; 'NO3a='NO3a+KMTNO3*NO3g-(KMTNO3/(HCNO3*R*TEMP))*NO3a; 'NO3g='NO3g-(KMTNO3*L)*NO3g+(KMTNO3*L/(HCNO3*R*TEMP))*NO3a; * ****************************; * N2O5 complete mechanism; 'N2O5a='N2O5a+KMTN2O5*N2O5g-(KMTN2O5/(HCN2O5*R*TEMP))*N2O5a; 'N2O5g='N2O5g-(KMTN2O5*L)*N2O5g+(KMTN2O5*L/(HCN2O5*R*TEMP))*N2O5a; *'NO2p='NO2p+KMTN2O5*N2O5g; *'NO3m='NO3m+KMTN2O5*N2O5g; *'N2O5g='N2O5g-(KMTN2O5*L)*N2O5g; * ********************; 'NO2a='NO2a+KMTNO2*NO2g-(KMTNO2/(HCNO2*R*TEMP))*NO2a; 'NO2g='NO2g-(KMTNO2*L)*NO2g+(KMTNO2*L/(HCNO2*R*TEMP))*NO2a; 'HNO3a='HNO3a+KMTHNO3*HNO3g-(KMTHNO3/(HCHNO3*R*TEMP))*HNO3a; 'HNO3g='HNO3g-(KMTHNO3*L)*HNO3g+(KMTHNO3*L/(HCHNO3*R*TEMP))*HNO3a; 'HNO4a='HNO4a+KMTHNO4*HNO4g-(KMTHNO4/(HCHNO4*R*TEMP))*HNO4a; 'HNO4g='HNO4g-(KMTHNO4*L)*HNO4g+(KMTHNO4*L/(HCHNO4*R*TEMP))*HNO4a; 'HNO2a='HNO2a+KMTHNO2*HNO2g-(KMTHNO2/(HCHNO2*R*TEMP))*HNO2a; 'HNO2g='HNO2g-(KMTHNO2*L)*HNO2g+(KMTHNO2*L/(HCHNO2*R*TEMP))*HNO2a; 'NH3a='NH3a+KMTNH3*NH3g-(KMTNH3/(HCNH3*R*TEMP))*NH3a; 'NH3g='NH3g-(KMTNH3*L)*NH3g+(KMTNH3*L/(HCNH3*R*TEMP))*NH3a; 'CH3O2Ha='CH3O2Ha+KMTOP1*OP1g- (KMTOP1/(HCOP1*R*TEMP))*CH3O2Ha; 'OP1g='OP1g- (KMTOP1*L)*OP1g+(KMTOP1*L/(HCOP1*R*TEMP))*CH3O2Ha; 'HCHOa='HCHOa+KMTHCHO*HCHOg-(KMTHCHO/(HCHCHO*R*TEMP))*HCHOa; 'HCHOg='HCHOg-(KMTHCHO*L)*HCHOg+(KMTHCHO*L/(HCHCHO*R*TEMP))*HCHOa; 'HCOOHa='HCOOHa+KMTORA1*ORA1g- (KMTORA1/(HCORA1*R*TEMP))*HCOOHa; 'ORA1g='ORA1g- (KMTORA1*L)*ORA1g+(KMTORA1*L/(HCORA1*R*TEMP))*HCOOHa; 'CO2a='CO2a+KMTCO2*CO2g-(KMTCO2/(HCCO2*R*TEMP))*CO2a; 'CO2g='CO2g-(KMTCO2*L)*CO2g+(KMTCO2*L/(HCCO2*R*TEMP))*CO2a; 'HCLa='HCLa+KMTHCL*HCLg-(KMTHCL/(HCHCL*R*TEMP))*HCLa; 'HCLg='HCLg-(KMTHCL*L)*HCLg+(KMTHCL*L/(HCHCL*R*TEMP))*HCLa; * ********************; 'CL2a='CL2a+KMTCL2*CL2g-(KMTCL2/(HCCL2*R*TEMP))*CL2a; 'CL2g='CL2g-(KMTCL2*L)*CL2g+(KMTCL2*L/(HCCL2*R*TEMP))*CL2a; * ********************; 'BR2a='BR2a+KMTBR2*BR2g-(KMTBR2/(HCBR2*R*TEMP))*BR2a; 'BR2g='BR2g-(KMTBR2*L)*BR2g+(KMTBR2*L/(HCBR2*R*TEMP))*BR2a; 'BRCLa='BRCLa+KMTBRCL*BRCLg-(KMTBRCL/(HCBRCL*R*TEMP))*BRCLa; 'BRCLg='BRCLg-(KMTBRCL*L)*BRCLg+(KMTBRCL*L/(HCBRCL*R*TEMP))*BRCLa; 'BRNO2a='BRNO2a+KMTBRNO2*BRNO2g- (KMTBRNO2/(HCBRNO2*R*TEMP))*BRNO2a; 'BRNO2g='BRNO2g- (KMTBRNO2*L)*BRNO2g+(KMTBRNO2*L/(HCBRNO2*R*TEMP))*BRNO2a; 'CLNO2a='CLNO2a+KMTCLNO2*CLNO2g-(KMTCLNO2/(HCCLNO2*R*TEMP))*CLNO2a; 'CLNO2g='CLNO2g- (KMTCLNO2*L)*CLNO2g+(KMTCLNO2*L/(HCCLNO2*R*TEMP))*CLNO2a; 'SO2a='SO2a+KMTSO2*SO2g-(KMTSO2/(HCSO2*R*TEMP))*SO2a; 'SO2g='SO2g-(KMTSO2*L)*SO2g+(KMTSO2*L/(HCSO2*R*TEMP))*SO2a; * ******************************************. * ********************; 'ACO3a='ACO3a+KMTACO3*ACO3g-(KMTACO3/(HCACO3*R*TEMP))*ACO3a; 'ACO3g='ACO3g-(KMTACO3*L)*ACO3g+(KMTACO3*L/(HCACO3*R*TEMP))*ACO3a; 'CH3CHOa='CH3CHOa+KMTALD*ALDg-(KMTALD/(HCALD*R*TEMP))*CH3CHOa; 'ALDg='ALDg-(KMTALD*L)*ALDg+(KMTALD*L/(HCALD*R*TEMP))*CH3CHOa; 'CH3OHa='CH3OHa+KMTCH3OH*CH3OHg-(KMTCH3OH/(HCCH3OH*R*TEMP))*CH3OHa; 'CH3OHg='CH3OHg-(KMTCH3OH*L)*CH3OHg+(KMTCH3OH*L/(HCCH3OH*R*TEMP)) *CH3OHa; 'ETOHa='ETOHa+KMTETOH*ETOHg-(KMTETOH/(HCETOH*R*TEMP))*ETOHa; 'ETOHg='ETOHg-(KMTETOH*L)*ETOHg+(KMTETOH*L/(HCETOH*R*TEMP))*ETOHa; 'CH3COOHa='CH3COOHa+KMTORA2*ORA2g-(KMTORA2/(HCORA2*R*TEMP))*CH3COOHa; 'ORA2g='ORA2g-(KMTORA2*L)*ORA2g+(KMTORA2*L/(HCORA2*R*TEMP))*CH3COOHa; 'H2SO4a='H2SO4a+KMTSULF*SULFg-(KMTSULF/(HCSULF*R*TEMP))*H2SO4a; 'SULFg='SULFg-(KMTSULF*L)*SULFg+(KMTSULF*L/(HCSULF*R*TEMP))*H2SO4a; 'GLYa='GLYa+KMTGLY*GLYg-(KMTGLY/(HCGLY*R*TEMP))*GLYa; 'GLYg='GLYg-(KMTGLY*L)*GLYg+(KMTGLY*L/(HCGLY*R*TEMP))*GLYa; 'CH3O2a='CH3O2a+KMTMO2*MO2g-(KMTMO2/(HCMO2*R*TEMP))*CH3O2a; 'MO2g='MO2g-(KMTMO2*L)*MO2g+(KMTMO2*L/(HCMO2*R*TEMP))*CH3O2a; 'PAAa='PAAa+KMTPAA*PAAg-(KMTPAA/(HCPAA*R*TEMP))*PAAa; 'PAAg='PAAg-(KMTPAA*L)*PAAg+(KMTPAA*L/(HCPAA*R*TEMP))*PAAa; 'ETHPa='ETHPa+KMTETHP*ETHPg-(KMTETHP/(HCETHP*R*TEMP))*ETHPa; 'ETHPg='ETHPg-(KMTETHP*L)*ETHPg+(KMTETHP*L/(HCETHP*R*TEMP))*ETHPa; 'NOa='NOa+KMTNO*NOg-(KMTNO/(HCNO*R*TEMP))*NOa; 'NOg='NOg-(KMTNO*L)*NOg+(KMTNO*L/(HCNO*R*TEMP))*NOa; * ******************************************; * ****************************************************************; * Emissions; * ****************************************************************; * ; % EMICH3OH : = CH3OHg ; % EMIETOH : = ETOHg ; % EMIHC3 : = HC3g ; % EMIXYL : = XYLg ; % EMIKET : = KETg ; % EMIHC5 : = HC5g ; % EMIETH : = ETHg ; % EMIOLT : = OLTg ; % EMIALD : = ALDg ; % EMIHCHO : = HCHOg ; % EMITOL : = TOLg ; % EMIDIEN : = DIENg ; % EMICSL : = CSLg ; % EMIETE : = ETEg ; % EMIHC8 : = HC8g ; % EMICO : = COg ; % EMIORA2 : = ORA2g ; % EMISO2 : = SO2g ; % EMINH3 : = NH3g ; % EMINO : = NOg ; % EMIAPI : = APIg ; % EMILIM : = LIMg ; % EMIDMS : = DMSg ; % EMIISO : = ISOg ; * ; * Gas phase reactions according to RACM; * ; % J1 : NO2g = NOg + O3Pg ; % J2 : O3g = O1Dg + O2g ; % J3 : O3g = O3Pg + O2g ; % J4 : HNO2g = NOg + OHg ; % J5 : HNO3g = NO2g + OHg ; % J6*0.65: HNO4g = NO2g + HO2g ; % J6*0.35 : HNO4g = NO3g + OHg ; % J7 : NO3g = NOg + O2g ; % J8 : NO3g = NO2g + O3Pg ; % J9 : H2O2g = OHg + OHg ; % J10 : HCHOg = COg + H2g ; % J11 : HCHOg = COg + HO2g + HO2g ; % J12 : ALDg = COg + HO2g + MO2g ; % J13 : OP1g = OHg + HO2g + HCHOg ; % J14 : OP2g = OHg + HO2g + ALDg ; % J15 : PAAg = OHg + MO2g ; % J16 : KETg = ETHPg + ACO3g ; % J17 : GLYg = ; % 1.87 * J17 * GLYg : = COg ; % 0.87 * J17 * GLYg : = H2g ; % 0.13 * J17 * GLYg : = HCHOg ; % J18 : GLYg = ; % 1.55 * J18 * GLYg : = COg ; % 0.80 * J18 * GLYg : = HO2g ; % 0.45 * J18 * GLYg : = HCHOg ; % J18 * 0.15 * GLYg : = H2g ; % J19 : MGLYg = COg + HO2g + ACO3g ; % J20 : DCBg = HO2g + TCO3g ; % J21 : ONITg = NO2g + HO2g ; % 0.2 * J21 * ONITg : = ALDg ; % 0.8 * J21 * ONITg : = KETg ; % J22 : MACRg = COg + HO2g + HCHOg +ACO3g ; % J23 : HKETg = HO2g + HCHOg + ACO3g ; % K24 : O3Pg + O2g = O3g ; % K25 : O3g + O3Pg = O2g + O2g ; % K26 : O1Dg + N2g = O3Pg + N2g ; % K27 : O1Dg + O2g = O3Pg + O2g ; % K28 : O1Dg + H2Og = OHg + OHg ; % K29 : O3g + OHg = HO2g + O2g ; % K30 : O3g + HO2g = OHg + O2g + O2g ; % K31 : OHg + HO2g = H2Og + O2g ; % K32 : H2O2g + OHg = HO2g + H2Og ; % K33 : HO2g + HO2g = H2O2g + O2g ; % K34 : HO2g + HO2g + H2Og = H2O2g + H2Og + O2g ; % K35 : NOg + O3Pg = NO2g ; % K36 : NO2g + O3Pg = NOg + O2g ; % K37 : NO2g + O3Pg = NO3g ; % K38 : NOg + OHg = HNO2g ; % K39 : NO2g + OHg = HNO3g ; % K40 : NO3g + OHg = NO2g + HO2g ; % K41 : NOg + HO2g = NO2g + OHg ; % K42 : NO2g + HO2g = HNO4g ; % K43 : HNO4g = NO2g + HO2g ; % 0.70 * K44 : NO3g + HO2g = NO2g + OHg + O2g; % 0.30 * K44 : NO3g + HO2g = HNO3g + O2g; % K45 : HNO2g + OHg = NO2g + H2Og ; % K46 : HNO3g + OHg = NO3g + H2Og ; % K47 : HNO4g + OHg = NO2g + H2Og + O2g ; % K48 : O3g + NOg = NO2g + O2g ; % K49 : O3g + NO2g = NO3g + O2g ; % K50 : NOg + NOg + O2g = NO2g + NO2g ; % K51 : NOg + NO3g = NO2g + NO2g ; % K52 : NO2g + NO3g = NOg + NO2g + O2g ; % K53 : NO2g + NO3g = N2O5g ; % K54 : N2O5g = NO2g + NO3g ; % K55 : NO3g + NO3g = NO2g + NO2g + O2g ; % K56 : H2g + OHg = HO2g + H2Og ; % K57 : SO2g + OHg = SULFg + HO2g ; % K58 : COg + OHg = HO2g + CO2g ; % K59 : ISOg + O3Pg = ; % 0.86 * K59 * ISOg * O3Pg : = OLTg ; % 0.05 * K59 * ISOg * O3Pg : = HCHOg ; % 0.02 * K59 * ISOg * O3Pg : = OHg ; % 0.01 * K59 * ISOg * O3Pg : = COg ; % 0.13 * K59 * ISOg * O3Pg : = DCBg ; % 0.28 * K59 * ISOg * O3Pg : = HO2g ; % 0.15 * K59 * ISOg * O3Pg : = XO2g ; % K60 : MACRg + O3Pg = ALDg ; % K61 : CH4g + OHg = MO2g + H2Og ; % K62 : ETHg + OHg = ETHPg + H2Og ; % 0.583 * K63 * HC3g * OHg : = HC3Pg ; % 0.381 * K63 * HC3g * OHg : = HO2g ; % 0.036 * K63 * HC3g * OHg : = COg+ OHg + ORA1g + GLYg ; % 0.010 * K63 * HC3g * OHg : = HCHOg ; % 0.335 * K63 * HC3g * OHg : = ALDg ; % K63 : HC3g + OHg = H2Og ; % 0.75 * K64 : HC5g + OHg = HC5Pg ; % 0.25 * K64 : HC5g + OHg = HO2g + KETg ; % K64 * HC5g * OHg : = H2Og ; % 0.951 * K65 * HC8g * OHg : = HC8Pg ; % 0.025 * K65 * HC8g * OHg : = ALDg ; % 0.024 * K65 * HC8g * OHg : = HKETg ; % 0.049 * K65 * HC8g * OHg : = HO2g ; % K65 : HC8g + OHg = H2Og ; % K66 : ETEg + OHg = ETEPg ; % K67 : OLTg + OHg = OLTPg ; % K68 : OLIg + OHg = OLIPg ; % K69 : DIENg + OHg = ISOPg ; % K70 : ISOg + OHg = ISOPg ; % K71 : APIg + OHg = APIPg ; % K72 : LIMg + OHg = LIMPg ; % 0.9 * K73 : TOLg + OHg = ADDTg ; % 0.1 * K73 : TOLg + OHg = XO2g + HO2g ; % 0.9 * K74 : XYLg + OHg = ADDXg ; % 0.1 * K74 : XYLg + OHg = XO2g + HO2g ; % 0.85 * K75 : CSLg + OHg = ADDCg ; % 0.1 * K75 : CSLg + OHg = PHOg ; % 0.05 * K75 : CSLg + OHg = HO2g+ XO2g ; % K76 : HCHOg + OHg = HO2g + COg + H2Og ; % K77 : ALDg + OHg = ACO3g+ H2Og ; % K78 : KETg + OHg = KETPg+ H2Og ; % K79 : HKETg + OHg = HO2g + MGLYg + H2Og ; % K80 : GLYg + OHg = HO2g + COg + COg +H2Og ; % K81 : MGLYg + OHg = ACO3g + COg + H2Og ; % K82 : MACRg + OHg = ; % 0.51 * K82 * MACRg * OHg : = TCO3g ; % 0.49 * K82 * MACRg * OHg : = XO2g + HO2g ; % 0.41 * K82 * MACRg * OHg : = COg + HKETg ; % 0.08 * K82 * MACRg * OHg : = HCHOg + MGLYg ; % 0.5 * K83 : OHg + DCBg = TCO3g + HO2g + XO2g ; % 0.35 * K83 : OHg + DCBg = UDDg ; % 0.15 * K83 : OHg + DCBg = GLYg + MGLYg ; % K84 : OHg + UDDg = HO2g ; % 0.88 * K84 * OHg * UDDg : = ALDg ; % 0.12 * K84 * OHg * UDDg : = KETg ; % 0.65 * K85 : OP1g + OHg = MO2g ; % 0.35 * K85 : OP1g + OHg = OHg + HCHOg ; % K86 : OP2g + OHg = ; % 0.44 * K86 * OP2g * OHg : = HC3Pg ; % 0.08 * K86 * OP2g * OHg : = ALDg ; % 0.41 * K86 * OP2g * OHg : = KETg ; % 0.49 * K86 * OP2g * OHg : = OHg ; % 0.07 * K86 * OP2g * OHg : = XO2g ; % 0.65 * K87 : PAAg + OHg = ACO3g ; % 0.35 * K87 : PAAg + OHg = XO2g + HO2g + HCHOg ; % K88 : PANg + OHg = XO2g + NO3g + HCHOg + H2Og ; % K89 : TPANg + OHg = XO2g ; % 0.6 * K89 * TPANg * OHg : = NO3g + HKETg ; % 0.4 * K89 * TPANg * OHg : = HO2g + HCHOg + PANg ; % K90 : ONITg + OHg = HC3Pg + NO2g + H2Og ; % K91 : HCHOg + NO3g = HO2g + COg + HNO3g ; % K92 : ALDg + NO3g = ACO3g + HNO3g ; % K93 : GLYg + NO3g = HO2g + COg + COg + HNO3g ; % K94 : MGLYg + NO3g = ACO3g + COg + HNO3g ; % 0.2 * K95 : MACRg + NO3g = TCO3g + HNO3g ; % 0.8 * K95 : MACRg + NO3g = OLNNg + COg ; % K96 : DCBg + NO3g = ; % 0.5 * K96 * DCBg * NO3g : = TCO3g + HO2g + XO2g+ NO2g+ HNO3g ; % 0.25 * K96 * DCBg * NO3g : = GLYg + ALDg + MGLYg ; % 0.03 * K96 * DCBg * NO3g : = KETg ; % K97 : CSLg + NO3g = HNO3g + PHOg ; % 0.8 * K98 : ETEg + NO3g = OLNNg ; % 0.2 * K98 : ETEg + NO3g = OLNDg ; % 0.43 * K99 : OLTg + NO3g = OLNNg ; % 0.57 * K99 : OLTg + NO3g = OLNDg ; % 0.11 * K100 : OLIg + NO3g = OLNNg ; % 0.89 * K100 : OLIg + NO3g = OLNDg ; % 0.9 * K101 : DIENg + NO3g = OLNNg + MACRg ; % 0.1 * K101 : DIENg + NO3g = OLNDg ; % 0.9 * K102 : ISOg + NO3g = OLNNg + MACRg ; % 0.1 * K102 : ISOg + NO3g = OLNDg ; % 0.9 * K103 : APIg + NO3g = OLNDg ; % 0.1 * K103 : APIg + NO3g = OLNNg ; % 0.13 * K104 : LIMg + NO3g = OLNNg ; % 0.87 * K104 : LIMg + NO3g = OLNDg ; % K105 : TPANg + NO3g = XO2g ; % 0.6 * K105 * TPANg * NO3g : = NO3g + ONITg ; % 0.4 * K105 * TPANg * NO3g : = PANg + HCHOg + NO2g ; % K106 : ETEg + O3g = HCHOg ; % 0.26 * K106 * ETEg * O3g : = HO2g ; % 0.43 * K106 * ETEg * O3g : = COg ; % 0.37 * K106 * ETEg * O3g : = ORA1g ; % 0.13 * K106 * ETEg * O3g : = H2g ; % 0.12 * K106 * ETEg * O3g : = OHg ; % K107 : OLTg + O3g = ; % 0.64 * K107 * OLTg * O3g : = HCHOg ; % 0.44 * K107 * OLTg * O3g : = ALDg ; % 0.37 * K107 * OLTg * O3g : = COg ; % 0.14 * K107 * OLTg * O3g : = ORA1g ; % 0.10 * K107 * OLTg * O3g : = ORA2g + ETHPg ; % 0.25 * K107 * OLTg * O3g : = HO2g ; % 0.40 * K107 * OLTg * O3g : = OHg ; % 0.03 * K107 * OLTg * O3g : = KETg + KETPg + ETHg ; % 0.06 * K107 * OLTg * O3g : = CH4g ; % 0.05 * K107 * OLTg * O3g : = H2g ; % 0.006 * K107 * OLTg * O3g : = H2O2g ; % 0.19 * K107 * OLTg * O3g : = MO2g ; % K108 : OLIg + O3g = ; % 0.02 * K108 * OLIg * O3g : = HCHOg ; % 0.99 * K108 * OLIg * O3g : = ALDg ; % 0.16 * K108 * OLIg * O3g : = KETg ; % 0.3 * K108 * OLIg * O3g : = COg ; % 0.011 * K108 * OLIg * O3g : = H2O2g ; % 0.14 * K108 * OLIg * O3g : = ORA2g ; % 0.07 * K108 * OLIg * O3g : = CH4g ; % 0.22 * K108 * OLIg * O3g : = HO2g ; % 0.63 * K108 * OLIg * O3g : = OHg ; % 0.23 * K108 * OLIg * O3g : = MO2g ; % 0.12 * K108 * OLIg * O3g : = KETPg ; % 0.06 * K108 * OLIg * O3g : = ETHg ; % 0.18 * K108 * OLIg * O3g : = ETHPg ; % K109 : DIENg + O3g = ; % 0.9 * K109 * DIENg * O3g : = HCHOg ; % 0.39 * K109 * DIENg * O3g : = MACRg ; % 0.36 * K109 * DIENg * O3g : = COg ; % 0.15 * K109 * DIENg * O3g : = ORA1g ; % 0.09 * K109 * DIENg * O3g : = O3Pg ; % 0.3 * K109 * DIENg * O3g : = HO2g ; % 0.35 * K109 * DIENg * O3g : = OLTg ; % 0.28 * K109 * DIENg * O3g : = OHg ; % 0.05 * K109 * DIENg * O3g : = H2g ; % 0.15 * K109 * DIENg * O3g : = ACO3g ; % 0.03 * K109 * DIENg * O3g : = MO2g ; % 0.02 * K109 * DIENg * O3g : = KETPg ; % 0.13 * K109 * DIENg * O3g : = XO2g ; % 0.001 * K109 * DIENg * O3g : = H2O2g ; % K110 : ISOg + O3g = ; % 0.9 * K110 * ISOg * O3g : = HCHOg ; % 0.39 * K110 * ISOg * O3g : = MACRg ; % 0.36 * K110 * ISOg * O3g : = COg ; % 0.15 * K110 * ISOg * O3g : = ORA1g ; % 0.09 * K110 * ISOg * O3g : = O3Pg ; % 0.3 * K110 * ISOg * O3g : = HO2g ; % 0.35 * K110 * ISOg * O3g : = OLTg ; % 0.28 * K110 * ISOg * O3g : = OHg ; % 0.05 * K110 * ISOg * O3g : = H2g ; % 0.15 * K110 * ISOg * O3g : = ACO3g ; % 0.03 * K110 * ISOg * O3g : = MO2g ; % 0.02 * K110 * ISOg * O3g : = KETPg ; % 0.13 * K110 * ISOg * O3g : = XO2g ; % 0.001 * K110 * ISOg * O3g : = H2O2g ; % K111 : APIg + O3g = ; % 0.65 * K111 * APIg * O3g : = ALDg ; % 0.53 * K111 * APIg * O3g : = KETg ; % 0.14 * K111 * APIg * O3g : = COg ; % 0.2 * K111 * APIg * O3g : = ETHPg ; % 0.42 * K111 * APIg * O3g : = KETPg ; % 0.85 * K111 * APIg * O3g : = OHg ; % 0.1 * K111 * APIg * O3g : = HO2g ; % 0.02 * K111 * APIg * O3g : = H2O2g ; % K112 : LIMg + O3g = ; % 0.04 * K112 * LIMg * O3g : = HCHOg ; % 0.46 * K112 * LIMg * O3g : = OLTg ; % 0.14 * K112 * LIMg * O3g : = COg ; % 0.16 * K112 * LIMg * O3g : = ETHPg ; % 0.42 * K112 * LIMg * O3g : = KETPg ; % 0.85 * K112 * LIMg * O3g : = OHg ; % 0.1 * K112 * LIMg * O3g : = HO2g ; % 0.02 * K112 * LIMg * O3g : = H2O2g ; % 0.79 * K112 * LIMg * O3g : = MACRg ; % 0.01 * K112 * LIMg * O3g : = ORA1g ; % 0.07 * K112 * LIMg * O3g : = ORA2g ; % K113 : MACRg + O3g = ; % 0.29 * K113 * MACRg * O3g : = HO2g ; % 0.54 * K113 * MACRg * O3g : = COg ; % 0.07 * K113 * MACRg * O3g : = OHg ; % 0.08 * K113 * MACRg * O3g : = H2g ; % 0.40 * K113 * MACRg * O3g : = HCHOg ; % 0.22 * K113 * MACRg * O3g : = ORA1g ; % 0.13 * K113 * MACRg * O3g : = ACO3g + OP2g + ORA2g ; % 0.6 * K113 * MACRg * O3g : = MGLYg ; % K114 : DCBg + O3g = ; % 0.21 * K114 * DCBg * O3g : = OHg + ORA2g ; % 0.29 * K114 * DCBg * O3g : = HO2g ; % 0.66 * K114 * DCBg * O3g : = COg ; % 0.5 * K114 * DCBg * O3g : = GLYg ; % 0.62 * K114 * DCBg * O3g : = MGLYg ; % 0.28 * K114 * DCBg * O3g : = ACO3g ; % 0.16 * K114 * DCBg * O3g : = ALDg ; % 0.11 * K114 * DCBg * O3g : = PAAg + ORA1g ; % K115 : TPANg + O3g = ; % 0.7 * K115 * TPANg * O3g : = NO2g + ACO3g + HCHOg ; % 0.3 * K115 * TPANg * O3g : = PANg ; % 0.13 * K115 * TPANg * O3g : = COg ; % 0.04 * K115 * TPANg * O3g : = H2g ; % 0.11 * K115 * TPANg * O3g : = ORA1g ; % 0.08 * K115 * TPANg * O3g : = HO2g ; % 0.036 * K115 * TPANg * O3g : = OHg ; % K116 : PHOg + NO2g = ONITg ; % 0.1 * K116 * PHOg * NO2g : = CSLg ; % K117 : PHOg + HO2g = CSLg ; % K118 : ADDTg + NO2g = CSLg + HNO2g ; % 0.98 * K119 : ADDTg + O2g = TOLPg ; % 0.02 * K119 : ADDTg + O2g = CSLg + HO2g ; % K120 : ADDTg + O3g = CSLg + OHg ; % K121 : ADDXg + NO2g = CSLg + HNO2g ; % 0.98 * K122 : ADDXg + O2g = XYLPg ; % 0.02 * K122 : ADDXg + O2g = CSLg + HO2g ; % K123 : ADDXg + O3g = CSLg + OHg ; % K124 : ADDCg + NO2g = CSLg + HNO2g ; % 0.98 * K125 : ADDCg + O2g = CSLPg ; % 0.02 * K125 : ADDCg + O2g = CSLg + HO2g ; % K126 : ADDCg + O3g = CSLg + OHg ; % K127 : ACO3g + NO2g = PANg ; % K128 : PANg = ACO3g + NO2g ; % K129 : TCO3g + NO2g = TPANg ; % K130 : TPANg = TCO3g + NO2g ; % K131 : MO2g + NOg = NO2g + HO2g + HCHOg ; % K132 : ETHPg + NOg = NO2g + HO2g + ALDg ; % K133 : HC3Pg + NOg = ; % 0.941 * K133 * HC3Pg * NOg : = NO2g ; % 0.742 * K133 * HC3Pg * NOg : = HO2g ; % 0.048 * K133 * HC3Pg * NOg : = XO2g + ETHPg ; % 0.047 * K133 * HC3Pg * NOg : = HCHOg ; % 0.15 * K133 * HC3Pg * NOg : = MO2g ; % 0.063 * K133 * HC3Pg * NOg : = GLYg ; % 0.623 * K133 * HC3Pg * NOg : = KETg ; % 0.059 * K133 * HC3Pg * NOg : = ONITg ; % 0.233 * K133 * HC3Pg * NOg : = ALDg ; % K134 : HC5Pg + NOg = ; % 0.876 * K134 * HC5Pg * NOg : = NO2g ; % 0.021 * K134 * HC5Pg * NOg : = HCHOg ; % 0.211 * K134 * HC5Pg * NOg : = ALDg ; % 0.722 * K134 * HC5Pg * NOg : = KETg ; % 0.599 * K134 * HC5Pg * NOg : = HO2g ; % 0.031 * K134 * HC5Pg * NOg : = MO2g ; % 0.245 * K134 * HC5Pg * NOg : = ETHPg ; % 0.334 * K134 * HC5Pg * NOg : = XO2g ; % 0.124 * K134 * HC5Pg * NOg : = ONITg ; % K135 : HC8Pg + NOg = ; % 0.15 * K135 * HC8Pg * NOg : = ALDg ; % 0.642 * K135 * HC8Pg * NOg : = KETg ; % 0.133 * K135 * HC8Pg * NOg : = ETHPg ; % 0.261 * K135 * HC8Pg * NOg : = ONITg ; % 0.739 * K135 * HC8Pg * NOg : = NO2g ; % 0.606 * K135 * HC8Pg * NOg : = HO2g ; % 0.416 * K135 * HC8Pg * NOg : = XO2g ; % K136 : ETEPg + NOg = NO2g + HO2g ; % 1.6 * K136 * ETEPg * NOg : = HCHOg ; % 0.2 * K136 * ETEPg * NOg : = ALDg ; % K137 : OLTPg + NOg = NO2g + HCHOg + HO2g ; % 0.94 * K137 * OLTPg * NOg : = ALDg ; % 0.06 * K137 * OLTPg * NOg : = KETg ; % K138 : OLIPg + NOg = NO2g + HO2g ; % 1.71 * K138 * OLIPg * NOg : = ALDg ; % 0.29 * K138 * OLIPg * NOg : = KETg ; % K139 : ISOPg + NOg = ; % 0.847 * K139 * ISOPg * NOg : = NO2g + HO2g ; % 0.446 * K139 * ISOPg * NOg : = MACRg ; % 0.354 * K139 * ISOPg * NOg : = OLTg ; % 0.606 * K139 * ISOPg * NOg : = HCHOg ; % 0.153 * K139 * ISOPg * NOg : = ONITg ; % 0.8 * K140 : APIPg + NOg = HO2g + ALDg + NO2g + KETg ; % 0.2 * K140 : APIPg + NOg = ONITg ; % K141 : LIMPg + NOg = ; % 0.65 * K141 * LIMPg * NOg : = HO2g + NO2g ; % 0.4 * K141 * LIMPg * NOg : = MACRg ; % 0.25 * K141 * LIMPg * NOg : = OLIg + HCHOg ; % 0.35 * K141 * LIMPg * NOg : = ONITg ; % K142 : TOLPg + NOg = ; % 0.95 * K142 * TOLPg * NOg : = NO2g + HO2g ; % 1.2 * K142 * TOLPg * NOg : = GLYg ; % 0.65 * K142 * TOLPg * NOg : = MGLYg ; % 0.5 * K142 * TOLPg * NOg : = DCBg ; % 0.05 * K142 * TOLPg * NOg : = ONITg ; % K143 : XYLPg + NOg = ; % 0.95 * K143 * XYLPg * NOg : = NO2g + HO2g + DCBg ; % 0.6 * K143 * XYLPg * NOg : = MGLYg ; % 0.35 * K143 * XYLPg * NOg : = GLYg ; % 0.05 * K143 * XYLPg * NOg : = ONITg ; % K144 : CSLPg + NOg = GLYg + MGLYg + NO2g+ HO2g ; % K145 : ACO3g + NOg = NO2g + MO2g ; % K146 : TCO3g + NOg = NO2g + HCHOg + ACO3g ; % K147 : KETPg + NOg = ; % 0.54 * K147 * KETPg * NOg : = MGLYg ; % 0.46 * K147 * KETPg * NOg : = ALDg ; % 0.23 * K147 * KETPg * NOg : = ACO3g ; % 0.77 * K147 * KETPg * NOg : = HO2g ; % 0.16 * K147 * KETPg * NOg : = XO2g ; % K147 * KETPg * NOg : = NO2g ; % K148 : OLNNg + NOg = NO2g + HO2g + ONITg ; % K149 : OLNDg + NOg = ; % 0.287 * K149 * OLNDg * NOg : = HCHOg ; % 1.24 * K149 * OLNDg * NOg : = ALDg ; % 0.464 * K149 * OLNDg * NOg : = KETg ; % K149 * OLNDg * NOg : = NO2g + NO2g ; % K150 : MO2g + HO2g = OP1g ; % K151 : ETHPg + HO2g = OP2g ; % K152 : HC3Pg + HO2g = OP2g ; % K153 : HC5Pg + HO2g = OP2g ; % K154 : HC8Pg + HO2g = OP2g ; % K155 : ETEPg + HO2g = OP2g ; % K156 : OLTPg + HO2g = OP2g ; % K157 : OLIPg + HO2g = OP2g ; % K158 : ISOPg + HO2g = OP2g ; % K159 : APIPg + HO2g = OP2g ; % K160 : LIMPg + HO2g = OP2g ; % K161 : TOLPg + HO2g = OP2g ; % K162 : XYLPg + HO2g = OP2g ; % K163 : CSLPg + HO2g = OP2g ; % K164 : ACO3g + HO2g = PAAg ; % K165 : ACO3g + HO2g = O3g + ORA2g ; % K166 : TCO3g + HO2g = OP2g ; % K167 : TCO3g + HO2g = O3g + ORA2g ; % K168 : KETPg + HO2g = OP2g ; % K169 : OLNNg + HO2g = ONITg ; % K170 : OLNDg + HO2g = ONITg ; % K171 : MO2g + MO2g = ; % 0.66 * K171 * MO2g * MO2g : = HO2g ; % 1.33 * K171 * MO2g * MO2g : = HCHOg ; % K172 : ETHPg + MO2g = HO2g ; % 0.75 * K172 * ETHPg * MO2g : = HCHOg + ALDg ; % K173 : HC3Pg + MO2g = ; % 0.81 * K173 * HC3Pg * MO2g : = HCHOg ; % 0.992 * K173 * HC3Pg * MO2g : = HO2g ; % 0.58 * K173 * HC3Pg * MO2g : = ALDg ; % 0.018 * K173 * HC3Pg * MO2g : = KETg ; % 0.007 * K173 * HC3Pg * MO2g : = MO2g ; % 0.005 * K173 * HC3Pg * MO2g : = MGLYg ; % 0.085 * K173 * HC3Pg * MO2g : = XO2g ; % 0.119 * K173 * HC3Pg * MO2g : = GLYg ; % K174 : HC5Pg + MO2g = ; % 0.829 * K174 * HC5Pg * MO2g : = HCHOg ; % 0.946 * K174 * HC5Pg * MO2g : = HO2g ; % 0.523 * K174 * HC5Pg * MO2g : = ALDg ; % 0.24 * K174 * HC5Pg * MO2g : = KETg ; % 0.014 * K174 * HC5Pg * MO2g : = ETHPg ; % 0.049 * K174 * HC5Pg * MO2g : = MO2g ; % 0.245 * K174 * HC5Pg * MO2g : = XO2g ; % K175 : HC8Pg + MO2g = ; % 0.753 * K175 * HC8Pg * MO2g : = HCHOg ; % 0.993 * K175 * HC8Pg * MO2g : = HO2g ; % 0.411 * K175 * HC8Pg * MO2g : = ALDg ; % 0.419 * K175 * HC8Pg * MO2g : = KETg ; % 0.322 * K175 * HC8Pg * MO2g : = XO2g ; % 0.013 * K175 * HC8Pg * MO2g : = ETHPg ; % K176 : ETEPg + MO2g = HO2g ; % 1.55 * K176 * ETEPg * MO2g : = HCHOg ; % 0.35 * K176 * ETEPg * MO2g : = ALDg ; % K177 : OLTPg + MO2g = HO2g ; % 1.25 * K177 * OLTPg * MO2g : = HCHOg ; % 0.669 * K177 * OLTPg * MO2g : = ALDg ; % 0.081 * K177 * OLTPg * MO2g : = KETg ; % K178 : OLIPg + MO2g = HO2g ; % 0.755 * K178 * OLIPg * MO2g : = HCHOg ; % 0.932 * K178 * OLIPg * MO2g : = ALDg ; % 0.313 * K178 * OLIPg * MO2g : = KETg ; % K179 : ISOPg + MO2g = HO2g ; % 1.09 * K179 * ISOPg * MO2g : = HCHOg ; % 0.55 * K179 * ISOPg * MO2g : = MACRg ; % 0.37 * K179 * ISOPg * MO2g : = OLTg ; % 0.08 * K179 * ISOPg * MO2g : = OLIg ; % K180 : APIPg + MO2g = HO2g + HO2g + ALDg + HCHOg + KETg ; % K181 : LIMPg + MO2g = HO2g + HO2g ; % 1.4 * K181 * LIMPg * MO2g : = HCHOg ; % 0.4 * K181 * LIMPg * MO2g : = OLIg ; % 0.6 * K181 * LIMPg * MO2g : = MACRg ; % K182 : TOLPg + MO2g = HCHOg + HO2g + DCBg ; % 0.65 * K182 * TOLPg * MO2g : = GLYg ; % 0.35 * K182 * TOLPg * MO2g : = MGLYg ; % K183 : XYLPg + MO2g = HCHOg + HO2g + DCBg ; % 0.37 * K183 * XYLPg * MO2g : = GLYg ; % 0.63 * K183 * XYLPg * MO2g : = MGLYg ; % K184 : CSLPg + MO2g = GLYg + MGLYg + HO2g + HO2g+ HCHOg ; % K185 : ACO3g + MO2g = HO2g + HCHOg + MO2g ; % K186 : ACO3g + MO2g = HCHOg + ORA2g ; % K187 : TCO3g + MO2g = HCHOg + HCHOg + ACO3g + HO2g ; % K188 : TCO3g + MO2g = HCHOg + ORA2g ; % K189 : KETPg + MO2g = ; % 0.88 * K189 * KETPg * MO2g : = HO2g ; % 0.08 * K189 * KETPg * MO2g : = XO2g ; % 0.75 * K189 * KETPg * MO2g : = HCHOg ; % 0.4 * K189 * KETPg * MO2g : = MGLYg ; % 0.3 * K189 * KETPg * MO2g : = ALDg + HKETg ; % 0.12 * K189 * KETPg * MO2g : = ACO3g ; % K190 : OLNNg + MO2g = HO2g + ONITg ; % 0.75 * K190 * OLNNg * MO2g : = HCHOg ; % K191 : OLNDg + MO2g = ; % 0.5 * K191 * OLNDg * MO2g : = NO2g + HO2g + ONITg ; % 0.96 * K191 * OLNDg * MO2g : = HCHOg ; % 0.64 * K191 * OLNDg * MO2g : = ALDg ; % 0.149 * K191 * OLNDg * MO2g : = KETg ; % K192 : ETHPg + ACO3g = ALDg ; % 0.5 * K192 * ETHPg * ACO3g : = HO2g + MO2g + ORA2g ; % K193 : HC3Pg + ACO3g = ; % 0.488 * K193 * HC3Pg * ACO3g : = HO2g ; % 0.071 * K193 * HC3Pg * ACO3g : = XO2g ; % 0.091 * K193 * HC3Pg * ACO3g : = HCHOg ; % 0.508 * K193 * HC3Pg * ACO3g : = MO2g ; % 0.1 * K193 * HC3Pg * ACO3g : = GLYg ; % 0.006 * K193 * HC3Pg * ACO3g : = ETHPg ; % 0.127 * K193 * HC3Pg * ACO3g : = KETg ; % 0.724 * K193 * HC3Pg * ACO3g : = ALDg ; % 0.004 * K193 * HC3Pg * ACO3g : = MGLYg ; % 0.499 * K193 * HC3Pg * ACO3g : = ORA2g ; % K194 : HC5Pg + ACO3g = ; % 0.677 * K194 * HC5Pg * ACO3g : = ALDg ; % 0.33 * K194 * HC5Pg * ACO3g : = KETg ; % 0.438* K194 * HC5Pg * ACO3g : = HO2g ; % 0.554 * K194 * HC5Pg * ACO3g : = MO2g ; % 0.495 * K194 * HC5Pg * ACO3g : = ORA2g ; % 0.018 * K194 * HC5Pg * ACO3g : = ETHPg ; % 0.237 * K194 * HC5Pg * ACO3g : = XO2g ; % 0.076 * K194 * HC5Pg * ACO3g : = HCHOg ; % K195 : HC8Pg + ACO3g =; % 0.497 * K195 * HC8Pg * ACO3g : = ALDg ; % 0.581 * K195 * HC8Pg * ACO3g : = KETg ; % 0.489 * K195 * HC8Pg * ACO3g : = HO2g ; % 0.507 * K195 * HC8Pg * ACO3g : = MO2g ; % 0.495 * K195 * HC8Pg * ACO3g : = ORA2g ; % 0.015 * K195 * HC8Pg * ACO3g : = ETHPg ; % 0.318 * K195 * HC8Pg * ACO3g : = XO2g ; % K196 : ETEPg + ACO3g = ; % 0.5 * K196 * ETEPg * ACO3g : = HO2g + MO2g+ ORA2g ; % 0.6 * K196 * ETEPg * ACO3g : = ALDg ; % 0.8 * K196 * ETEPg * ACO3g : = HCHOg ; % K197 : OLTPg + ACO3g = ; % 0.501 * K197 * OLTPg * ACO3g : = HO2g + HCHOg + MO2g ; % 0.859 * K197 * OLTPg * ACO3g : = ALDg ; % 0.499 * K197 * OLTPg * ACO3g : = ORA2g ; % 0.141 * K197 * OLTPg * ACO3g : = KETg ; % K198 : OLIPg + ACO3g = ; % 0.941 * K198 * OLIPg * ACO3g : = ALDg ; % 0.569 * K198 * OLIPg * ACO3g : = KETg ; % 0.51 * K198 * OLIPg * ACO3g : = HO2g + MO2g ; % 0.49 * K198 * OLIPg * ACO3g : = ORA2g ; % K199 : ISOPg + ACO3g = ; % 0.34 * K199 * ISOPg * ACO3g : = HCHOg ; % 0.771 * K199 * ISOPg * ACO3g : = MACRg ; % 0.229 * K199 * ISOPg * ACO3g : = OLTg ; % 0.506 * K199 * ISOPg * ACO3g : = HO2g + MO2g ; % 0.494 * K199 * ISOPg * ACO3g : = ORA2g ; % K200 : APIPg + ACO3g = HO2g + ALDg + MO2g+ KETg ; % K201 : LIMPg + ACO3g = HO2g + MO2g ; % 0.4 * K201 * LIMPg * ACO3g : = OLIg+ HCHOg ; % 0.6 * K201 * LIMPg * ACO3g : = MACRg ; % K202 : TOLPg + ACO3g = MO2g + HO2g + DCBg ; % 0.65 * K202 * TOLPg * ACO3g : = GLYg ; % 0.35 * K202 * TOLPg * ACO3g : = MGLYg ; % K203 : XYLPg + ACO3g = MO2g + HO2g + DCBg ; % 0.37 * K203 * XYLPg * ACO3g : = GLYg ; % 0.63 * K203 * XYLPg * ACO3g : = MGLYg ; % K204 : CSLPg + ACO3g = GLYg + MGLYg + HO2g+ MO2g ; % K205 : ACO3g + ACO3g = MO2g + MO2g ; % K206 : TCO3g + ACO3g = HCHOg + MO2g + ACO3g ; % K207 : KETPg + ACO3g = ; % 0.54 * K207 * KETPg * ACO3g : = MGLYg ; % 0.35 * K207 * KETPg * ACO3g : = ALDg ; % 0.11 * K207 * KETPg * ACO3g : = KETg ; % 0.12 * K207 * KETPg * ACO3g : = ACO3g ; % 0.38 * K207 * KETPg * ACO3g : = HO2g ; % 0.08 * K207 * KETPg * ACO3g : = XO2g ; % 0.5 * K207 * KETPg * ACO3g : = MO2g + ORA2g ; % K208 : OLNNg + ACO3g = ONITg ; % 0.5 * K208 * OLNNg * ACO3g : = HO2g + MO2g + ORA2g ; % K209 : OLNDg + ACO3g = ; % 0.207 * K209 * OLNDg * ACO3g : = HCHOg ; % 0.65 * K209 * OLNDg * ACO3g : = ALDg ; % 0.167 * K209 * OLNDg * ACO3g : = KETg ; % 0.484 * K209 * OLNDg * ACO3g : = ORA2g + ONITg ; % 0.516 * K209 * OLNDg * ACO3g : = NO2g + MO2g ; % K210 : OLNNg + OLNNg = ONITg + ONITg + HO2g ; % K211 : OLNNg + OLNDg = ; % 0.202 * K211 * OLNNg * OLNDg : = HCHOg ; % 0.64 * K211 * OLNNg * OLNDg : = ALDg ; % 0.149 * K211 * OLNNg * OLNDg : = KETg ; % 0.5 * K211 * OLNNg * OLNDg : = HO2g + NO2g ; % 1.5 * K211 * OLNNg * OLNDg : = ONITg ; % K212 : OLNDg + OLNDg = ; % 0.504 * K212 * OLNDg * OLNDg : = HCHOg ; % 1.21 * K212 * OLNDg * OLNDg : = ALDg ; % 0.285 * K212 * OLNDg * OLNDg : = KETg ; % K212 * OLNDg * OLNDg : = ONITg + NO2g ; % K213 : MO2g + NO3g = NO2g + HO2g + HCHOg ; % K214 : ETHPg + NO3g = NO2g + HO2g + ALDg ; % K215 : HC3Pg + NO3g = ; % 0.048 * K215 * HC3Pg * NO3g : = HCHOg ; % 0.243 * K215 * HC3Pg * NO3g : = ALDg ; % 0.67 * K215 * HC3Pg * NO3g : = KETg ; % 0.063 * K215 * HC3Pg * NO3g : = GLYg ; % 0.792 * K215 * HC3Pg * NO3g : = HO2g ; % 0.155 * K215 * HC3Pg * NO3g : = MO2g ; % 0.053 * K215 * HC3Pg * NO3g : = ETHPg ; % 0.051 * K215 * HC3Pg * NO3g : = XO2g ; % K215 * HC3Pg * NO3g : = NO2g ; % 0.021 * K216 * HC5Pg * NO3g : = HCHOg ; % 0.239 * K216 * HC5Pg * NO3g : = ALDg ; % 0.828 * K216 * HC5Pg * NO3g : = KETg ; % 0.699 * K216 * HC5Pg * NO3g : = HO2g ; % 0.04 * K216 * HC5Pg * NO3g : = MO2g ; % 0.262 * K216 * HC5Pg * NO3g : = ETHPg ; % 0.391 * K216 * HC5Pg * NO3g : = XO2g ; % K216 : HC5Pg + NO3g = NO2g ; % 0.187 * K217 * HC8Pg * NO3g : = ALDg ; % 0.88 * K217 * HC8Pg * NO3g : = KETg ; % 0.845 * K217 * HC8Pg * NO3g : = HO2g ; % 0.155 * K217 * HC8Pg * NO3g : = ETHPg ; % 0.587 * K217 * HC8Pg * NO3g : = XO2g ; % K217 : HC8Pg + NO3g = NO2g ; % K218 : ETEPg + NO3g = NO2g + HO2g ; % 1.6 * K218 * ETEPg * NO3g : = HCHOg ; % 0.2 * K218 * ETEPg * NO3g : = ALDg ; % K219 : OLTPg + NO3g = NO2g + HO2g + HCHOg ; % 0.94 * K219 * OLTPg * NO3g : = ALDg ; % 0.06 * K219 * OLTPg * NO3g : = KETg ; % K220 : OLIPg + NO3g = NO2g + HO2g ; % 1.71 * K220 * OLIPg * NO3g : = ALDg ; % 0.29 * K220 * OLIPg * NO3g : = KETg ; % K221 : ISOPg + NO3g = NO2g + HO2g ; % 0.6 * K221 * ISOPg * NO3g : = MACRg ; % 0.4 * K221 * ISOPg * NO3g : = OLTg ; % 0.686 * K221 * ISOPg * NO3g : = HCHOg ; % K222 : APIPg + NO3g = HO2g + ALDg + NO2g+ KETg ; % K223 : LIMPg + NO3g = HO2g + NO2g ; % 0.6 * K223 * LIMPg * NO3g : = MACRg ; % 0.4 * K223 * LIMPg * NO3g : = OLIg + HCHOg ; % K224 : TOLPg + NO3g = NO2g + HO2g ; % 0.5 * K224 * TOLPg * NO3g : = DCBg ; % 1.3 * K224 * TOLPg * NO3g : = GLYg ; % 0.7 * K224 * TOLPg * NO3g : = MGLYg ; % K225 : XYLPg + NO3g = NO2g + HO2g + DCBg ; % 0.74 * K225 * XYLPg * NO3g : = GLYg ; % 1.26 * K225 * XYLPg * NO3g : = MGLYg ; % K226 : CSLPg + NO3g = GLYg + MGLYg + NO2g+ HO2g ; % K227 : ACO3g + NO3g = NO2g + MO2g ; % K228 : TCO3g + NO3g = NO2g + HCHOg + ACO3g ; % K229 : KETPg + NO3g = NO2g ; % 0.77 * K229 * KETPg * NO3g : = HO2g ; % 0.16 * K229 * KETPg * NO3g : = XO2g ; % 0.23 * K229 * KETPg * NO3g : = ACO3g ; % 0.46 * K229 * KETPg * NO3g : = ALDg ; % 0.54 * K229 * KETPg * NO3g : = MGLYg ; % K230 : OLNNg + NO3g = NO2g + HO2g + ONITg ; % 0.28 * K231 * OLNDg * NO3g : = HCHOg ; % 1.24 * K231 * OLNDg * NO3g : = ALDg ; % 0.469 * K231 * OLNDg * NO3g : = KETg ; % K231 : OLNDg + NO3g = NO2g + NO2g ; % K232 : XO2g + HO2g = OP2g ; % K233 : XO2g + MO2g = HO2g + HCHOg ; % K234 : XO2g + ACO3g = MO2g ; % K235 : XO2g + XO2g = ; % K236 : XO2g + NOg = NO2g ; % K237 : XO2g + NO3g = NO2g ; * ; * Aqueous phase chemical mechanism according to CAPRAM; * ; * ; KCAP1r % KCAP1 : H2O2a + Fepp = OHa + OHm + Feppp + KCAP1t ; KCAP2r % KCAP2 : H2O2a + Mnppp = HO2a + Hp + Mnpp + KCAP2t ; KCAP3r % KCAP3 : MnOHpp + H2O2a = MnO2p + Hp + KCAP3t ; KCAP4r % KCAP4 : Mnpppp + H2O2a = Mnpp + Hp + Hp + O2a + KCAP4t ; KCAP5r % KCAP5 : H2O2a + Cup = OHa + OHm + Cupp + KCAP5t ; KCAP6r % KCAP6 : O2m + Feppp = O2a + Fepp + KCAP6t ; KCAP7r % KCAP7 : HO2a + FeOHpp = Fepp + O2a + KCAP7t ; KCAP8r % KCAP8 : O2m + FeOHpp = O2a + Fepp + OHm + KCAP8t ; KCAP9r % KCAP9 : O2m + FeOH2p = O2a + Fepp + OHm + OHm + KCAP9t ; KCAP10r % KCAP10/Hp/Hp : O2m + Fepp + Hp + Hp = H2O2a + Feppp + KCAP10t ; KCAP11r % KCAP11/Hp : HO2a + Fepp + Hp= H2O2a + Feppp + KCAP11t ; KCAP12r % KCAP12 : OHa + Fepp = FeOHpp + KCAP12t ; KCAP13r % KCAP13/Hp : MnO2p + HO2a +Hp = Mnpp + H2O2a + O2a + KCAP13t ; KCAP14r % KCAP14 : OHa + Mnpp = OHm + Mnppp + KCAP14t ; KCAP15r % KCAP15/Hp/Hp : O2m + Cup + Hp+ Hp = H2O2a + Cupp + KCAP15t ; KCAP16r % KCAP16/Hp : HO2a + Cup+ Hp = H2O2a + Cupp + KCAP16t ; KCAP17r % KCAP17 : OHa + Cup = OHm + Cupp + KCAP17t ; KCAP18r % KCAP18 : HO2a + Cupp = O2a + Cup + Hp + KCAP18t ; KCAP19r % KCAP19 : O2m + Cupp = O2a + Cup + KCAP19t ; KCAP20r % KCAP20 : Feppp + Cup = Fepp + Cupp + KCAP20t ; KCAP21r % KCAP21 : FeOHpp + Cup = Fepp + Cupp + OHm + KCAP21t ; KCAP22r % KCAP22 : FeOH2p + Cup = Fepp + Cupp + OHm + OHm + KCAP22t ; KCAP23r % KCAP23 : Fepp + Mnppp = Feppp + Mnpp + KCAP23t ; KCAP24r % KCAP24 : O3a + O2m = O3m + O2a + KCAP24t ; KCAP25r % KCAP25 : HO3a = OHa + O2a + KCAP25t ; KCAP26r % KCAP26 : HO2a + HO2a = O2a + H2O2a + KCAP26t ; KCAP27r % KCAP27/Hp : HO2a + O2m+ Hp = H2O2a + O2a + KCAP27t ; KCAP28r % KCAP28 : HO2a + OHa = O2a + KCAP28t ; KCAP29r % KCAP29 : O2m + OHa = OHm + O2a + KCAP29t ; KCAP30r % KCAP30 : H2O2a + OHa = HO2a + KCAP30t ; KCAP31r % KCAP31 : CH3O2Ha + OHa = CH3O2a + KCAP31t ; KCAP32r % KCAP32 : CH3O2Ha + OHa = HCOOHa + HO2a + KCAP32t ; KCAP33r % KCAP33 : HSO3m + OHa = SO3m + KCAP33t ; KCAP34r % KCAP34 : SO3mm + OHa = OHm + SO3m + KCAP34t ; KCAP35r % KCAP35/Hp/Hp : MnO2p + MnO2p + Hp +Hp = Mnpp + Mnpp + H2O2a + KCAP35t ; KCAP36r % KCAP36 : Cup + O2a = Cupp + O2m + KCAP36t ; KCAP37r % KCAP37 : Fepp + O3a = FeOpp + O2a + KCAP37t ; KCAP38r % KCAP38 : FeOpp + H2O2a = Feppp + HO2a + OHm + KCAP38t ; KCAP39r % KCAP39 : FeOpp + H2Oa = Feppp + OHa + OHm + KCAP39t ; KCAP40r % KCAP40 : FeOpp + HO2a = Feppp + O2a + OHm + KCAP40t ; KCAP41r % KCAP41/Hp : FeOpp + OHa+Hp = Feppp + H2O2a + KCAP41t ; KCAP42r % KCAP42 : FeOpp + HNO2a = Feppp + NO2a + OHm + KCAP42t ; KCAP43r % KCAP43 : FeOpp + HSO3m = Feppp + SO3m + OHm + KCAP43t ; KCAP44r % KCAP44/Hp : FeOpp + Clm+ Hp = Feppp + ClOHm + KCAP44t ; KCAP45r % KCAP45/Hp/Hp : FeOpp + Mnpp+ Hp+ Hp = Mnppp + Feppp + KCAP45t ; KCAP46r % KCAP46/Hp : FeOpp + HCOOHa+ Hp = Feppp + CO2Ha + KCAP46t ; KCAP47r % KCAP47 : FeOpp + HCOOm = Feppp + CO2Ha + OHm + KCAP47t ; KCAP48r % KCAP48 : FeOpp + ETOHa = Feppp + CH3CHOHa + OHm + KCAP48t ; KCAP49r % KCAP49 : FeOpp + CH2OH2a = Feppp + OHm + CHOH2a + KCAP49t ; KCAP50r % KCAP50/Hp : FeOpp + NO2m+ Hp = Feppp + NO2a + OHm + KCAP50t ; KCAP51r % KCAP51 : FeOpp + Fepp = Feppp + Feppp + OHm + OHm + KCAP51t ; KCAP52r % KCAP52 : FeOpp + Fepp = HFe4p + KCAP52t ; KCAP53r % KCAP53/Hp : HFe4p + Hp + Hp = Feppp + Feppp + KCAP53t ; KCAP54r % KCAP54 : HFe4p = Feppp + Feppp + OHm + OHm + KCAP54t ; KCAP55r % KCAP55/Hp/Hp : MnOpp + Mnpp+ Hp+ Hp = Mnppp + Mnppp + KCAP55t ; KCAP56r % KCAP56 : Mnpp + O3a = MnOpp + O2a + KCAP56t ; KCAP57r % KCAP57/Hp : Cup + O3a+ Hp = Cupp + OHa + O2a + KCAP57t ; * KCAP58r % KCAP58 : N2O5a + H2Oa = Hp + Hp + NO3m + NO3m + KCAP58t ; KCAP59r % KCAP59 : N2O5a = NO2p + NO3m + KCAP59t ; KCAP60r % KCAP60 : NO2p + H2Oa = NO3m + Hp + Hp + KCAP60t ; KCAP61r % KCAP61 : NO3a + OHm = NO3m + OHa + KCAP61t ; KCAP62r % KCAP62 : NO3a + Fepp = NO3m + Feppp + KCAP62t ; KCAP63r % KCAP63 : NO3a + Mnpp = NO3m + Mnppp + KCAP63t ; KCAP64r % KCAP64 : NO3a + H2O2a = NO3m + Hp + HO2a + KCAP64t ; KCAP65r % KCAP65 : NO3a + CH3O2Ha = NO3m + Hp + CH3O2a + KCAP65t ; KCAP66r % KCAP66 : NO3a + HO2a = NO3m + Hp + O2a + KCAP66t ; KCAP67r % KCAP67 : NO3a + O2m = NO3m + O2a + KCAP67t ; KCAP68r % KCAP68 : NO3a + HSO3m = NO3m + Hp + SO3m + KCAP68t ; KCAP69r % KCAP69 : NO3a + SO3mm = NO3m + SO3m + KCAP69t ; KCAP70r % KCAP70 : NO3a + HSO4m = NO3m + Hp + SO4m + KCAP70t ; KCAP71r % KCAP71 : NO3a + SO4mm = NO3m + SO4m + KCAP71t ; KCAP72r % KCAP72 : NO2a + OHa = HOONOa + KCAP72t ; KCAP73r % KCAP73 : NO2a + O2m = NO2m + O2a + KCAP73t ; KCAP74r % KCAP74 : NO2a + NO2a = HNO2a + NO3m + Hp + KCAP74t ; KCAP75r % KCAP75 : O2NO2m = NO2m + O2a + KCAP75t ; KCAP76r % KCAP76 : NO2m + OHa = NO2a + OHm + KCAP76t ; KCAP77r % KCAP77 : NO2m + SO4m = SO4mm + NO2a + KCAP77t ; KCAP78r % KCAP78 : NO2m + NO3a = NO3m + NO2a + KCAP78t ; KCAP79r % KCAP79 : NO2m + Cl2m = Clm + Clm + NO2a + KCAP79t ; KCAP80r % KCAP80 : NO2m + Br2m = Brm + Brm + NO2a + KCAP80t ; KCAP81r % KCAP81 : NO2m + CO3m = CO3mm + NO2a + KCAP81t ; KCAP82r % KCAP82 : NO2m + O3a = NO3m + O2a + KCAP82t ; KCAP83r % KCAP83 : HNO2a + OHa = NO2a + KCAP83t ; *NO3m pour NO3a cidessous; KCAP84r % KCAP84 : HNO4a + HSO3m = HSO4m + NO3m + Hp + KCAP84t ; KCAP85r % KCAP85 : NO2p + Clm = ClNO2a + KCAP85t ; KCAP86r % KCAP86 : NO2p + Brm = BrNO2a + KCAP86t ; KCAP87r % KCAP87 : ClNO2a + Brm = NO2m + BrCla + KCAP87t ; KCAP88r % KCAP88 : BrNO2a + Brm = NO2m + Br2a + KCAP88t ; KCAP89r % KCAP89 : BrNO2a + Clm = NO2m + BrCla + KCAP89t ; KCAP90r % KCAP90 : HMSm + OHa = CHOHSO3m + KCAP90t ; KCAP91r % KCAP91 : HMSm + SO4m = SO4mm + OHCH2SO3a + KCAP91t ; KCAP92r % KCAP92 : HMSm + NO3a = NO3m + OHCH2SO3a + KCAP92t ; KCAP93r % KCAP93 : HMSm + Cl2m = Clm + Clm + OHCH2SO3a + KCAP93t ; KCAP94r % KCAP94 : HMSm + Br2m = Brm + Brm + OHCH2SO3a + KCAP94t ; KCAP95r % KCAP95 : OHCH2SO3a = SO3a + CH2OHa + KCAP95t ; KCAP96r % KCAP96 : SO3a + H2Oa = SO4mm + Hp + Hp + KCAP96t ; KCAP97r % KCAP97 : CHOHSO3m + O2a = HS2m + KCAP97t ; KCAP98r % KCAP98 : HS2m = HO2a + CHOSO3m + KCAP98t ; KCAP99r % KCAP99 : HS2m = O2CHOa + HSO3m + KCAP99t ; KCAP100r % KCAP100 : O2CHOa + H2Oa = HCOOHa + HO2a + KCAP100t ; KCAP101r % KCAP101 : CHOSO3mm + O2a = CHOSO3m + O2m + KCAP101t ; KCAP102r % KCAP102 : HSO3m + H2O2a + Hp = SO4mm + Hp + Hp + KCAP102t ; KCAP103r % KCAP103 : HSO3m + CH3O2Ha + Hp = SO4mm + Hp + Hp + CH3OHa + KCAP103t ; KCAP104r % KCAP104 : HSO3m + PAAa + Hp = SO4mm + Hp + Hp + CH3COOHa + KCAP104t ; KCAP105r % KCAP105 : SO2a + O3a = HSO4m + O2a + Hp + KCAP105t ; KCAP106r % KCAP106 : HSO3m + O3a = SO4mm + Hp + O2a + KCAP106t ; KCAP107r % KCAP107 : SO3mm + O3a = SO4mm + O2a + KCAP107t ; KCAP108r % KCAP108 : FeOHpp + HSO3m = Fepp + SO3m + KCAP108t ; KCAP109r % KCAP109 : Fepp + SO5m = FeOHpp + HSO5m + KCAP109t ; KCAP110r % KCAP110 : Fepp + HSO5m = FeOHpp + SO4m + KCAP110t ; KCAP111r % KCAP111 : Mnpp + SO5m = Mnppp + HSO5m + OHm + KCAP111t ; * KCAP112r % KCAP112 : Fepp + SO4m = FeOHpp + SO4mm + KCAP112t ; KCAP112r % KCAP112 : Fepp + SO4m = Feppp + SO4mm + Hp + KCAP112t ; KCAP113r % KCAP113 : Fepp + S2O8mm = Feppp + SO4mm + SO4m + KCAP113t ; KCAP114r % KCAP114 : SO5m + SO5m = S2O8mm + O2a + KCAP114t ; KCAP115r % KCAP115 : SO5m + SO5m = SO4m + SO4m + O2a + KCAP115t ; KCAP116r % KCAP116 : SO5m + HO2a = SO5O2Hm + KCAP116t ; KCAP117r % KCAP117/Hp : SO5m + O2m+ Hp = HSO5m + O2a + KCAP117t ; KCAP118r % KCAP118 : SO5O2Hm + HSO3m = SO4mm + SO4mm + Hp + Hp + KCAP118t ; KCAP119r % KCAP119 : SO5O2mm + HSO3m = SO4mm + SO4mm + Hp + KCAP119t ; KCAP120r % KCAP120/Hp : SO5O2mm+ Hp = HSO5m + O2a + KCAP120t ; KCAP121r % KCAP121 : SO3m + O2a = SO5m + KCAP121t ; KCAP122r % KCAP122 : SO5m + HSO3m = HSO5m + SO3m + KCAP122t ; KCAP123r % KCAP123 : SO5m + HSO3m = SO4mm + SO4m + Hp + KCAP123t ; KCAP124r % KCAP124/Hp : SO5m + SO3mm+ Hp = HSO5m + SO3m + KCAP124t ; KCAP125r % KCAP125 : SO5m + SO3mm = SO4m + SO4mm + KCAP125t ; KCAP126r % KCAP126 : OHa + HSO4m = SO4m + KCAP126t ; KCAP127r % KCAP127 : SO4m + SO4m = S2O8mm + KCAP127t ; KCAP128r % KCAP128 : SO4m + HSO3m = SO4mm + SO3m + Hp + KCAP128t ; KCAP129r % KCAP129 : SO4m + SO3mm = SO4mm + SO3m + KCAP129t ; KCAP130r % KCAP130 : SO4m + Mnpp = SO4mm + Mnppp + KCAP130t ; KCAP131r % KCAP131 : SO4m + Cup = SO4mm + Cupp + KCAP131t ; KCAP132r % KCAP132 : SO4m + H2O2a = SO4mm + Hp + HO2a + KCAP132t ; KCAP133r % KCAP133 : SO4m + CH3O2Ha = SO4mm + Hp + CH3O2a + KCAP133t ; KCAP134r % KCAP134 : SO4m + HO2a = SO4mm + Hp + O2a + KCAP134t ; KCAP135r % KCAP135 : SO4m + O2m = SO4mm + O2a + KCAP135t ; KCAP136r % KCAP136 : SO4m + NO3m = SO4mm + NO3a + KCAP136t ; KCAP137r % KCAP137 : SO4m + OHm = SO4mm + OHa + KCAP137t ; KCAP138r % KCAP138 : SO4m + H2Oa = SO4mm + Hp + OHa + KCAP138t ; KCAP139r % KCAP139 : HSO5m + HSO3m + Hp = SO4mm + SO4mm + Hp + Hp + Hp + KCAP139t ; KCAP140r % KCAP140 : HSO5m + SO3mm + Hp = SO4mm + SO4mm + Hp + Hp + KCAP140t ; KCAP141r % KCAP141 : HSO5m + OHa = SO5m + KCAP141t ; KCAP142r % KCAP142 : OHa + SO4m = HSO5m + KCAP142t ; KCAP143r % KCAP143 : MnHSO3p + Mnppp = Mnpp + Mnpp + SO3m + Hp + KCAP143t ; KCAP144r % KCAP144 : CH3OHa + OHa = CH2OHa + KCAP144t ; KCAP145r % KCAP145 : CH3OHa + SO4m = SO4mm + Hp + CH2OHa + KCAP145t ; KCAP146r % KCAP146 : CH3OHa + NO3a = NO3m + Hp + CH2OHa + KCAP146t ; KCAP147r % KCAP147 : CH3OHa + Cl2m = Clm + Clm + Hp + CH2OHa + KCAP147t ; KCAP148r % KCAP148 : CH3OHa + Br2m = Brm + Brm + Hp + CH2OHa + KCAP148t ; KCAP149r % KCAP149 : CH3OHa + CO3m = CO3mm + Hp + CH2OHa + KCAP149t ; KCAP150r % KCAP150 : CH2OHa + O2a = O2CH2OHa + KCAP150t ; KCAP151r % KCAP151 : O2CH2OHa + OHm = HCHOa + O2m + KCAP151t ; KCAP152r % KCAP152 : O2CH2OHa + O2CH2OHa = CH3OHa + O2a + HCHOa + KCAP152t ; KCAP153r % KCAP153 : ETOHa + OHa = CH3CHOHa + KCAP153t ; KCAP154r % KCAP154 : ETOHa + SO4m = SO4mm + Hp + CH3CHOHa + KCAP154t ; KCAP155r % KCAP155 : ETOHa + NO3a = NO3m + Hp + CH3CHOHa + KCAP155t ; KCAP156r % KCAP156 : ETOHa + Cl2m = Clm + Clm + Hp + CH3CHOHa + KCAP156t ; KCAP157r % KCAP157 : ETOHa + Br2m = Brm + Brm + Hp + CH3CHOHa + KCAP157t ; KCAP158r % KCAP158 : ETOHa + CO3m = CO3mm + Hp + CH3CHOHa + KCAP158t ; KCAP159r % KCAP159 : CH3CHOHa + O2a = HOP3a + KCAP159t ; KCAP160r % KCAP160 : HOP3a = CH3CHOa + HO2a + KCAP160t ; KCAP161r % KCAP161 : HOP3a + OHm = CH3CHOa + O2m + KCAP161t ; KCAP162r % KCAP162 : CH2OH2a + OHa = CHOH2a + KCAP162t ; KCAP163r % KCAP163 : CH2OH2a + SO4m = SO4mm + Hp + CHOH2a + KCAP163t ; KCAP164r % KCAP164 : CH2OH2a + NO3a = NO3m + Hp + CHOH2a + KCAP164t ; KCAP165r % KCAP165 : CH2OH2a + Cl2m = Clm + Clm + Hp + CHOH2a + KCAP165t ; KCAP166r % KCAP166 : CH2OH2a + Br2m = Brm + Brm + Hp + CHOH2a + KCAP166t ; KCAP167r % KCAP167 : CH2OH2a + CO3m = CO3mm + Hp + CHOH2a + KCAP167t ; KCAP168r % KCAP168 : CHOH2a + O2a = HO2a + HCOOHa + KCAP168t ; KCAP169r % KCAP169 : ALCOOL1a + OHa = CH3COH2a + KCAP169t ; KCAP170r % KCAP170 : CH3CHOa + OHa = CH3COa + KCAP170t ; KCAP171r % KCAP171 : CH3CHOa + SO4m = SO4mm + Hp + CH3COa + KCAP171t ; KCAP172r % KCAP172 : ALCOOL1a + SO4m = SO4mm + Hp + CH3COH2a + KCAP172t ; KCAP173r % KCAP173 : ALCOOL1a + NO3a = NO3m + Hp + CH3COH2a + KCAP173t ; KCAP174r % KCAP174 : CH3CHOa + NO3a = NO3m + Hp + CH3COa + KCAP174t ; KCAP175r % KCAP175 : CH3CHOa + Cl2m = Clm + Clm + Hp + CH3COa + KCAP175t ; KCAP176r % KCAP176 : ALCOOL1a + Cl2m = Clm + Clm + Hp + CH3COH2a + KCAP176t ; KCAP177r % KCAP177 : CH3CHOa + Br2m = Brm + Brm + Hp + CH3COa + KCAP177t ; KCAP178r % KCAP178 : ALCOOL1a + Br2m = Brm + Brm + Hp + CH3COH2a + KCAP178t ; KCAP179r % KCAP179 : CH3CHOa + CO3m = CO3mm + Hp + CH3COa + KCAP179t ; KCAP180r % KCAP180 : ALCOOL1a + CO3m = CO3mm + Hp + CH3COH2a + KCAP180t ; KCAP181r % KCAP181 : HCOOHa + OHa = CO2m + Hp + KCAP181t ; KCAP182r % KCAP182 : HCOOm + OHa = OHm + CO2m + Hp + KCAP182t ; KCAP183r % KCAP183 : HCOOHa + SO4m = SO4mm + Hp + CO2Ha + KCAP183t ; KCAP184r % KCAP184 : HCOOm + SO4m = SO4mm + CO2Ha + KCAP184t ; KCAP185r % KCAP185 : HCOOHa + NO3a = NO3m + Hp + CO2Ha + KCAP185t ; KCAP186r % KCAP186 : HCOOm + NO3a = NO3m + CO2Ha + KCAP186t ; KCAP187r % KCAP187 : HCOOHa + Cl2m = Clm + Clm + Hp + CO2Ha + KCAP187t ; KCAP188r % KCAP188 : HCOOm + Cl2m = Clm + Clm + CO2Ha + KCAP188t ; KCAP189r % KCAP189 : HCOOHa + Br2m = Brm + Brm + Hp + CO2Ha + KCAP189t ; KCAP190r % KCAP190 : HCOOm + Br2m = Brm + Brm + CO2Ha + KCAP190t ; KCAP191r % KCAP191 : HCOOm + CO3m = CO3mm + CO2Ha + KCAP191t ; KCAP192r % KCAP192 : CO2Ha + O2a = CO2a + HO2a + KCAP192t ; KCAP193r % KCAP193 : CH3COOHa + OHa = CH2COOHa + KCAP193t ; KCAP194r % KCAP194 : CH3COOm + OHa = CH2COOm + KCAP194t ; KCAP195r % KCAP195 : CH3COOHa + SO4m = SO4mm + Hp + CH2COOHa + KCAP195t ; KCAP196r % KCAP196 : CH3COOm + SO4m = SO4mm + CH3a + CO2a + KCAP196t ; KCAP197r % KCAP197 : CH3COOHa + NO3a = NO3m + Hp + CH2COOHa + KCAP197t ; KCAP198r % KCAP198 : CH3COOm + NO3a = NO3m + CH3a + CO2a + KCAP198t ; KCAP199r % KCAP199 : CH3COOHa + Cl2m = Clm + Clm + Hp + CH2COOHa + KCAP199t ; KCAP200r % KCAP200 : CH3COOm + Cl2m = Clm + Clm + CH3a + CO2a + KCAP200t ; KCAP201r % KCAP201 : CH3COOHa + Br2m = Brm + Brm + Hp + CH2COOHa + KCAP201t ; KCAP202r % KCAP202 : CH2COOHa + O2a = ACO3a + KCAP202t ; KCAP203r % KCAP203 : CH3COOm + Br2m = Brm + Brm + CH3a + CO2a + KCAP203t ; KCAP204r % KCAP204 : CH3COOm + CO3m = CO3mm + CH3a + CO2a + KCAP204t ; KCAP205r % KCAP205 : CH3a + O2a = CH3O2a + KCAP205t ; KCAP206r % KCAP206 : CH3O2a + CH3O2a = CH3OHa + HCHOa + O2a + KCAP206t ; KCAP207r % KCAP207 : CH3O2a + CH3O2a = CH3Oa + CH3Oa + O2a + KCAP207t ; KCAP208r % KCAP208 : ACO3a + ACO3a = CH3O2a + CH3O2a + CO2a + CO2a + O2a + KCAP208t ; KCAP209r % KCAP209 : CH3O2a + HSO3m = CH3O2Ha + SO3m + KCAP209t ; KCAP210r % KCAP210 : ETHPa + ETHPa = ETOHa + CH3CHOa + O2a + KCAP210t ; KCAP211r % KCAP211 : ETHPa + ETHPa = CH3CH2Oa + CH3CH2Oa + O2a + KCAP211t ; KCAP212r % KCAP212 : CH3CH2Oa + O2a = HO2a + CH3CHOa + KCAP212t ; KCAP213r % KCAP213 : CH3CH2Oa = CH3CHOHa + KCAP213t ; KCAP214r % KCAP214 : Fepp + CH3O2a = FeMO2pp + KCAP214t ; KCAP215r % KCAP215 : FeMO2pp + Hp = Feppp + CH3O2Ha + KCAP215t ; KCAP216r % KCAP216 : FeMO2pp = Feppp + CH3O2Ha + OHm + KCAP216t ; KCAP217r % KCAP217 : OHa + HC2O4m = C2O4m + KCAP217t ; KCAP218r % KCAP218 : NO3a + HC2O4m = NO3m + Hp + C2O4m + KCAP218t ; KCAP219r % KCAP219 : Cl2m + HC2O4m = Clm + Clm + Hp + C2O4m + KCAP219t ; KCAP220r % KCAP220 : Br2m + HC2O4m = Brm + Brm + Hp + C2O4m + KCAP220t ; KCAP221r % KCAP221 : SO4m + HC2O4m = SO4mm + Hp + C2O4m + KCAP221t ; KCAP222r % KCAP222 : OHa + C2O4mm = OHm + C2O4m + KCAP222t ; KCAP223r % KCAP223 : NO3a + C2O4mm = NO3m + C2O4m + KCAP223t ; KCAP224r % KCAP224 : Cl2m + C2O4mm = Clm + Clm + C2O4m + KCAP224t ; KCAP225r % KCAP225 : Br2m + C2O4mm = Brm + Brm + C2O4m + KCAP225t ; KCAP226r % KCAP226 : SO4m + C2O4mm = SO4mm + C2O4m + KCAP226t ; KCAP227r % KCAP227 : C2O4m + O2a = CO2a + CO2a + O2m + KCAP227t ; KCAP228r % KCAP228 : OHa + CHOH22a = COH22a + KCAP228t ; KCAP229r % KCAP229 : NO3a + CHOH22a = Hp + NO3m + COH22a + KCAP229t ; KCAP230r % KCAP230 : Cl2m + CHOH22a = Hp + Clm + Clm + COH22a + KCAP230t ; KCAP231r % KCAP231 : Br2m + CHOH22a = Hp + Brm + Brm + COH22a + KCAP231t ; KCAP232r % KCAP232 : SO4m + CHOH22a = Hp + SO4mm + COH22a + KCAP232t ; KCAP233r % KCAP233 : COH22a + O2a = HOP2a + KCAP233t ; KCAP234r % KCAP234 : HOP2a = HO2a + HOA2a + KCAP234t ; KCAP235r % KCAP235 : OHa + HOA2a = HOA1a + KCAP235t ; KCAP236r % KCAP236 : NO3a + HOA2a = Hp + NO3m + HOA1a + KCAP236t ; KCAP237r % KCAP237 : Cl2m + HOA2a = Hp + Clm + Clm + HOA1a + KCAP237t ; KCAP238r % KCAP238 : Br2m + HOA2a = Hp + Brm + Brm + HOA1a + KCAP238t ; KCAP239r % KCAP239 : SO4m + HOA2a = Hp + SO4mm + HOA1a + KCAP239t ; KCAP240r % KCAP240 : HOA1a + O2a = HOA3a + KCAP240t ; KCAP241r % KCAP241 : HOA3a = HO2a + Hp + HC2O4m + KCAP241t ; KCAP242r % KCAP242 : CH3COa + O2a = ACO3a + KCAP242t ; KCAP243r % KCAP243 : ACO3a + O2m = CH3COOOm + KCAP243t ; KCAP244r % KCAP244 : CH3COH2a + O2a = HOP1a + KCAP244t ; KCAP245r % KCAP245 : HOP1a = Hp + Hp + CH3COOm + O2m + KCAP245t ; KCAP246r % KCAP246 : CH3Oa + O2a = HCHOa + HO2a + KCAP246t ; KCAP247r % KCAP247 : CH3Oa = CH2OHa + KCAP247t ; KCAP248r % KCAP248 : HOA4m + HOA4m = HOA2m + HOA2m + H2O2a + KCAP248t ; KCAP249r % KCAP249 : HOA4m + HOA4m = HCHOa + HCHOa + H2O2a+ OHm + OHm + CO2a + CO2a + KCAP249t ; KCAP250r % KCAP250 : HOA4m + HOA4m = HOA2m + HOA3m + O2a + KCAP250t ; KCAP251r % KCAP251 : HOA4m + HOA4m = O2m + O2m + HOA2m + KCAP251t ; KCAP252r % KCAP252 : Cl2m + Cl2m = Cl2a + Clm + Clm + KCAP252t ; KCAP253r % KCAP253 : Cl2m + Fepp = Clm + Clm + Feppp + KCAP253t ; KCAP254r % KCAP254 : Cl2m + Fepp = FeClpp + Clm + KCAP254t ; KCAP255r % KCAP255 : Cl2m + Mnpp = MnCl2p + KCAP255t ; KCAP256r % KCAP256 : MnCl2p = Mnpp + Cl2m + KCAP256t ; KCAP257r % KCAP257 : MnCl2p = Mnppp + Clm + Clm + KCAP257t ; KCAP258r % KCAP258 : Cl2m + Cup = Clm + Clm + Cupp + KCAP258t ; KCAP259r % KCAP259 : Cl2m + H2O2a = Clm + Clm + Hp + HO2a + KCAP259t ; KCAP260r % KCAP260 : Cl2m + CH3O2Ha = Clm + Clm + Hp + CH3O2a + KCAP260t ; KCAP261r % KCAP261 : Cl2m + OHm = Clm + Clm + OHa + KCAP261t ; KCAP262r % KCAP262 : Cl2m + HO2a = Clm + Clm + Hp + O2a + KCAP262t ; KCAP263r % KCAP263 : Cl2m + O2m = Clm + Clm + O2a + KCAP263t ; KCAP264r % KCAP264 : Cl2m + HSO3m = Clm + Clm + Hp + SO3m + KCAP264t ; KCAP265r % KCAP265 : Cl2m + SO3mm = Clm + Clm + SO3m + KCAP265t ; KCAP266r % KCAP266 : Cl2a + H2Oa = Hp + Clm + HOCla + KCAP266t ; KCAP267r % KCAP267 : Cl2m + H2Oa = Hp + Clm + Clm + OHa + KCAP267t ; KCAP268r % KCAP268 : SO4m + Brm = SO4mm + Bra + KCAP268t ; KCAP269r % KCAP269 : NO3a + Brm = NO3m + Bra + KCAP269t ; KCAP270r % KCAP270 : Br2m + Br2m = Br2a + Brm + Brm + KCAP270t ; KCAP271r % KCAP271 : Br2m + Fepp = Brm + Brm + Feppp + KCAP271t ; KCAP272r % KCAP272 : Br2m + Mnpp = MnBr2p + KCAP272t ; KCAP273r % KCAP273 : MnBr2p = Br2m + Mnpp + KCAP273t ; KCAP274r % KCAP274 : MnBr2p = Brm + Brm + Mnppp + KCAP274t ; KCAP275r % KCAP275 : Br2m + Cup = Brm + Brm + Cupp + KCAP275t ; KCAP276r % KCAP276 : Br2m + H2O2a = Brm + Brm + Hp + HO2a + KCAP276t ; KCAP277r % KCAP277 : Br2m + CH3O2Ha = Brm + Brm + Hp + CH3O2a + KCAP277t ; KCAP278r % KCAP278 : Br2m + OHm = Brm + Brm + OHa + KCAP278t ; KCAP279r % KCAP279 : Br2m + HO2a = Brm + Brm + Hp + O2a + KCAP279t ; KCAP280r % KCAP280 : Br2m + O2m = Brm + Brm + O2a + KCAP280t ; KCAP281r % KCAP281 : Br2m + HSO3m = Brm + Brm + Hp + SO3m + KCAP281t ; KCAP282r % KCAP282 : Br2m + SO3mm = Brm + Brm + SO3m + KCAP282t ; KCAP283r % KCAP283 : Br2a + H2Oa = Brm + Hp + HOBra + KCAP283t ; KCAP284r % KCAP284 : BrOHm = Bra + OHm + KCAP284t ; KCAP285r % KCAP285 : HCO3m + OHa = CO3m + KCAP285t ; KCAP286r % KCAP286 : CO3mm + OHa = OHm + CO3m + KCAP286t ; KCAP287r % KCAP287 : CO3mm + SO4m = SO4mm + CO3m + KCAP287t ; KCAP288r % KCAP288 : HCO3m + SO4m = SO4mm + CO3m + Hp + KCAP288t ; KCAP289r % KCAP289 : CO3mm + NO3a = NO3m + CO3m + KCAP289t ; KCAP290r % KCAP290 : HCO3m + NO3a = NO3m + CO3m + Hp + KCAP290t ; KCAP291r % KCAP291 : CO3mm + Cl2m = Clm + Clm + CO3m + KCAP291t ; KCAP292r % KCAP292 : HCO3m + Cl2m = Clm + Clm + CO3m + Hp + KCAP292t ; KCAP293r % KCAP293 : HCO3m + Br2m = Brm + Brm + CO3m + Hp + KCAP293t ; KCAP294r % KCAP294 : CO3m + CO3m = O2m + O2m + CO2a + CO2a + KCAP294t ; KCAP295r % KCAP295 : CO3m + Fepp = CO3mm + Feppp + KCAP295t ; KCAP296r % KCAP296 : CO3m + Mnpp = CO3mm + Mnppp + KCAP296t ; KCAP297r % KCAP297 : CO3m + Cup = CO3mm + Cupp + KCAP297t ; KCAP298r % KCAP298 : CO3m + H2O2a = HCO3m + HO2a + KCAP298t ; KCAP299r % KCAP299 : CO3m + CH3O2Ha = HCO3m + CH3O2a + KCAP299t ; KCAP300r % KCAP300 : CO3m + HO2a = HCO3m + O2a + KCAP300t ; KCAP301r % KCAP301 : CO3m + O2m = CO3mm + O2a + KCAP301t ; KCAP302r % KCAP302 : CO3m + HSO3m = HCO3m + SO3m + KCAP302t ; KCAP303r % KCAP303 : CO3m + SO3mm = CO3mm + SO3m + KCAP303t ; KCAP304r % KCAP304 : CO3m + NO2a = CO2a + NO3m + KCAP304t ; KCAP305r % KCAP305 : CO3m + O3a = CO2a + O2a + O2m + KCAP305t ; KCAP306r % KCAP306 : SO5m + HC2O4m = HSO5m + C2O4m + KCAP306t ; KCAP307r % KCAP307/Hp : SO5m + C2O4mm + Hp = HSO5m + C2O4m + KCAP307t ; KCAP308r % KCAP308 : SO5m + CHOH22a = HSO5m + COH22a + KCAP308t ; KCAP309r % KCAP309 : CH2COOm + O2a = HOA4m + KCAP309t ; KCAP310r % KCAP310 : OHa + O3a = O2a + HO2a + KCAP310t ; KCAP312r % KCAP312 : Mnpp + HSO5m = Mnppp + SO4m + OHm + KCAP312t ; KCAP313r % KCAP313 : CHOSO3m + H2Oa = HSO3m + HCOOHa + KCAP313t ; % E1FWD*H2Oa % E1BCK : = Hp + OHm ; % E2FWD % E2BCK : CO2a + H2Oa = H2CO3a ; % E3FWD % E3BCK : H2CO3a = Hp + HCO3m ; % E4FWD % E4BCK : HCO3m = Hp + CO3mm ; % E5FWD % E5BCK : HCla = Hp + Clm ; % E6FWD % E6BCK : NH3a + H2Oa = NH4p + OHm ; % E7FWD % E7BCK : HO2a = Hp + O2m ; % E8FWD % E8BCK : HNO3a = Hp + NO3m ; % E9FWD % E9BCK : HNO2a = Hp + NO2m ; % E10FWD % E10BCK : HNO4a = Hp + O2NO2m ; % E11FWD % E11BCK : NO2a + HO2a = HNO4a ; % E12FWD % E12BCK : SO2a + H2Oa = HSO3m + Hp ; % E13FWD % E13BCK : HSO3m = SO3mm + Hp ; % E14FWD % E14BCK : H2SO4a = HSO4m + Hp ; % E15FWD % E15BCK : HSO4m = SO4mm + Hp ; % E16FWD % E16BCK : HCOOHa = HCOOm + Hp ; % E17FWD % E17BCK : CH3COOHa = CH3COOm + Hp ; % E18FWD % E18BCK : Feppp + H2Oa = FeOHpp + Hp ; % E19FWD % E19BCK : FeOHpp + H2Oa = FeOH2p + Hp ; % E20FWD % E20BCK : Feppp + SO4mm = FeSO4p ; % E21FWD % E21BCK : HCHOa + H2Oa = CH2OH2a ; % E22FWD % E22BCK : CH3CHOa + H2Oa = ALCOOL1a ; % E23FWD % E23BCK : CH2OH2a + HSO3m = HMSm + H2Oa ; % E24FWD % E24BCK : CH2OH2a + SO3mm = HMSm + OHm ; % E25FWD % E25BCK : Cla + Clm = Cl2m ; % E26FWD % E26BCK : Bra + Brm = Br2m ; % E27FWD % E27BCK : Clm + OHa = ClOHm ; % E28FWD % E28BCK : ClOHm + Hp = H2Oa + Cla ; % E29FWD % E29BCK : ClOHm + Clm = Cl2m + OHm ; % E30FWD % E30BCK : Brm + OHa = BrOHm ; % E31FWD % E31BCK : BrOHm + Hp = Bra + H2Oa ; % E32FWD % E32BCK : BrOHm + Brm = Br2m + OHm ; % E33FWD % E33BCK : Mnppp + H2Oa = MnOHpp + Hp ; % E34FWD % E34BCK : O2m + Mnpp = MnO2p ; % E35FWD % E35BCK : HO2a + Mnpp = MnO2p + Hp ; % E36FWD % E36BCK : Mnppp + Mnppp =Mnpp + Mnpppp ; % E37FWD % E37BCK : MnOHpp + H2Oa = MnOH2p + Hp ; % E38FWD % E38BCK : Mnpp + HSO3m = MnHSO3p ; % E39FWD % E39BCK : Cupp + OHa =CuOHpp ; % E40FWD % E40BCK : HO3a = Hp + O3m ; % E41FWD % E41BCK : HOONOa = Hp + OONOm ; % E42FWD % E42BCK : CHOHSO3m = CHOSO3mm + Hp ; % E43FWD % E43BCK : SO5O2Hm= SO5O2mm + Hp ; % E44FWD % E44BCK : H2C2O4a = Hp + HC2O4m ; % E45FWD % E45BCK : HC2O4m= Hp + C2O4mm ; % E46FWD % E46BCK : HOA2a =Hp + HOA2m ; % E47FWD % E47BCK : GLYa + H2Oa = CHOH22a ; % E48FWD % E48BCK : FeC2O4p = Feppp + C2O4mm ; % E49FWD % E49BCK : FeC2O42m = FeC2O4p + C2O4mm ; % E50FWD % E50BCK : FeORG3 = FeC2O42m + C2O4mm ; % E51FWD % E51BCK : SO4m + Clm = SO4mm + Cla ; % E52FWD % E52BCK : NO3a + Clm = NO3m + Cla ; % E53FWD % E53BCK : Clm + Feppp = FeClpp ; % E54FWD % E54BCK : PAAa = CH3COOOm + Hp ; % E55FWD % E55BCK : CH3COa + H2Oa = CH3COH2a ; % E56FWD % E56BCK : Brm + Brm + CO3m = Br2m + CO3mm ; % E57FWD % E57BCK : ACO3a =HOA4m+ Hp ; * **% E58FWD % E58BCK : ETHPa = HOA3m + Hp ; % P1 : H2O2a = OHa + OHa ; % P2 : FeOHpp = Fepp + OHa ; % P3 : FeOH2p = Fepp + OHa + OHm ; % P4 : FeSO4p = Fepp + SO4m ; *% P5 : NO2m = NOa + OHa + OHm ; *% P6 : NO3m = NO2a + OHa + OHm ; % P5 : NO2m = NOa + OHa ; % P6 : NO3m = NO2a + OHa ; % P7 : FeC2O42m = Fepp + C2O4mm + CO2a + CO2m ; % P8 : FeORG3 = Fepp + C2O4mm + C2O4mm + CO2a + CO2m ; % 4E9*1000 : CO2m + O2a = CO2a + O2m ; % P9 : CH3O2Ha = CH3Oa + OHa ; P10r % P10 : NO3a = NOa + O2a + P10t; P11r % P11 : NO3a = NO2a + O3Pa + P11t; % 4E9*1000 : O3Pa + O2a = O3a ; * ; * ******************************************************************; * Depositions; * ******************************************************************; * ; % VDCO : COg = ; % VDHOBr : HOBrg = ; % VDHOCl : HOClg = ; % VDNO2 : NO2g = ; % VDO3 : O3g = ; % VDETOH : ETOHg = ; % VDOP1 : OP1g = ; % VDCH3OH : CH3OHg = ; % VDH2O2 : H2O2g = ; % VDHBr : HBrg = ; % VDHCHO : HCHOg = ; % VDHCl : HClg = ; % VDNH3 : NH3g = ; % VDORA1 : ORA1g = ; % VDSO2 : SO2g = ; % VDHNO3 : HNO3g = ; % VDN2O5 : N2O5g = ; % VDSULF : SULFg = ; * ; * ****************************************************************; * Duration ; * ****************************************************************; * ; % 1/3600 : = hours ; **; * *******************************************************************; * Definition of output; * *******************************************************************; SETPSTREAM 1 1 20; Hours O3g H2O2g NOg NO2g NO3g N2O5g HNO2g HNO3g HNO4g SO2g SULFg OHg HO2g ETHg HC3g HC5g HC8g ETEg OLTg ; **; SETPSTREAM 2 2 20; ISOg APIg LIMg TOLg XYLg CSLg HCHOg ALDg KETg GLYg MGLYg DCBg MACRg UDDg HKETg ONITg PANg TPANg OP1g OP2g ; **; SETPSTREAM 3 3 20; ACO3g TCO3g KETPg OLNNg OLNDg XO2g HC3Pg HC5Pg HC8Pg KHETg XOg CH3OHg OL2Pg ETOHg SO2g ORA1g ORA2g; **; SETPSTREAM 4 4 20; * hours Cup Cupp SO2a HSO3m NO3a Cla NO3m clm Hp OHm Fepp Feppp GLYa H2O2a OHa HO2a H2SO4a HSO4m SO4m ; hours NO3a NO3g OHa OHg O2a Clm SO4m NO3m ON HET SZA P1 P10; **; COMPILE PRINT; PSTREAM 1; PSTREAM 2; PSTREAM 3; PSTREAM 4; **; WHENEVER TIME = 0+2160*160 % CALL PRINT ; TIME = 14400+86400 CALL DAY RESTART; TIME = 72000+86400 CALL NIGHT RESTART; TIME = 1 CALL HETERO RESTART; **; * MONITOR 3 ; * KMAX 3 ; * HMAX 3; * error 1E-4; BEGIN; STOP; * *********************************************************************; * END *; * *********************************************************************;