########################################################
# Melchior Mechanism: Adjusted from Hov et al (1985)   #
# by M. Lattuati, M. Beekmann (Service d'Aeronomie,    #
# CNRS, Paris, France) (Lattuati, 1997), simplified    #
# and updated by M. Beekmann and R. Vautard            #
# Last modification: February 1998                     #
########################################################

#############################
#    Inorganic chemistry    #
#############################

O3+NO->NO2                k(T)=Aexp(-B/T),A=1.8e-12,B=1370
O3+NO2->NO3               k(T)=Aexp(-B/T),A=1.2e-13,B=2450
O3+OH->HO2                k(T)=Aexp(-B/T),A=1.9e-12,B=1000
O3+HO2->OH                k(T)=Aexp(-B/T),A=1.4e-14,B=600
NO+HO2->OH+NO2            k(T)=Aexp(-B/T),A=3.7e-12,B=-240
NO2+OH+M->HNO3            k(T,M)=mtroe(3.4e-30,0,3.2,4.77e-11,0,1.4,0.30)
HO2+OH->H2O               k(T)=Aexp(-B/T),A=4.8e-11,B=-250
H2O2+OH->HO2              k(T)=Aexp(-B/T),A=2.9e-12,B=160
HNO3+OH->NO3              k(T)=Aexp(-B/T),A=7.2e-15,B=-785
HNO3+OH+M->NO3            k(T,M)=troe(1.9e-33,-725,0,4.1e-16,-1440,0,1)
CO+OH->HO2+CO2            k(T)=Aexp(-B/T)(300/T)**N,A=1.3e-13,B=0,N=1
CO+OH+M->HO2+CO2          k=3.21e-33
HO2+HO2->H2O2             k(T)=Aexp(-B/T),A=2.2e-13,B=-600
HO2+HO2+M->H2O2           k(T)=Aexp(-B/T),A=1.7e-33,B=-1000
HO2+HO2+H2O->H2O2         k(T)=Aexp(-B/T),A=3.22e-34,B=-2800
HO2+HO2+H2O+M->H2O2       k(T)=Aexp(-B/T),A=2.38e-54,B=-3200
NO3+HO2->NO2+OH           k=4e-12
NO3+H2O2->HNO3+HO2        k=2e-15
NO3+O3->NO2               k=1e-17
NO3+CH3O2->HCHO+HO2+NO2   k=1.3e-12
NO3+NO->2*NO2             k(T)=Aexp(-B/T),A=1.8e-11,B=-110
NO3+NO3->2*NO2            k(T)=Aexp(-B/T),A=8.5e-13,B=2450
NO2+NO3->NO+NO2           k(T)=Aexp(-B/T),A=7.2e-14,B=1414
NO2+NO3+M->N2O5           k(T,M)=troe(2.7e-30,0,3.4,2e-12,0,-0.2,0.33)
N2O5+M->NO3+NO2           k(T,M)=troe(1e-3,11000,3.5,9.7e14,11080,-0.1,0.33)
N2O5+H2O->2*HNO3          k=2.6e-22
N2O5+H2O+H2O->2*HNO3      k=2e-39
NO2+HO2+M->HNO4           k(T,M)=troe(1.8e-31,0,3.2,4.7e-12,0,0,0.6)
HNO4+M->HO2+NO2           k(T,M)=troe(5e-6,10000,0,2.6e15,10900,0,0.6)
HNO4+OH->NO2              k(T)=Aexp(-B/T),A=1.5e-12,B=-360

#################
# SOx chemistry #
#################

SO2+CH3O2->H2SO4+HCHO+HO2 k=4e-17
SO2+OH+M->H2SO4+HO2       k(T,M)=troe(4e-31,0,3.3,2e-12,0,0,0.45)

########################################
#      OH attack to hydrocarbons       #
########################################

CH4+OH->CH3O2             k(T)=Aexp(-B/T),A=2.3e-12,B=1765
C2H6+OH->C2H5O2           k(T)=Aexp(-B/T),A=7.9e-12,B=1030
NC4H10+OH->SECC4H         k(T)=Aexp(-B/T)(300/T)**N,A=1.36e-12,B=-190,N=-2
C2H4+OH+M->CH2O2C         k(T,M)=troe(7e-29,0,3.1,9e-12,0,0,0.7)
C3H6+OH+M->CH3CHX         k(T,M)=troe(8e-27,0,3.5,3e-11,0,0,0.5)
OXYL+OH->OXYL1            k=1.37e-11
C5H8+OH->RO2IP1           k(T)=Aexp(-B/T),A=2.55e-11,B=-410
CH3OH+OH->HCHO+HO2        k(T)=Aexp(-B/T),A=3.1e-12,B=360
C2H5OH+OH->CH3CHO+HO2     k(T)=Aexp(-B/T),A=4.1e-12,B=70

