Changeset 3135

Timestamp:

Nov 21, 2023, 4:07:58 PM (19 months ago)

Author:

flefevre

Message:

• mise a jour des reactions O + O2 + M et O + O + M
• un peu de cosmetique

File:

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

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

@@ -17 +17 @@
 !*****************************************************************
 
-subroutine photochemistry_venus(nz, n_lon, zlocal, ptimestep,                 &
-                                    ok_jonline,ok_ionchem,tuneupperatm,       &
-                                    nb_reaction_3_max,nb_reaction_4_max,      &
-                                    nb_phot_max,nphotion,ig,                  &
-                                    p, t, t_elect, tr, vmr_dens_euv, mumean,  &
-                                    sza_input, lon, lat, nesp, nespeuv, iter, &
-                                    prod_tr, loss_tr, em_no, em_o2)  
+subroutine photochemistry_venus(nz, n_lon, zlocal, ptimestep,             &
+                                ok_jonline, ok_ionchem, tuneupperatm,     &
+                                nb_reaction_3_max, nb_reaction_4_max,     &
+                                nb_phot_max, nphotion, ig,                &
+                                p, t, t_elect, tr, vmr_dens_euv, mumean,  &
+                                sza_input, lon, lat, nesp, nespeuv, iter, &
+                                prod_tr, loss_tr, em_no, em_o2)  
 
 use chemparam_mod
@@ -137 +137 @@
 
 ! reaction rates
-
-!integer, parameter :: nb_phot_max       = 35
-!integer, parameter :: nb_phot           = 22
-!integer, parameter :: nb_reaction_3_max = 12
-!integer, parameter :: nb_reaction_4_max = 98
 
 real (kind = 8), dimension(nz,      nb_phot_max) :: v_phot
@@ -2581 +2576 @@
 
 integer :: iz
-real    :: ak0, ak1, xpo, rate, rate1, rate2, pi, gam
+real    :: ak0, ak1, xpo, rate, rate1, rate2, pi, gam, epsil
 real    :: k1a0, k1b0, k1ainf, k1a, k1b, k0, kinf, kf, kint, kca, fc, fx, x, y
 integer :: nb_phot, nb_reaction_3, nb_reaction_4
@@ -2638 +2633 @@
 !----------------------------------------------------------------------
 
-!     jpl 2003
-!     ------ BEFORE VCD 2.1 -------
-!      a001(:) = 2.5*6.0E-34*(t(:)/300.)**(-2.4)*conc(:)
-!     -----------------------------
-
-!    ! nist expression + take account of the N2 and CO2 factor
-!     ------ VCD 2.1 update -------
-      a001(:) = 6.0E-34*(t(:)/300.)**(-2.4) *    &
-              ( 2.5 * c(:,i_co2) +               &
-                1.0 * (conc(:)-c(:,i_co2))  ) 
-              nb_reaction_4 = nb_reaction_4 + 1
-              v_4(:,nb_reaction_4) = a001(:)
-
-
-!     Tsang and Hampson, J. Chem. Phys. Ref. Data, 15, 1087, 1986
-!     ------ BEFORE VCD 2.1 -------
-!      a002(:) = 2.5*9.46E-34*exp(485./t(:))*conc(:) ! nist expression
-!     -----------------------------
-
-!     ! nist expression + take account of the N2 and CO2 factor
-!       a002(:) = 9.46E-34*exp(485./t(:)) *     &
-!              ( 2.5 * c(:,i_co2) +            &
-!                1.0 * (conc(:)-c(:,i_co2)) )
-
-!     Baulch et al., 1976 (confirmed by smith and robertson, 2008)
-!     ------ VCD 2.1 update -------
-
-       a002(:) = 2.76E-34*exp(720./t(:)) *    &
-              ( 2.5 * c(:,i_co2) +            & 
-                1.0 * (conc(:)-c(:,i_co2)) )
-!     -----------------------------
+!---  a001: o + o2 + m -> o3 + m
+
+!     jpl 2019
+ 
+!     co2/n2 efficiency as a third body = 2.075
+!     from sehested et al., j. geophys. res., 100, 1995.
+
+      a001(:) = 6.1e-34*(t(:)/298.)**(-2.4)     &
+              *(2.075*c(:,i_co2) + 1.0*(conc(:) - c(:,i_co2))) 
+
+      nb_reaction_4 = nb_reaction_4 + 1
+      v_4(:,nb_reaction_4) = a001(:)
+
+!---  a002: o + o + m -> o2(delta_g) + m
+ 
+!     baulch et al., 1976 (confirmed by smith and robertson, 2008)
+
+!     epsil : net effective yield 
+
+      epsil = 0.75 ! (crisp et al., 1996; krasnopolsky, 1991)
+
+      a002(:) = 2.76e-34*exp(720./t(:))  &
+              *(2.5*c(:,i_co2) + 1.0*(conc(:) - c(:,i_co2)))*epsil
       
       nb_reaction_3 = nb_reaction_3 + 1
       v_3(:,nb_reaction_3) = a002(:)
+
       ind_orec = nb_reaction_3
 
-!     a003: o + o3 -> o2 + o2
+!---  a003: o + o3 -> o2 + o2
 
 !     jpl 2003