Changeset 2862 for trunk/LMDZ.VENUS

Timestamp:

Jan 4, 2023, 3:40:26 PM (23 months ago)

Author:

flefevre

Message:

Cosmetic changes and some more info on reaction SO2 + O + CO2 -> SO3 + CO2

File:

• trunk/LMDZ.VENUS/libf/phyvenus/photochemistry_venus.F90 (modified) (1 diff)

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

@@ -3630 +3630 @@
 !---  g009: so2 + o + co2 -> so3 + co2
 
-!     Naido 2005
-
-     do iz = 1,nz
-        ak0 = 5.*9.5*1.E-23*(t(iz)**(-3.0))*EXP(-2400./t(iz))
-        ak1 = 6.1*1.E-13*EXP(-850./t(iz))
-        rate = (ak0*conc(iz))/(1. + ak0*conc(iz)/ak1)
-        xpo = 1./(1. + alog10((ak0*conc(iz))/ak1)**2)
-        fc = 0.558*EXP(-t(iz)/316.)+0.442*EXP(-t(iz)/7442.)
-        g009(iz) = rate*fc**xpo
-     end do
-
-!     jpl 2019
-
-!      do iz = 1,nz
-!         ak0 = 2.5*1.8e-33*(t(iz)/298)**(2.0)
-!         ak1 = 4.1e-14*(t(iz)/298)**(1.8)
-
-!         rate = (ak0*conc(iz))/(1. + ak0*conc(iz)/ak1)
-!         xpo = 1./(1. + log10((ak0*conc(iz))/ak1)**2)
-!         g009(iz) = rate*0.6**xpo
-!      end do
+!     Naidoo et al., Proceedings of the Combustion Institute,
+!                    30, 1219-1225, 2005
+!     Also recommended by jpl 2019 with a simpler expression (fc = 0.6)
+!
+!     Factor of 5 for third-body efficiency of co2 to be confirmed!
+
+      do iz = 1,nz
+         ak0 = 5.*9.5e-23*(t(iz)**(-3.0))*exp(-2400./t(iz))
+         ak1 = 6.1e-13*exp(-850./t(iz))
+         rate = (ak0*conc(iz))/(1. + ak0*conc(iz)/ak1)
+         xpo = 1./(1. + alog10((ak0*conc(iz))/ak1)**2)
+         fc = 0.558*exp(-t(iz)/316.) + 0.442*exp(-t(iz)/7442.)
+         g009(iz) = rate*fc**xpo
+      end do
 
       nb_reaction_4 = nb_reaction_4 + 1