Changeset 2851 for trunk/LMDZ.VENUS/libf/phyvenus/photolysis_online.F

Timestamp:

Dec 21, 2022, 10:55:05 AM (2 years ago)

Author:

streel

Message:

LMDZ.VENUS - chemical reaction rates update with JPL 2019

• Update of the main 1D subroutine with nitrogen species + generalisation of VMR output every 1/5 venusian day
• Added the H2 photolysis
• Added the traceur.def for all neutral species only
• Added the VMR def files for 1D model

NS

File:

• trunk/LMDZ.VENUS/libf/phyvenus/photolysis_online.F (modified) (7 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/photolysis_online.F

@@ -3 +3 @@
       subroutine photolysis_online(nlayer, nb_phot_max,
      $           alt, press, temp, zmmean, 
-     $           i_co2, i_co, i_o, i_o1d, i_o2, i_o3,
+     $           i_co2, i_co, i_o, i_o1d, i_o2, i_o3,i_h2,
      $           i_oh, i_ho2, i_h2o2, i_h2o,i_h,i_hcl,
      $           i_cl2, i_hocl, i_so2, i_so, i_so3,
@@ -23 +23 @@
      $                       i_cl2, i_hocl, i_so2, i_so, i_so3, i_clo,
      $                       i_ocs, i_cl, i_cocl2, i_h2so4, i_no2,
-     $                       i_no, i_n2, i_n2d 
+     $                       i_no, i_n2, i_n2d, i_h2 
 
       real, dimension(nlayer), intent(in) :: press, temp,  zmmean    ! pressure (hpa)/temperature (k)/mean molecular mass (g.mol-1)
@@ -67 +67 @@
      $           j_h2o, j_h2o2, j_ho2, j_h, j_hcl, j_cl2, j_hocl, j_so2,
      $           j_so, j_so3, j_clo, j_ocs, j_cocl2, j_h2so4, j_no2,
-     $           j_no, j_n2
+     $           j_no, j_n2, j_h2
 
       integer :: a_o2, a_co2, a_o3, a_h2o, a_h2o2, a_ho2, a_hcl, a_cl2,
      $           a_hocl, a_so2, a_so, a_so3, a_clo, a_ocs, a_cocl2,
-     $           a_h2so4, a_no2, a_no, a_n2
+     $           a_h2so4, a_no2, a_no, a_n2, a_h2
       integer :: i, ilay, iw, ialt
       real    :: deltaj
@@ -81 +81 @@
       a_co2     = 2      ! co2
       a_o3      = 3      ! o3
-      a_h2o     = 4      ! h2o
-      a_h2o2    = 5      ! h2o2
-      a_ho2     = 6      ! ho2
-      a_hcl     = 7      ! hcl
-      a_cl2     = 8      ! cl2
-      a_hocl    = 9      ! hocl
-      a_so2     = 10     ! so2
-      a_so      = 11     ! so
-      a_so3     = 12     ! so3
-      a_clo     = 13     ! clo
-      a_ocs     = 14     ! ocs
-      a_cocl2   = 15     ! cocl2
-      a_h2so4   = 16     ! h2so4
-      a_no2     = 17     ! no2
-      a_no      = 18     ! no
-      a_n2      = 19     ! n2
+      a_h2      = 4      ! h2
+      a_h2o     = 5      ! h2o
+      a_h2o2    = 6      ! h2o2
+      a_ho2     = 7      ! ho2
+      a_hcl     = 8      ! hcl
+      a_cl2     = 9      ! cl2
+      a_hocl    = 10     ! hocl
+      a_so2     = 11     ! so2
+      a_so      = 12     ! so
+      a_so3     = 13     ! so3
+      a_clo     = 14     ! clo
+      a_ocs     = 15     ! ocs
+      a_cocl2   = 16     ! cocl2
+      a_h2so4   = 17     ! h2so4
+      a_no2     = 18     ! no2
+      a_no      = 19     ! no
+      a_n2      = 20     ! n2
 
 !     photodissociation rates numbering. 
@@ -107 +108 @@
       j_o3_o1d  = 5      ! o3 + hv     -> o2 + o(1d)
       j_o3_o    = 6      ! o3 + hv     -> o2 + o
-      j_h2o     = 7      ! h2o + hv    -> h + oh
-      j_ho2     = 8      ! ho2 + hv    -> oh + o
-      j_h2o2    = 9      ! h2o2 + hv   -> oh + oh
-      j_hcl     = 10     ! hcl + hv    -> h + cl
-      j_cl2     = 11     ! cl2 + hv    -> cl + cl 
-      j_hocl    = 12     ! hocl + hv   -> oh + cl
-      j_so2     = 13     ! so2 + hv    -> so + o
-      j_so      = 14     ! so + hv     -> s + o
-      j_so3     = 15     ! so3 + hv    -> so2 + o
-      j_clo     = 16     ! clo + hv    -> cl + o
-      j_ocs     = 17     ! ocs + hv    -> co + s
-      j_cocl2   = 18     ! cocl2 + hv  -> 2cl + co 
-      j_h2so4   = 19     ! h2so4 + hv  -> so3 + h2o
-      j_no2     = 20     ! no2 + hv    -> no + o
-      j_no      = 21     ! no + hv     -> n + o
-      j_n2      = 22     ! n2 + hv     -> n(2d) + n
+      j_h2      = 7      ! h2 + hv     -> h + h
+      j_h2o     = 8      ! h2o + hv    -> h + oh
+      j_ho2     = 9      ! ho2 + hv    -> oh + o
+      j_h2o2    = 10     ! h2o2 + hv   -> oh + oh
+      j_hcl     = 11     ! hcl + hv    -> h + cl
+      j_cl2     = 12     ! cl2 + hv    -> cl + cl 
+      j_hocl    = 13     ! hocl + hv   -> oh + cl
+      j_so2     = 14     ! so2 + hv    -> so + o
+      j_so      = 15     ! so + hv     -> s + o
+      j_so3     = 16     ! so3 + hv    -> so2 + o
+      j_clo     = 17     ! clo + hv    -> cl + o
+      j_ocs     = 18     ! ocs + hv    -> co + s
+      j_cocl2   = 19     ! cocl2 + hv  -> 2cl + co 
+      j_h2so4   = 20     ! h2so4 + hv  -> so3 + h2o
+      j_no2     = 21     ! no2 + hv    -> no + o
+      j_no      = 22     ! no + hv     -> n + o
+      j_n2      = 23     ! n2 + hv     -> n(2d) + n
 
 !     j_hdo_od  =        ! hdo + hv    -> od + h
@@ -174 +176 @@
       do ilay = 1,nlayer
          do iw = 1,nw-1
+            dtgas(ilay,iw,a_h2) = colinc(ilay)*rm(ilay,i_h2)*xsh2(iw)
             dtgas(ilay,iw,a_h2o) = colinc(ilay)*rm(ilay,i_h2o)*xsh2o(iw)
             dtgas(ilay,iw,a_ho2) = colinc(ilay)*rm(ilay,i_ho2)*xsho2(iw)
@@ -209 +212 @@
       do ilay = 1,nlayer
          do iw = 1,nw-1
+            sj(ilay,iw,j_h2) = xsh2(iw)                  ! h2
             sj(ilay,iw,j_h2o) = xsh2o(iw)                ! h2o
             sj(ilay,iw,j_ho2) = xsho2(iw)                ! ho2