##################################################
#      OH attack to carbonyls and peroxides      #
##################################################

HCHO+OH->CO+HO2             k(T)=Aexp(-B/T),A=8.6e-12,B=-20
CH3CHO+OH->CH3COO           k(T)=Aexp(-B/T),A=5.6e-12,B=-310
CH3COE+OH->CH3COX           k(T)=Aexp(-B/T)(300/T)**N,A=2.92e-13,B=-414,N=-2
MEMALD+OH->MEMAL1           k=5.6e-11
GLYOX+OH->2*CO+HO2          k=1.1e-11
MGLYOX+OH->CH3COO+CO        k=1.5e-11
MAC+OH->0.5*MACO2           k(T)=Aexp(-B/T),A=1.86e-11,B=-175
MVK+OH->MVKO2               k(T)=Aexp(-B/T),A=4.1e-12,B=-453

CH3O2H+OH->CH3O2            k(T)=Aexp(-B/T),A=1.9e-12,B=-190
ETO2H+OH->C2H5O2            k(T)=Aexp(-B/T),A=1.9e-12,B=-190
PPA+OH->CH3COO              k(T)=Aexp(-B/T),A=1.9e-12,B=-190
BUO2H+OH->SECC4H            k(T)=Aexp(-B/T),A=1.9e-12,B=-190
CH4COX+OH->CH3COX           k(T)=Aexp(-B/T),A=1.9e-12,B=-190
CH3O2C+OH->CH2O2C           k(T)=Aexp(-B/T),A=1.9e-12,B=-190
CH4CHX+OH->CH3CHX           k(T)=Aexp(-B/T),A=1.9e-12,B=-190
OXYLH1+OH->OXYL1            k(T)=Aexp(-B/T),A=1.9e-12,B=-190
MEMALH+OH->MEMAL1           k(T)=Aexp(-B/T),A=1.9e-12,B=-190
PEROX+OH->RO2IP1            k(T)=Aexp(-B/T),A=1.9e-12,B=-190
MACO2H+OH->MACO2            k(T)=Aexp(-B/T),A=1.9e-12,B=-190
MVKO2H+OH->MVKO2            k(T)=Aexp(-B/T),A=1.9e-12,B=-190

