Changeset 3764 for trunk/LMDZ.VENUS

Timestamp:

May 15, 2025, 2:01:57 PM (6 weeks ago)

Author:

flefevre

Message:

Implementation of ClSO2 + Cl -> Cl2SO2

Location:

trunk/LMDZ.VENUS/libf/phyvenus

Files:

• photochemistry_venus.F90 (modified) (4 diffs)
• phytrac_chimie.F (modified) (1 diff)

trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90 ¶

@@ -1517 +1517 @@
 
 !===========================================================
+!      f040 : ClSO2 + Cl -> Cl2SO2
+!===========================================================
+
+nb_reaction_4 = nb_reaction_4 + 1
+
+indice_4(nb_reaction_4) = z4spec(1.0, i_clso2, 1.0, i_cl, 1.0, i_cl2so2, 0.0, i_dummy)
+
+!===========================================================
 !      g001 : S + O2 -> SO + O
 !===========================================================
@@ -2673 +2681 @@
                        f022, f023, f024, f025, f026, f027, f028,   &
                        f029, f030, f031, f032, f033, f034, f035,   &
-                       f036, f037, f038, f039,                     &
+                       f036, f037, f038, f039, f040,               &
                        g001, g002, g003, g004, g005, g006, g007,   &
                        g008, g009, g010, g011, g012, g013, g014,   &
@@ -3566 +3574 @@
       end do
 
+      f030(:) = 0.  ! temporary FL 15 May 2025
+
       nb_reaction_3 = nb_reaction_3 + 1
       v_3(:,nb_reaction_3) = f030(:)
@@ -3649 +3659 @@
       nb_reaction_4 = nb_reaction_4 + 1
       v_4(:,nb_reaction_4) = f039(:)
+
+!---  f040: clso2 + cl -> cl2so2
+
+!     croce and cobos, 2018
+
+      f040(:) = 1.44e-11*(t(:)/250.)**0.47
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = f040(:)
 
 !----------------------------------------------------------------------

trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F ¶

@@ -143 +143 @@
            ! default numbers if no ion and no deuterium chemistry
     
-           nb_reaction_4_max = 98     ! set number of bimolecular reactions
+           nb_reaction_4_max = 99     ! set number of bimolecular reactions
            nb_reaction_3_max = 14     ! set number of quadratic reactions
            nquench           = 15     ! set number of first-order reactions: