Changeset 2273 for trunk/LMDZ.MARS/libf/aeronomars

Timestamp:

Apr 2, 2020, 4:23:47 PM (5 years ago)

Author:

emillour

Message:

Mars GCM:
A first step towards water ion chemistry. Add a new reaction and modify reaction
rates for HCO2+
FGG

Location:

trunk/LMDZ.MARS/libf/aeronomars

Files:

• calchim_mod.F90 (modified) (1 diff)
• photochemistry.F90 (modified) (4 diffs)

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

@@ -630 +630 @@
             ! set number of reactions, depending on ion chemistry or not
             if (ionchem) then
-               nb_reaction_4_max = 60   ! set number of bimolecular reactions
+               nb_reaction_4_max = 61   ! set number of bimolecular reactions
                nb_reaction_3_max = 6    ! set number of quadratic reactions
                nquench           = 9    ! set number of quenching + heterogeneous reactions

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

@@ -592 +592 @@
                            i013, i014, i015, i016, i017, i018, i019,   &
                            i020, i021, i022, i023, i024, i025, i026,   &
-                           i027, i028, i029,                           &
+                           i027, i028, i029, i030,                     &
                            h001, h002, h003, h004, h005
 
@@ -1338 +1338 @@
          !i028: 0.5 (HCO2+ + e-) -> O + CO
 
-         i027(:) = 2.38e-7*((300./t_elect(:))**0.5)
+         i027(:) = 8.1e-7*((300./t_elect(:))**0.64)
 
          nb_reaction_4 = nb_reaction_4 + 1
@@ -1346 +1346 @@
          v_4(:,nb_reaction_4) = i027(:)
 
-!---     i029: HCO2+ + e -> H + CO2
+!---     i029: HCO2+ + e -> OH + CO
 
 !        UMIST
 
-         i029(:) = 9.45e-8*((300./t_elect(:))**0.5)
+         i029(:) = 3.2e-7*((300./t_elect(:))**0.64)
 
          nb_reaction_4 = nb_reaction_4 + 1
          v_4(:,nb_reaction_4) = i029(:)
+
+!---     i030: HCO2+ + e -> H + CO2
+
+         i030(:) = 6.0e-8*((300./t_elect(:))**0.64)
+         nb_reaction_4 = nb_reaction_4 + 1
+         v_4(:,nb_reaction_4) = i030(:)
 
       end if   !ionchem
@@ -2425 +2431 @@
 
    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 
+!===========================================================
+
+   nb_reaction_4 = nb_reaction_4 + 1
+
+   indice_4(nb_reaction_4) = z4spec(1.0, i_hco2plus, 1.0, i_elec, 1.0, i_h, 1.0, i_co2)
 
 end if    !ionchem