CH3O2H+OH->HCHO+OH          k(T)=Aexp(-B/T),A=1.e-12,B=-190
ETO2H+OH->CH3CHO+OH         k(T)=Aexp(-B/T)(300/T)**N,A=4.02e-13,B=227,N=-2
ETO2H+OH->CH3CHO+OH         k(T)=Aexp(-B/T)(300/T)**N,A=1.39e-13,B=-1168,N=-2
BUO2H+OH->CH3COE+OH         k(T)=Aexp(-B/T)(300/T)**N,A=8.05e-13,B=227,N=-2
BUO2H+OH->CH3COE+OH         k(T)=Aexp(-B/T)(300/T)**N,A=3.89e-13,B=-309,N=-2
BUO2H+OH->CH3COE+OH         k(T)=Aexp(-B/T)(300/T)**N,A=6.08e-14,B=-1722,N=-2
CH4COX+OH->CH3COY+OH        k(T)=Aexp(-B/T)(300/T)**N,A=4.02e-13,B=385,N=-2
CH4COX+OH->CH3COY+OH        k(T)=Aexp(-B/T)(300/T)**N,A=4.02e-13,B=-139,N=-2
CH4COX+OH->CH3COY+OH        k(T)=Aexp(-B/T)(300/T)**N,A=6.08e-14,B=-1564,N=-2
CH3O2C+OH->2*HCHO+OH        k(T)=Aexp(-B/T)(300/T)**N,A=1.39e-13,B=-1244,N=-2
CH3O2C+OH->2*HCHO+OH        k(T)=Aexp(-B/T)(300/T)**N,A=3.89e-13,B=-674,N=-2
CH4CHX+OH->HCHO+CH3CHO+OH   k(T)=Aexp(-B/T)(300/T)**N,A=4.02e-13,B=227,N=-2
CH4CHX+OH->HCHO+CH3CHO+OH   k(T)=Aexp(-B/T)(300/T)**N,A=3.89e-13,B=-674,N=-2
CH4CHX+OH->HCHO+CH3CHO+OH   k(T)=Aexp(-B/T)(300/T)**N,A=6.08e-14,B=-1722,N=-2
OXYLH1+OH->MEMALD+MGLYOX+OH k(T)=Aexp(-B/T)(300/T)**N,A=6.08e-14,B=-1722,N=-2
OXYLH1+OH->MEMALD+MGLYOX+OH k(T)=Aexp(-B/T)(300/T)**N,A=6.08e-14,B=-1798,N=-2
OXYLH1+OH->MEMALD+MGLYOX+OH k(T)=Aexp(-B/T)(300/T)**N,A=1.70e-13,B=-1228,N=-2
MEMALH+OH->GLYOX+MGLYOX+OH  k(T)=Aexp(-B/T)(300/T)**N,A=4.02e-13,B=385,N=-2
MEMALH+OH->GLYOX+MGLYOX+OH  k(T)=Aexp(-B/T)(300/T)**N,A=1.70e-13,B=-1436,N=-2
MEMALH+OH->GLYOX+MGLYOX+OH  k(T)=Aexp(-B/T)(300/T)**N,A=6.08e-14,B=-2006,N=-2
MEMALH+OH->GLYOX+MGLYOX+OH  k(T)=Aexp(-B/T)(300/T)**N,A=1.70e-13,B=-1435,N=-2
PEROX+OH->0.43*MAC+0.57*MVK+HCHO+HO2                             k=6.1e-11
MACO2H+OH->0.5*CH3COE+0.5*MGLYOX+HO2+0.5*CO                      k=5.5e-11
MVKO2H+OH->0.28*HCHO+0.28*MGLYOX+0.28*OH+0.72*CH3CHO+0.72*CH3COO k(T)=Aexp(-B/T)(300/T)**N,A=4.02e-13,B=385,N=-2
MVKO2H+OH->0.28*HCHO+0.28*MGLYOX+0.28*OH+0.72*CH3CHO+0.72*CH3COO k(T)=Aexp(-B/T)(300/T)**N,A=3.89e-13,B=-1040,N=-2
MVKO2H+OH->0.28*HCHO+0.28*MGLYOX+0.28*OH+0.72*CH3CHO+0.72*CH3COO k(T)=Aexp(-B/T)(300/T)**N,A=6.08e-14,B=-1640,N=-2

########################################################
#      NO3 attack to hydrocarbons and aldehydes        #
########################################################

CH4+NO3->CH3O2+HNO3        k=4e-19
C2H6+NO3->C2H5O2+HNO3      k=8e-18
NC4H10+NO3->SECC4H+HNO3    k=5.5e-17
C2H4+NO3->RNC2H4           k=2e-16
C3H6+NO3->RNC3H6           k=9.45e-15
OXYL+NO3->ROXYL1           k=3.7e-16
C5H8+NO3->RNC5H8           k=7.8e-13
CH3OH+NO3->CH3O2+HNO3      k(T)=Aexp(-B/T),A=1.3e-12,B=2560
C2H5OH+NO3->C2H5O2+HNO3    k=9e-16
HCHO+NO3->CO+HNO3+HO2      k=6.8e-16
CH3CHO+NO3->CH3COO+HNO3    k=2.8e-15

########################################################
#      O3 attack to unsaturated hydrocarbons           #
########################################################

C2H4+O3->HCHO+0.12*HO2+0.13*H2+0.44*CO                                            k(T)=Aexp(-B/T),A=9.1e-15,B=2580
C3H6+O3->0.53*HCHO+0.5*CH3CHO+0.31*CH3O2+0.28*HO2+0.15*OH+0.065*H2+0.4*CO+0.7*CH4 k(T)=Aexp(-B/T),A=5.5e-15,B=1880
C5H8+O3->0.67*MAC+0.26*MVK+0.55*OH+0.07*C3H6+0.8*HCHO+0.06*HO2+0.05*CO            k(T)=Aexp(-B/T),A=1.2e-14,B=2013
MAC+O3->0.8*MGLYOX+0.7*HCHO+0.215*OH+0.275*HO2+0.2*CO                             k(T)=Aexp(-B/T),A=5.3e-15,B=2520
MVK+O3->0.82*MGLYOX+0.8*HCHO+0.04*CH3CHO+0.08*OH+0.06*HO2+0.05*CO                 k(T)=Aexp(-B/T),A=4.3e-15,B=2016
PEROX+O3->0.7*HCHO                                                                k=8e-18

#########################################
#      organic radical conversion       #
#########################################

CH3O2+NO->HCHO+NO2+HO2          k(T)=Aexp(-B/T),A=4.2e-12,B=-180
C2H5O2+NO->CH3CHO+NO2+HO2       k=8.7e-12
CH3COO+NO->CH3O2+NO2+CO2        k=2e-11
CH2O2C+NO->2*HCHO+NO2+HO2       k=9e-12
SECC4H+NO->CH3COE+HO2+NO2       k=4e-12
SECC4H+NO+M->CH3COE+HO2+NO2     k(T,M)=troe(-3.76e-32,0,0,-3.3e-12,0,8.1,0.41)
CH3COX+NO->CH3COY+NO2+HO2       k=4e-12
CH3CHX+NO->HCHO+CH3CHO+NO2+HO2  k=4e-12
OXYL1+NO->MEMALD+MGLYOX+NO2+HO2 k=4e-12
MEMAL1+NO->GLYOX+MGLYOX+NO2+HO2 k=4e-12
MACO2+NO->CH3COZ+NO2+CO2        k=1.4e-11
RO2IP1+NO->0.32*MAC+0.42*MVK+0.74*HCHO+0.14*ISNI+0.12*RO2IP1+0.78*HO2+0.86*NO2          k(T)=Aexp(-B/T),A=1.4e-11,B=180
MVKO2+NO->0.266*MGLYOX+0.266*HCHO+0.684*CH3CHO+0.684*CH3COO+0.05*ISNI+0.95*HO2+0.95*NO2 k(T)=Aexp(-B/T),A=1.4e-11,B=180
CH3COZ+NO->0.25*CH3COE+0.75*MGLYOX+HO2+NO2                                              k(T)=Aexp(-B/T),A=1.4e-11,B=180

############################################
#      organic radical recombination       #
############################################

CH3O2+HO2->CH3O2H            k(T)=Aexp(-B/T),A=4.1e-13,B=-790
C2H5O2+HO2->ETO2H            k(T)=Aexp(-B/T),A=6.5e-13,B=-650
CH3COO+HO2->0.67*PPA+0.33*O3 k(T)=Aexp(-B/T),A=4.3e-13,B=-1040
SECC4H+HO2->BUO2H            k(T)=Aexp(-B/T),A=2.6e-13,B=-1245
CH3COX+HO2->CH4COX           k(T)=Aexp(-B/T),A=3.5e-13,B=-1000
CH2O2C+HO2->CH3O2C           k(T)=Aexp(-B/T),A=3.5e-13,B=-1000
CH3CHX+HO2->CH4CHX           k(T)=Aexp(-B/T),A=3.5e-13,B=-1000
OXYL1+HO2->OXYLH1            k(T)=Aexp(-B/T),A=3.5e-13,B=-1000
MEMAL1+HO2->MEMALH           k(T)=Aexp(-B/T),A=3.5e-13,B=-1000
RO2IP1+HO2->PEROX            k(T)=Aexp(-B/T),A=3.5e-13,B=-1000
MACO2+HO2->MACO2H            k(T)=Aexp(-B/T),A=3.5e-13,B=-1000
MVKO2+HO2->MVKO2H            k(T)=Aexp(-B/T),A=3.5e-13,B=-1000

CH3O2+CH3O2->2*HCHO+2*HO2              k(T)=SPECIAL_1(9.1e-14,-416,25,1165)
CH3O2+CH3O2->HCHO+CH3OH                k(T)=SPECIAL_2(9.1e-14,-416,25,1165)
C2H5O2+CH3O2->CH3CHO+HCHO+2*HO2        k(T)=SPECIAL_3(9.8e-14,100,9.1e-14,-416,10.2,533,25,1165)
C2H5O2+CH3O2->CH3CHO+CH3OH             k(T)=SPECIAL_4(9.8e-14,100,9.1e-14,-416,10.2,533,25,1165)
C2H5O2+CH3O2->C2H5OH+HCHO              k(T)=SPECIAL_4(9.8e-14,100,9.1e-14,-416,25,1165,10.2,533)
CH3COO+CH3O2->CH3O2+CO2+HCHO+HO2       k(T)=SPECIAL_1(5.1e-12,-272,4.4e5,3910)
CH3COO+CH3O2->MECOOH+HCHO              k(T)=SPECIAL_2(5.1e-12,-272,4.4e5,3910)
SECC4H+CH3O2->CH3COE+HCHO+2*HO2        k(T)=SPECIAL_3(1.7e-12,2188,9.1e-14,-416,50,1165,25,1165)
SECC4H+CH3O2->CH3COE+CH3OH             k(T)=SPECIAL_4(1.7e-12,2188,9.1e-14,-416,50,1165,25,1165)
SECC4H+CH3O2->CH3COE+HCHO              k(T)=SPECIAL_4(1.7e-12,2188,9.1e-14,-416,25,1165,50,1165)
CH3COX+CH3O2->CH3COY+HCHO+2*HO2        k(T)=SPECIAL_3(1.2e-13,0,9.1e-14,-416,50,1165,25,1165)
CH3COX+CH3O2->CH3COY+CH3OH             k(T)=SPECIAL_4(1.2e-13,0,9.1e-14,-416,50,1165,25,1165)
CH3COX+CH3O2->CH3COY+HCHO              k(T)=SPECIAL_4(1.2e-13,0,9.1e-14,-416,25,1165,50,1165)
CH2O2C+CH3O2->3*HCHO+2*HO2             k(T)=SPECIAL_3(1.9e-13,-745,9.1e-14,-416,50,1165,25,1165)
CH2O2C+CH3O2->2*HCHO+CH3OH             k(T)=SPECIAL_4(1.9e-13,-745,9.1e-14,-416,50,1165,25,1165)
CH2O2C+CH3O2->3*HCHO                   k(T)=SPECIAL_4(1.9e-13,-745,9.1e-14,-416,25,1165,50,1165)
CH3CHX+CH3O2->CH3CHO+2*HCHO+2*HO2      k(T)=SPECIAL_3(3.5e-14,0,9.1e-14,-416,50,1165,25,1165)
CH3CHX+CH3O2->CH3CHO+HCHO+CH3OH        k(T)=SPECIAL_4(3.5e-14,0,9.1e-14,-416,50,1165,25,1165)
CH3CHX+CH3O2->CH3CHO+2*HCHO            k(T)=SPECIAL_4(3.5e-14,0,9.1e-14,-416,25,1165,50,1165)
OXYL1+CH3O2->MEMALD+MGLYOX+HCHO+2*HO2  k(T)=SPECIAL_3(1.2e-13,-745,9.1e-14,-416,50,1165,25,1165)
OXYL1+CH3O2->MEMALD+MGLYOX+CH3OH       k(T)=SPECIAL_4(1.2e-13,-745,9.1e-14,-416,50,1165,25,1165)
OXYL1+CH3O2->MEMALD+MGLYOX+HCHO        k(T)=SPECIAL_4(1.2e-13,-745,9.1e-14,-416,25,1165,50,1165)
MEMAL1+CH3O2->GLYOX+MGLYOX+HCHO+2*HO2  k(T)=SPECIAL_3(1.2e-13,0,9.1e-14,-416,50,1165,25,1165)
MEMAL1+CH3O2->GLYOX+MGLYOX+CH3OH       k(T)=SPECIAL_4(1.2e-13,0,9.1e-14,-416,50,1165,25,1165)
MEMAL1+CH3O2->GLYOX+MGLYOX+HCHO        k(T)=SPECIAL_4(1.2e-13,0,9.1e-14,-416,25,1165,50,1165)
RO2IP1+CH3O2->0.3*MVK+0.3*MAC+1.6*HCHO+2*HO2 k(T)=SPECIAL_3(3.0e-13,0,9.1e-14,-416,50,1165,25,1165)
RO2IP1+CH3O2->0.3*MVK+0.3*MAC+0.6*HCHO+CH3OH k(T)=SPECIAL_4(3.0e-13,0,9.1e-14,-416,50,1165,25,1165)
RO2IP1+CH3O2->0.3*MVK+0.3*MAC+1.6*HCHO k(T)=SPECIAL_4(3.0e-13,0,9.1e-14,-416,25,1165,50,1165)
MACO2+CH3O2->CH3COZ+CO2+HCHO+HO2       k(T)=SPECIAL_1(5.1e-12,-272,4.4e5,3910)
MACO2+CH3O2->MAA+HCHO                  k(T)=SPECIAL_2(5.1e-12,-272,4.4e5,3910)
MVKO2+CH3O2->MGLYOX+2*HCHO+2*HO2       k(T)=SPECIAL_3(3.8e-14,0,9.1e-14,-416,50,1165,25,1165)
MVKO2+CH3O2->MGLYOX+HCHO+CH3OH         k(T)=SPECIAL_4(3.8e-14,0,9.1e-14,-416,50,1165,25,1165)
MVKO2+CH3O2->MGLYOX+2*HCHO             k(T)=SPECIAL_4(3.8e-14,0,9.1e-14,-416,25,1165,50,1165)

CH3COO+CH3COO->2.*CH3O2+2.*CO2            k(T)=Aexp(-B/T),A=2.8e-12,B=-530
C2H5O2+CH3COO->CH3CHO+CH3O2+CO2+HO2       k(T)=SPECIAL_3(9.8e-14,100,2.8e-12,-530,10.2,533,1e9,0)
C2H5O2+CH3COO->CH3CHO+MECOOH              k(T)=SPECIAL_4(9.8e-14,100,2.8e-12,-530,1e9,0,10.2,533)
SECC4H+CH3COO->CH3COE+CH3O2+CO2+HO2       k(T)=SPECIAL_3(1.7e-12,2188,2.8e-12,-530,1e9,0,50,1165)
SECC4H+CH3COO->CH3COE+MECOOH              k(T)=SPECIAL_4(1.7e-12,2188,2.8e-12,-530,1e9,0,50,1165)
CH3COX+CH3COO->CH3COY+CH3O2+CO2+HO2       k(T)=SPECIAL_3(1.2e-13,0,2.8e-12,-530,1e9,0,50,1165)
CH3COX+CH3COO->CH3COY+MECOOH              k(T)=SPECIAL_4(1.2e-13,0,2.8e-12,-530,1e9,0,50,1165)
CH2O2C+CH3COO->2*HCHO+CH3O2+CO2+HO2       k(T)=SPECIAL_3(1.9e-13,-745,2.8e-12,-530,1e9,0,50,1165)
CH2O2C+CH3COO->2*HCHO+MECOOH              k(T)=SPECIAL_4(1.9e-13,-745,2.8e-12,-530,1e9,0,50,1165)
CH3CHX+CH3COO->CH3CHO+HCHO+CH3O2+CO2+HO2  k(T)=SPECIAL_3(3.5e-14,0,2.8e-12,-530,1e9,0,50,1165)
CH3CHX+CH3COO->CH3CHO+HCHO+MECOOH         k(T)=SPECIAL_4(3.5e-14,0,2.8e-12,-530,1e9,0,50,1165)
OXYL1+CH3COO->MEMALD+MGLYOX+CH3O2+CO2+HO2 k(T)=SPECIAL_3(1.2e-13,-745,2.8e-12,-530,1e9,0,50,1165)
OXYL1+CH3COO->MEMALD+MGLYOX+MECOOH        k(T)=SPECIAL_4(1.2e-13,-745,2.8e-12,-530,1e9,0,50,1165)
MEMAL1+CH3COO->GLYOX+MGLYOX+CH3O2+CO2+HO2 k(T)=SPECIAL_3(1.2e-13,0,2.8e-12,-530,1e9,0,50,1165)
MEMAL1+CH3COO->GLYOX+MGLYOX+MECOOH        k(T)=SPECIAL_4(1.2e-13,0,2.8e-12,-530,1e9,0,50,1165)
RO2IP1+CH3COO->0.3*MVK+0.3*MAC+0.6*HCHO+HO2+CH3O2+CO2  k(T)=SPECIAL_3(3e-13,0,2.8e-12,-530,1e9,0,50,1165)
RO2IP1+CH3COO->0.3*MVK+0.3*MAC+0.6*HCHO+MECOOH         k(T)=SPECIAL_4(3e-13,0,2.8e-12,-530,1e9,0,50,1165)
MACO2+CH3COO->CH3O2+CH3COZ+2*CO2          k(T)=Aexp(-B/T),A=5.6e-12,B=-530
MVKO2+CH3COO->MGLYOX+HCHO+HO2+CH3O2+CO2   k(T)=SPECIAL_3(3.8e-14,0,2.8e-12,-530,1e9,0,50,1165)
MVKO2+CH3COO->MGLYOX+HCHO+MECOOH          k(T)=SPECIAL_4(3.8e-14,0,2.8e-12,-530,1e9,0,50,1165)

