02 jan 2015
+---------+
NOTE ON THE SULFUR_AQ CHEM:

only aqueous phase reactions implemented, no gas-phase chemistry, it is assumed CB05 solves the chemistry

all the reaction implemented are from Takemura et al. (JGR2000) [SPRINTARS]
[JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 105, NO. D14, PAGES 17,853-17,873,J ULY 27, 2000]

with the exception of reactions:
DMS -> SO2, SO4 (which was R5 and R6)
replaced by:
DMS + OH -> 0.75 SO2 + 0.25 MSA (Tie et al., JGR2001) (R5) [MOZART, ECHAM-HAM (Feichter et al., AE2006)]
[JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 106, NO. D19, PAGES 22,931-22,964, OCTOBER 16, 2001]

03 nov 2016
+---------+
NOTE ON THE SULFUR_AQ CHEM:

implemented dms gas-phase (Emmons et al. 2010), sulfur aqueous chemistry (Takemura et al., 2000), and heterogeneous hydrolisis of N2O5 (caaba/mecca)

     '            OH + DMS --> SO2                                                                        ', &
     '            OH + DMS --> 0.5 HO2 + 0.5 SO2                                                          ', &
     '  O3_a01 + SO3mm_a01 --> SO4mm_a01                                                                  ', &
     '  O3_a01 + HSO3m_a01 --> SO4mm_a01 + Hp_a01                                                         ', &
     'H2O2_a01 + HSO3m_a01 --> SO4mm_a01 + Hp_a01                                                         ', &
     '                  O3 --> O3_a01                                                                     ', &
     '              O3_a01 --> O3                                                                         ', &
     '                H2O2 --> H2O2_a01                                                                   ', &
     '            H2O2_a01 --> H2O2                                                                       ', &
     '                N2O5 --> 2 HNO3_a01                                                                 ', &
     '                 SO2 --> SO2_a01                                                                    ', &
     '             SO2_a01 --> SO2                                                                        ', &
     '             SO2_a01 --> HSO3m_a01 + Hp_a01                                                         ', &
     '  HSO3m_a01 + Hp_a01 --> SO2_a01                                                                    ', &
     '           HSO3m_a01 --> Hp_a01 + SO3mm_a01                                                         ', &
     '  Hp_a01 + SO3mm_a01 --> HSO3m_a01                                                                  ' /)

  CHARACTER(LEN=32), PARAMETER, DIMENSION(16) :: EQN_TAGS = (/ &
     'G9400aG9400aOJ                  ','G9400bG9400bOJ                  ','A9101_a01                       ', &
     'A9206_a01                       ','A9209_a01                       ','H1001f_a01                      ', &
     'H1001b_a01                      ','H2102f_a01                      ','H2102b_a01                      ', &
     'H3201_a01                       ','H9100f_a01                      ','H9100b_a01                      ', &
     'EQ90f_a01                       ','EQ90b_a01                       ','EQ91f_a01                       ', &
     'EQ91b_a01                       ' /)

27 mar 2017
+---------+
NOTE ON THE SULFUR_AQ CHEM:

implemented dms gas-phase (Emmons et al. 2010), sulfur aqueous chemistry (Takemura et al., 2000), and heterogeneous hydrolisis of N2O5 (caaba/mecca)
correct error in previous implementation, removed H3201_a01 and added HET201

     '            OH + DMS --> SO2                                                                        ', &
     '            OH + DMS --> 0.5 HO2 + 0.5 SO2                                                          ', &
     '                N2O5 --> 2 NO3m_cs + 2 Hp_cs                                                        ', &
     '  O3_a01 + SO3mm_a01 --> SO4mm_a01                                                                  ', &
     '  HSO3m_a01 + O3_a01 --> SO4mm_a01 + Hp_a01                                                         ', &
     'H2O2_a01 + HSO3m_a01 --> SO4mm_a01 + Hp_a01                                                         ', &
     '                  O3 --> O3_a01                                                                     ', &
     '              O3_a01 --> O3                                                                         ', &
     '                H2O2 --> H2O2_a01                                                                   ', &
     '            H2O2_a01 --> H2O2                                                                       ', &
     '                 SO2 --> SO2_a01                                                                    ', &
     '             SO2_a01 --> SO2                                                                        ', &
     '             SO2_a01 --> Hp_a01 + HSO3m_a01                                                         ', &
     '  Hp_a01 + HSO3m_a01 --> SO2_a01                                                                    ', &
     '           HSO3m_a01 --> Hp_a01 + SO3mm_a01                                                         ', &
     '  Hp_a01 + SO3mm_a01 --> HSO3m_a01                                                                  ' /)

  CHARACTER(LEN=32), PARAMETER, DIMENSION(16) :: EQN_TAGS = (/ &
     'G9400aG9400aOJ                  ','G9400bG9400bOJ                  ','HET201                          ', &
     'A9101_a01                       ','A9206_a01                       ','A9209_a01                       ', &
     'H1001f_a01                      ','H1001b_a01                      ','H2102f_a01                      ', &
     'H2102b_a01                      ','H9100f_a01                      ','H9100b_a01                      ', &
     'EQ90f_a01                       ','EQ90b_a01                       ','EQ91f_a01                       ', &
     'EQ91b_a01                       ' /)

23 jan 2018
+---------+
NOTE ON THE SULFUR_AQ CHEM:
change OH from variable to FIX
