Changeset 2843 for trunk/LMDZ.VENUS

Timestamp:

Dec 16, 2022, 2:46:20 PM (2 years ago)

Author:

flefevre

Message:

Venus chemical kinetics: cosmetic changes prior to JPL 2019 update

File:

• trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90 (modified) (13 diffs)

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

@@ -2700 +2700 @@
 !---  b003: o(1d) + h2  -> oh + h
 
-!     jpl 2003
-!      b003(:) = 1.1E-10
-
-!      jpl 2016      
+!     jpl 2011      
+
       b003(:) = 1.2E-10
 
@@ -2751 +2749 @@
 !---  c002: o + oh -> o2 + h
 
+!     jpl 2011  
+
+      c002(:) = 1.8E-11*exp(180./t(:))
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = c002(:)
+
+!---  c003: h + o3 -> oh + o2
+
 !     jpl 2003
-!      c002(:) = 2.2E-11*exp(120./t(:))
-
-!     jpl 2016  
-      c002(:) = 1.8E-11*exp(180./t(:))
-
-      nb_reaction_4 = nb_reaction_4 + 1
-      v_4(:,nb_reaction_4) = c002(:)
-
-!---  c003: h + o3 -> oh + o2
+
+      c003(:) = 1.4E-10*exp(-470./t(:))
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = c003(:)
+
+!---  c004: h + ho2 -> oh + oh
+
+!     jpl 2006
+
+      c004(:) = 7.2E-11
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = c004(:)
+
+!---  c005: h + ho2 -> h2 + o2
+
+!     jpl 2006
+
+      c005(:) = 6.9E-12
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = c005(:)
+
+!---  c006: h + ho2 -> h2o + o
+
+!     jpl 2006
+
+      c006(:) = 1.6E-12
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = c006(:)
+
+!---  c007: oh + ho2 -> h2o + o2
 
 !     jpl 2003
 
-      c003(:) = 1.4E-10*exp(-470./t(:))
-
-      nb_reaction_4 = nb_reaction_4 + 1
-      v_4(:,nb_reaction_4) = c003(:)
-
-!---  c004: h + ho2 -> oh + oh
-
-!     jpl 2006
-
-      c004(:) = 7.2E-11
-
-      nb_reaction_4 = nb_reaction_4 + 1
-      v_4(:,nb_reaction_4) = c004(:)
-
-!---  c005: h + ho2 -> h2 + o2
-
-!     jpl 2006
-
-      c005(:) = 6.9E-12
-
-      nb_reaction_4 = nb_reaction_4 + 1
-      v_4(:,nb_reaction_4) = c005(:)
-
-!---  c006: h + ho2 -> h2o + o
-
-!     jpl 2006
-
-      c006(:) = 1.6E-12
-
-      nb_reaction_4 = nb_reaction_4 + 1
-      v_4(:,nb_reaction_4) = c006(:)
-
-!---  c007: oh + ho2 -> h2o + o2
-
-!     jpl 2003
-
       c007(:) = 4.8E-11*exp(250./t(:))
 
@@ -2807 +2803 @@
 !---  c008: ho2 + ho2 -> h2o2 + o2
 
-!     jpl 2006
-
-!     c008(:) = 3.5E-13*exp(430./t(:))
+!     jpl 2015
+
+      c008(:) = 3.0E-13*exp(460./t(:))
 
 !     christensen et al., grl, 13, 2002
-!      c008(:) = 1.5E-12*exp(19./t(:))
-
-!     jpl 2016
-      c008(:) = 3.0E-13*exp(460./t(:))
+
+!     c008(:) = 1.5E-12*exp(19./t(:))
 
       nb_reaction_3 = nb_reaction_3 + 1
@@ -2892 +2886 @@
 !---  c016: ho2 + ho2 + co2 -> h2o2 + o2 + co2
 
-!     jpl 2003
-!      c016(:) = 2.5*1.7E-33*exp(1000./t(:))*conc(:)
-
-!     jpl 2016
+!     jpl 2011
+
       c016(:) = 2.5*2.1E-33*exp(920./t(:))*conc(:)
 
@@ -3135 +3127 @@
 !---  f004: hcl + oh -> h2o + cl
 
+!     jpl 2009
+
+      f004(:) = 1.8E-12*exp(-250./t(:))
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = f004(:)
+
+!---  f005: clo + o -> cl + o2
+
 !     jpl 2006
-!      f004(:) = 2.6E-12*exp(-350./t(:))
-
-!     jpl 2016
-      f004(:) = 1.8E-12*exp(-250./t(:))
-
-      nb_reaction_4 = nb_reaction_4 + 1
-      v_4(:,nb_reaction_4) = f004(:)
-
-!---  f005: clo + o -> cl + o2
+
+      f005(:) = 2.8E-11*exp(85./t(:))
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = f005(:)
+
+!---  f006: clo + oh -> cl + ho2
 
 !     jpl 2006
 
-      f005(:) = 2.8E-11*exp(85./t(:))
-
-      nb_reaction_4 = nb_reaction_4 + 1
-      v_4(:,nb_reaction_4) = f005(:)
-
-!---  f006: clo + oh -> cl + ho2
+      f006(:) = 7.4E-12*exp(270./t(:))
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = f006(:)
+
+!---  f007: clo + oh -> hcl + o2
 
 !     jpl 2006
 
-      f006(:) = 7.4E-12*exp(270./t(:))
-
-      nb_reaction_4 = nb_reaction_4 + 1
-      v_4(:,nb_reaction_4) = f006(:)
-
-!---  f007: clo + oh -> hcl + o2
+      f007(:) = 6.0E-13*exp(230./t(:))
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = f007(:)
+
+!---  f008: cl + h2 -> hcl + h
 
 !     jpl 2006
 
-      f007(:) = 6.0E-13*exp(230./t(:))
-
-      nb_reaction_4 = nb_reaction_4 + 1
-      v_4(:,nb_reaction_4) = f007(:)
-
-!---  f008: cl + h2 -> hcl + h
+      f008(:) = 3.05E-11*exp(-2270./t(:))
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = f008(:)
+
+!---  f009: cl + o3 -> clo + o2
 
 !     jpl 2006
 
-      f008(:) = 3.05E-11*exp(-2270./t(:))
-
-      nb_reaction_4 = nb_reaction_4 + 1
-      v_4(:,nb_reaction_4) = f008(:)
-
-!---  f009: cl + o3 -> clo + o2
-
-!     jpl 2006
-
       f009(:) = 2.3E-11*exp(-200./t(:))
 
@@ -3191 +3181 @@
 !---  f010: cl + ho2 -> clo + oh
 
-!     jpl 2006
-!      f010(:) = 4.1E-11*exp(-450./t(:))
-
-!     jpl 2016
+!     jpl 2009
+
       f010(:) = 3.6E-11*exp(-375./t(:)) 
       
@@ -3202 +3190 @@
 !---  f011: cl + ho2 -> hcl + o2
 
-!     jpl 2006
-!      f011(:) = 1.8E-11*exp(170./t(:))
-
-!     jpl 2016
+!     jpl 2009
+
       f011(:) = 1.4E-11*exp(270./t(:))
 
@@ -3493 +3479 @@
 !---  g001: s + o2 -> so + o
 
-!      g001(:) = 2.3E-12
-
-!     jpl 2016
+!     jpl 2015
+
       g001(:) = 1.6E-12*exp(100./t(:))
 
@@ -3503 +3488 @@
 !---  g002: s + o3 -> so + o2
 
+!     jpl 2015
+
       g002(:) = 1.2E-11
 
@@ -3510 +3497 @@
 !---  g003: so + o2 -> so2 + o
 
-!      g003(:) = 1.25E-13*exp(-2190./t(:))
-
-!     jpl 2016
-       g003(:) = 1.6E-13*exp(-2280./t(:))
+!     jpl 2015
+
+      g003(:) = 1.6E-13*exp(-2280./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -3520 +3506 @@
 !---  g004: so + o3 -> so2 + o2
 
+!     jpl 2015
+
       g004(:) = 3.4E-12*exp(-1100./t(:))
 
@@ -3527 +3515 @@
 !---  g005: so + oh -> so2 + h
 
+!     jpl 2015
+
       g005(:) = 2.7E-11*exp(335./t(:))
 
@@ -3533 +3523 @@
 
 !---  g006: s + oh -> so + h
+
+!     jpl 2015
 
       g006(:) = 6.6E-11