C2H5O2+C2H5O2->2*CH3CHO+2*HO2                   k(T)=SPECIAL_1(9.8e-14,100,10.2,533)
C2H5O2+C2H5O2->CH3CHO+C2H5OH                    k(T)=SPECIAL_2(9.8e-14,100,10.2,533)
SECC4H+SECC4H->2*CH3COE+2*HO2                   k(T)=SPECIAL_1(1.7e-12,2188,50,1165)
SECC4H+SECC4H->2*CH3COE                         k(T)=SPECIAL_2(1.7e-12,2188,50,1165)
CH3COX+CH3COX->2*CH3COY+2*HO2                   k(T)=SPECIAL_1(1.2e-13,0,50,1165)
CH3COX+CH3COX->2*CH3COY                         k(T)=SPECIAL_2(1.2e-13,0,50,1165)
CH2O2C+CH2O2C->4*HCHO+2*HO2                     k(T)=SPECIAL_1(1.9e-13,-745,50,1165)
CH2O2C+CH2O2C->4*HCHO                           k(T)=SPECIAL_2(1.9e-13,-745,50,1165)
CH3CHX+CH3CHX->2*HCHO+2*CH3CHO+2*HO2            k(T)=SPECIAL_1(3.5e-14,0,50,1165)
CH3CHX+CH3CHX->2*HCHO+2*CH3CHO                  k(T)=SPECIAL_2(3.5e-14,0,50,1165)
OXYL1+OXYL1->2*MEMALD+2*MGLYOX+2*HO2            k(T)=SPECIAL_1(1.2e-13,-745.,50,1165)
OXYL1+OXYL1->2*MEMALD+2*MGLYOX                  k(T)=SPECIAL_2(1.2e-13,-745.,50,1165)
MEMAL1+MEMAL1->2*GLYOX+2*MGLYOX+2*HO2           k(T)=SPECIAL_1(1.2e-13,0,50,1165)
MEMAL1+MEMAL1->2*GLYOX+2*MGLYOX                 k(T)=SPECIAL_2(1.2e-13,0,50,1165)
RO2IP1+RO2IP1->0.6*MVK+0.6*MAC+1.2*HCHO+1.2*HO2 k(T)=SPECIAL_1(3.0e-13,0,50,1165)
RO2IP1+RO2IP1->0.6*MVK+0.6*MAC+1.2*HCHO         k(T)=SPECIAL_2(3.0e-13,0,50,1165)
MACO2+MACO2->2*CH3COZ+2*CO2                     k(T)=Aexp(-B/T),A=2.8e-12,B=-530
MVKO2+MVKO2->2*MGLYOX+2*HCHO+2*HO2              k(T)=SPECIAL_1(3.8e-14,0,50,1165)
MVKO2+MVKO2->2*MGLYOX+2*HCHO                    k(T)=SPECIAL_2(3.8e-14,0,50,1165)

##############################################
#      organic nitrates and pernitrates      #
##############################################

CH3O2+NO->CH3NO3          k(T)=Aexp(-B/T),A=2.1e-14,B=-180
CH3NO3+OH->HCHO+HNO3+HO2  k(T)=Aexp(-B/T),A=1.e-14,B=-1060
C2H5O2+NO->ETNO3          k=1.25e-13
ETNO3+OH->CH3CHO+HNO3+HO2 k(T)=Aexp(-B/T),A=4.4e-14,B=-720
SECC4H+NO+M->BUNO3        k(T,M)=troe(3.76e-32,0,0,3.3e-12,0,8.1,0.41)
BUNO3+OH->RNC3H6          k=1.63e-12
RNC2H4+NO->ORGNIT         k=8.9e-12
RNC3H6+NO->ORGNIT         k=4e-12
ROXYL1+NO->ORGNIT         k=4e-12
ISNI+OH->ISNIR            k=3.4e-11
ISNI+O3->0.5*ISNI+0.5*HCHO+0.5*MGLYOX+0.5*NO2+0.08*OH   k=5e-18
ISNIR+NO->0.95*CH3CHO+0.475*CH3COE+0.475*MGLYOX+0.05*ISNI+0.05*HO2+1.9*NO2   k(T)=Aexp(-B/T),A=1.4e-11,B=180
RNC5H8+NO->0.85*ISNI+0.1*MAC+0.05*MVK+0.15*HCHO+1.1*NO2+0.8*HO2   k(T)=Aexp(-B/T),A=1.4e-11,B=180

CH3COO+NO2+M->PAN     k(T,M)=troe(2.7e-28,0,7.1,1.2e-11,0,0.9,0.3)
PAN+M->CH3COO+NO2     k(T,M)=troe(4.9e-3,12100,0,5.4e16,13830,0,0.3)
PAN+OH->HCOOH+NO3+CO2 k(T)=Aexp(-B/T),A=9.5e-13,B=650
MACO2+NO2->MACPAN     k=8.4e-12
MACPAN->MACO2+NO2     k(T)=Aexp(-B/T),A=1.6e17,B=13507
CH3O2+NO2+M->CH3NO4   k(T,M)=troe(2.5e-30,0,5.5,7.5e-12,0,0,0.36)
CH3NO4+M->CH3O2+NO2   k(T,M)=troe(9e-5,9690,0,1.1e16,10560,0,0.36)

##############################
#    photolytic reactions    #
##############################

O3->2*OH                                                       J(T,Z,H2O)=photorate
NO2->NO+O3                                                     J(Z)=photorate
NO3->NO2+O3                                                    J(Z)=photorate
H2O2->2*OH                                                     J(Z)=photorate
HNO3->NO2+OH                                                   J(Z)=photorate
HCHO->CO+2*HO2                                                 J(Z)=photorate
HCHO->CO+H2                                                    J(Z)=photorate
CH3CHO->CH3O2+HO2+CO                                           J(Z)=photorate
CH3COE->CH3COO+C2H5O2                                          J(Z)=photorate
CH3COY->2*CH3COO                                               J(Z)=photorate
MGLYOX->CH3COO+HO2+CO                                          J(Z)=photorate
GLYOX->0.6*HO2+2*CO+1.4*H2                                     J(Z)=photorate
NO3->NO                                                        J(Z)=photorate
N2O5->NO2+NO3                                                  J(Z)=photorate
CH3O2H->HCHO+OH+HO2                                            J(Z)=photorate
ETO2H->CH3CHO+OH+HO2                                           J(Z)=photorate
BUO2H->CH3COE+OH+HO2                                           J(Z)=photorate
CH4COX->CH3COY+OH+HO2                                          J(Z)=photorate
CH3O2C->2*HCHO+OH+HO2                                          J(Z)=photorate
CH4CHX->HCHO+CH3CHO+OH+HO2                                     J(Z)=photorate
OXYLH1->MEMALD+MGLYOX+OH+HO2                                   J(Z)=photorate
MEMALH->GLYOX+MGLYOX+OH+HO2                                    J(Z)=photorate
PEROX->0.14*RO2IP1+0.488*CH3COE+0.372*MAC+0.86*HCHO+0.9*HO2+OH J(Z)=photorate
MACO2H->CH3COZ+OH+CO2                                          J(Z)=photorate
MVKO2H->0.28*HCHO+0.28*MGLYOX+0.72*CH3CHO+0.72*CH3COO          J(Z)=photorate
PPA->CH3O2+CO2+OH                                              J(Z)=photorate
HNO4->HO2+NO2                                                  J(Z)=photorate
CH3NO4->CH3O2+NO2                                              J(Z)=photorate
CH3NO3->HCHO+HO2+NO2                                           J(Z)=photorate
ETNO3->CH3CHO+HO2+NO2                                          J(Z)=photorate
BUNO3->CH3COE+HO2+NO2                                          J(Z)=photorate
PAN->CH3COO+NO2                                                J(Z)=photorate

############
# Families #
############

NOX=NO+NO2
NOX->X     k=1e-20

OX=NO2+O3
OX->X      k=1e-20

HOX=OH+HO2+CH3O2+C2H5O2+CH3COO+CH2O2C+SECC4H+CH3COX+CH3CHX+OXYL1+MEMAL1+MACO2+RO2IP1+MVKO2+CH3COZ+RNC2H4+RNC3H6+ROXYL1+ISNIR
HOX->X     k=1e-20
