Changeset 3141

Timestamp:

Nov 27, 2023, 11:41:11 AM (13 months ago)

Author:

flefevre

Message:

Correction of a duplicated ionic reaction (i026 and i030) + a bit of cosmetics

File:

• trunk/LMDZ.MARS/libf/aeronomars/photochemistry.F90 (modified) (58 diffs)

trunk/LMDZ.MARS/libf/aeronomars/photochemistry.F90

@@ -1869 +1869 @@
          v_4(:,nb_reaction_4) = i025(:)
 
-!---     i026: HCO2+ + e -> H + CO2
-
-!        UMIST
-
-         i026(:) = 1.75e-8*((300./t_elect(:))**0.5)
+!---     i026: spare slot
+
+         i026(:) = 0.
 
          nb_reaction_4 = nb_reaction_4 + 1
@@ -1903 +1901 @@
 !---     i030: HCO2+ + e -> H + CO2
 
+!        UMIST
+
          i030(:) = 6.0e-8*((300./t_elect(:))**0.64)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i030(:)
@@ -1912 +1913 @@
 
          i031(:) = 1.e-9
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i031(:)
@@ -1920 +1922 @@
 
          i032(:) = 7.8e-10
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i032(:)
@@ -1928 +1931 @@
 
          i033(:) = 3.5e-9
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i033(:)
 
-
 !---     i034: CO2+ + H -> HCO+ + O
 
@@ -1945 +1948 @@
 
          i035(:) = 7.5e-10
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i035(:)
@@ -1953 +1957 @@
 
          i036(:) = 2.4e-7 *((300./t_elect(:))**0.69)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i036(:)
@@ -1961 +1966 @@
 
          i037(:) = 2.04e-9 *((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i037(:)
@@ -1969 +1975 @@
 
          i038(:) = 1.72e-9*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i038(:)
@@ -1977 +1984 @@
 
          i039(:) = 3.2e-9*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i039(:)
@@ -1985 +1993 @@
 
          i040(:) = 2.3e-9*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i040(:)
@@ -1993 +2002 @@
 
          i041(:) = 2.8e-9*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i041(:)
 
-
 !---     i042: H+ + H2O -> H2O+ + H
 
@@ -2002 +2011 @@
 
          i042(:) = 6.9e-9*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i042(:)
@@ -2010 +2020 @@
 
          i043(:) = 4.6e-10
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i043(:)
@@ -2018 +2029 @@
 
          i044(:) = 5.0e-10
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i044(:)
@@ -2026 +2038 @@
 
          i045(:) = 4.0e-11
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i045(:)
@@ -2034 +2047 @@
 
          i046(:) = 2.7e-10
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i046(:)
@@ -2042 +2056 @@
          
          i047(:) = 3.05e-7*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i047(:)
@@ -2050 +2065 @@
 
          i048(:) = 8.6e-8*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i048(:)
@@ -2058 +2074 @@
 
          i049(:) = 3.9e-8*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i049(:)
@@ -2066 +2083 @@
 
          i050(:) = 2.1e-9*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i050(:)
 
-
 !---     i051: H2O+ + H2 -> H3O+ + H
 
@@ -2075 +2092 @@
 
          i051(:) = 6.4e-10
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i051(:)
@@ -2083 +2101 @@
 
          i052(:) = 2.5e-9*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i052(:)
@@ -2091 +2110 @@
 
          i053(:) = 3.05e-7*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i053(:)
@@ -2099 +2119 @@
          
          i054(:) = 7.09e-8*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i054(:)
@@ -2107 +2128 @@
 
          i055(:) = 5.37e-8*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i055(:)
@@ -2115 +2137 @@
 
          i056(:) = 5.6e-9*((300./t_elect(:))**0.5)
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i056(:)
@@ -2126 +2149 @@
 
          i057(:) = 1.7e-9
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i057(:)
@@ -2134 +2158 @@
 
          i058(:) = 7.1e-10
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i058(:)
@@ -2142 +2167 @@
 
          i059(:) = 1.44e-9
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i059(:)
@@ -2150 +2176 @@
 
          i060(:) = 1.05e-9
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i060(:)
@@ -2158 +2185 @@
 
          i061(:) = 3.59e-10
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i061(:)
@@ -2166 +2194 @@
 
          i062(:) = 1.01e-9
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i062(:)
@@ -2174 +2203 @@
 
          i063(:) = 5.9e-10
+
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i063(:)
@@ -3189 +3219 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_oh, 1.0, i_hd, 1.0, i_hdo, 1.0, i_h)
 
-
 !===========================================================
 !      f016 : OH + HD -> H2O + D 
@@ -3198 +3227 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_oh, 1.0, i_hd, 1.0, i_h2o, 1.0, i_d)
 
-
 !===========================================================
 !      f017 : D + O2 + CO2 -> DO2 + CO2
