Changeset 2186 for trunk/LMDZ.VENUS

Timestamp:

Dec 6, 2019, 6:07:58 PM (5 years ago)

Author:

flefevre

Message:

Venus photochemistry: updated kinetics (Lisboa)

File:

• trunk/LMDZ.VENUS/libf/phyvenus/new_photochemistry_venus.F90 (modified) (11 diffs)

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

@@ -29 +29 @@
 
 integer :: ij, i, iesp, jesp, iz, iiz
+real :: fact_SF        !  Increasin/Decreasing Solar Flux
 
 ! matrix
@@ -436 +437 @@
          read(30,'(7e11.4)') (jphot(iso2,isza,iz,ij), ij = 1,nj)
          do ij = 1,nj
+            !!!    
+            ! GG TESTS increasing/decreasing Solar Flux 
+            !fact_SF = 5.
+            ! jphot(iso2,isza,iz,ij)=jphot(iso2,isza,iz,ij)*0.2
+
+            !   print*,'!!! WARNING !!!: You are Modifying Solar flux by a factor:'
+            !   print*, ' Please check new_chemistry* and sw_venus_rh ' 
+            !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
+
             if (jphot(iso2,isza,iz,ij) == 1.E-30) then 
                jphot(iso2,isza,iz,ij) = 0.
@@ -2738 +2748 @@
 
 !     jpl 2003
-
-      b003(:) = 1.1E-10
+!      b003(:) = 1.1E-10
+
+!      jpl 2016      
+      b003(:) = 1.2E-10
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -2787 +2799 @@
 
 !     jpl 2003
-
-      c002(:) = 2.2E-11*exp(120./t(:))
+!      c002(:) = 2.2E-11*exp(120./t(:))
+
+!     jpl 2016  
+      c002(:) = 1.8E-11*exp(180./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -2845 +2859 @@
 
 !     christensen et al., grl, 13, 2002
-
-      c008(:) = 1.5E-12*exp(19./t(:))
+!      c008(:) = 1.5E-12*exp(19./t(:))
+
+!     jpl 2016
+      c008(:) = 3.0E-13*exp(460./t(:))
 
       nb_reaction_3 = nb_reaction_3 + 1
@@ -2924 +2940 @@
 
 !     jpl 2003
-
-      c016(:) = 2.5*1.7E-33*exp(1000./t(:))*conc(:)
+!      c016(:) = 2.5*1.7E-33*exp(1000./t(:))*conc(:)
+
+!     jpl 2016
+      c016(:) = 2.5*2.1E-33*exp(920./t(:))*conc(:)
 
       nb_reaction_3 = nb_reaction_3 + 1
@@ -3072 +3090 @@
 
 !     jpl 2006
-
-      f004(:) = 2.6E-12*exp(-350./t(:))
+!      f004(:) = 2.6E-12*exp(-350./t(:))
+
+!     jpl 2016
+      f004(:) = 1.8E-12*exp(-250./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -3126 +3146 @@
 
 !     jpl 2006
-
-      f010(:) = 4.1E-11*exp(-450./t(:))
-
+!      f010(:) = 4.1E-11*exp(-450./t(:))
+
+!     jpl 2016
+      f010(:) = 3.6E-11*exp(-375./t(:)) 
+      
       nb_reaction_4 = nb_reaction_4 + 1
       v_4(:,nb_reaction_4) = f010(:)
@@ -3135 +3157 @@
 
 !     jpl 2006
-
-      f011(:) = 1.8E-11*exp(170./t(:))
+!      f011(:) = 1.8E-11*exp(170./t(:))
+
+!     jpl 2016
+      f011(:) = 1.4E-11*exp(270./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -3423 +3447 @@
 !---  g001: s + o2 -> so + o
 
-      g001(:) = 2.3E-12
+!      g001(:) = 2.3E-12
+
+!     jpl 2016
+      g001(:) = 1.6E-12*exp(100./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1
@@ -3437 +3464 @@
 !---  g003: so + o2 -> so2 + o
 
-      g003(:) = 1.25E-13*exp(-2190./t(:))
+!      g003(:) = 1.25E-13*exp(-2190./t(:))
+
+!     jpl 2016
+       g003(:) = 1.6E-13*exp(-2280./t(:))
 
       nb_reaction_4 = nb_reaction_4 + 1