@@ -3207 +3235 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_d, 1.0, i_o2, 1.0, i_do2, 0.0, i_dummy)
 
-
 !===========================================================
 !      f018 : OD + O3 -> DO2 + O2 
@@ -3216 +3243 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_od, 1.0, i_o3, 1.0, i_do2, 1.0, i_o2)
 
-
 !===========================================================
 !      f019 : D + HO2 -> DO2 + H 
@@ -3225 +3251 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_d, 1.0, i_ho2, 1.0, i_do2, 1.0, i_h)
 
-
 !===========================================================
 !      f020 : O + DO2 -> OD + O2
@@ -3234 +3259 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_o, 1.0, i_do2, 1.0, i_od, 1.0, i_o2)
 
-
 !===========================================================
 !      f021 : H + DO2 -> OH + OD
@@ -3243 +3267 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_h, 1.0, i_do2, 1.0, i_od, 1.0, i_oh)
 
-
 !===========================================================
 !      f022 : H + DO2 -> HD + O2
@@ -3252 +3275 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_h, 1.0, i_do2, 1.0, i_hd, 1.0, i_o2)
 
-
 !===========================================================
 !      f023 : H + DO2 -> HDO + O
@@ -3261 +3283 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_h, 1.0, i_do2, 1.0, i_hdo, 1.0, i_o)
 
-
 !===========================================================
 !      f024 : H + DO2 -> HO2 + D
@@ -3270 +3291 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_h, 1.0, i_do2, 1.0, i_ho2, 1.0, i_d)
 
-
 !===========================================================
 !      f025 : OH + DO2 -> HDO + O2
@@ -3279 +3299 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_oh, 1.0, i_do2, 1.0, i_hdo, 1.0, i_o2)
 
-
 !===========================================================
 !      f026 : DO2 + O3 -> OD + O2 + O2 
@@ -3288 +3307 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_do2, 1.0, i_o3, 1.0, i_od, 2.0, i_o2)
 
-
 !===========================================================
 !      f027 : OD + OH + CO2 -> HDO2 + CO2 
@@ -3297 +3315 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_od, 1.0, i_oh, 1.0, i_hdo2, 0.0, i_dummy)
 
-
 !===========================================================
 !      f028 : DO2 + HO2 -> HDO2 + O2 
@@ -3314 +3331 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_o, 1.0, i_hdo2, 1.0, i_od, 1.0, i_ho2)
 
-
 !===========================================================
 !      f030 : O + HDO2 -> OH + DO2 
@@ -3323 +3339 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_o, 1.0, i_hdo2, 1.0, i_oh, 1.0, i_do2)
 
-
 !===========================================================
 !      f031 : OH + HDO2 -> HDO + HO2 
@@ -3332 +3347 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_oh, 1.0, i_hdo2, 1.0, i_hdo, 1.0, i_ho2)
 
-
 !===========================================================
 !      f032 : OH + HDO2 -> H2O + DO2 
@@ -3340 +3354 @@
 
    indice_4(nb_reaction_4) = z4spec(1.0, i_oh, 1.0, i_hdo2, 1.0, i_h2o, 1.0, i_do2)
-
 
 endif !deutchem
@@ -3548 +3561 @@
 
 !===========================================================
-!      i026 : HCO2+ + e -> H + CO2 
+!      i026 : spare slot (reaction rate set to zero)
 !===========================================================
 
@@ -3580 +3593 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_hco2plus, 1.0, i_elec, 1.0, i_oh, 1.0, i_co)
 
-
 !===========================================================
 !      i030 : HCO2+ + e -> H + CO2 
@@ -3589 +3601 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_hco2plus, 1.0, i_elec, 1.0, i_h, 1.0, i_co2)
 
-
 !===========================================================
 !      i031 : HCO2+ + O -> HCO+ + O2 
@@ -3598 +3609 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_hco2plus, 1.0, i_o, 1.0, i_hcoplus, 1.0, i_o2)
 
-
 !===========================================================
 !      i032 : HCO2+ + CO -> HCO+ + CO2 
@@ -3606 +3616 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_hco2plus, 1.0, i_co, 1.0, i_hcoplus, 1.0, i_co2)
 
-
 !===========================================================
 !      i033 : H+ + CO2 -> HCO+ + O 
@@ -3614 +3623 @@
    indice_4(nb_reaction_4) = z4spec(1.0, i_hplus, 1.0, i_co2, 1.0, i_hcoplus, 1.0, i_o)
 
-
 !===========================================================
 !      i034 : CO2+ + H -> HCO+ + O 
@@ -3621 +3629 @@
    nb_reaction_4 = nb_reaction_4 + 1
    indice_4(nb_reaction_4) = z4spec(1.0, i_co2plus, 1.0, i_h, 1.0, i_hcoplus, 1.0, i_o)
-
 
 !===========================================================